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Abstract 

The quantum well solar cell (QWSC) is a p - i - n solar cell with quantum 
wells in the intrinsic region. Previous work has shown that QWSCs have a greater 
open circuit voltage (Kc) than would be provided by a cell with the quantum 
well effective bandgap. This suggests that the fundamental efficiency limits of 
QWSCs are greater than those of single bandgap solar cells. 

The following work investigates QWSCs in the GaAs/AlxGai_xAs materials 
system. The design and optimisation of a QWSC in this system requires studies 
of the voltage and current dependencies on the aluminium fraction. QWSCs with 
different aluminium fractions have been studied and show an increasing Vgc with 
increasing barrier aluminium composition. 

The QE however decreases with increasing aluminium composition. We de- 
velop a model of the QE to test novel QWSC designs with a view to minimising 
this problem. 

This work concentrates on two design changes. The first deals with com- 
positionally graded structures in which the bandgap varies with position. This 
bandgap variation introduces an quasi electric field which can be used to increase 
minority carrier collection in the low efficiency p and n layers. This technique 
also increases the hght flux reaching the highly efficient depletion regions. 

The second design change consists of coating the back of the cell with a mirror 
to exploit the portion of light which is not absorbed on the first pass. 

A model of the QE of compositionally graded QWSC solar cells with back 
surface mirrors is developed in order to analyse the effect of these design changes. 
These changes are implemented separately in a number of QWSC designs and 
the resulting experimental data compared with the model. An optimised design 
is then presented. 
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Chapter 1 



Introduction 

1.1 Current Energy Use and Alternatives 

The sun irradiates the planet with approximately 4x10^ EJ per year p]. When we 
compare this enormous energy flux with current human consumption of around 
400 EJ per year, we find that we receive somewhat more than this amount every 
hour in direct solar radiation. 

With the exception of such sources as tidal, geothermal and nuclear energy, all 
the energy we consume is ultimately derived from the sun. It is largely harvested 
in a roundabout manner. The most important is the burning of fossil fuels, 
followed by nuclear and hydroelectric power. 

The limited nature of the fossil fuel reserves however makes the search for 
alternatives important. This is exacerbated by recent concerns about global 
warming through the continued emission of greenhouse gasses. The recently 
published results of the United Nations Advisory Group on Climate Change 
included the conclusions of a wide range of expert opinion in this field. The 
conclusion was that some level of global warming is expected, and indeed, is 
already measurable. 

Currently exploited alternatives present problems. Nuclear power, once hailed 
as the power source which would be too cheap to meter, has experienced well 
publicised difficulties related to the disposal of nuclear waste. 
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1.2. Reducing the Cost of Photovoltaic Power 
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Hydroelectric power has been heavily exploited for decades, and in many coun- 
tries is close to capacity. Moreover, the environmental impact of such schemes 
has, with hindsight, been found to be both significant and not always beneficial. 

Solar or photovoltaic (PV) power is an attractive alternative for a number of 
reasons. PV installations can be designed on a small scale, point-of-use basis. 
Such installations have a low land requirement and little environmental impact. 
They do not require mechanical moving parts, and are therefore potentially ro- 
bust and reliable, and subject to low maintainance costs. Photovoltaic module 
manufacturers already quote guaranteed lifetimes measured in decades. This 
compares favourably with diesel generators, for example, or even other renew- 
able energy sources such thermal or wind schemes which require moving parts 
and regular servicing. 

Despite all these advantages, we see that photovoltaic power has yet to make 
a significant impact worldwide. The high capital cost associated with the pro- 
duction of the modules is a significant problem. This is partly due to inertia in 
the electricity generating industries and to institutional barriers resulting from 
the interests of established industries. Lowering the cost of alternative energy 
sources however clearly remains a viable route to achieving sustainable energy 
consumption patterns. 

1.2 Reducing the Cost of Photovoltaic Power 

There are three ways to reduce the cost of photovoltaic power. One is to con- 
centrate on designing low cost PV modules, at some expense to the conversion 
efficiency. Larger areas will therefore be necessary, but the reduction in final cost 
per kilowatt may outweigh other considerations if the module cost can be brought 
down far enough. Examples are amorphous or poly-crystalline solar cells, and 
thin film solar cells made from II- VI materials such as CdTe. This strategy is 
currently the most successful for terrestrial PV applications. Another example 
of this technique is the electrochemical solar cell which has been reported by 
O'Regan and Gratzel [3]. 
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The second, related method is to reduce the cost of PV systems. The cost 
of the solar cells in such a module may constitute less than half the cost of the 
finished PV module. For example the multicrystalline solar cells in Photowatt's 
PV modules [2] comprise about 40% of the cost of the finished module. 

Another approach is improve the efficiency of the solar cells dramatically at 
the expense of cost. Again, the final module cost may be attractive in dollars per 
watt terms if the increase in output power more than compensates the increased 
module cost. In some cases, it may be economically viable to increase the light 
input to individual solar cells in a module by means of light concentrators. This 
approach typically involves exotic structures, complex processing, and generally 
requires high quality material. 

To date the low cost approach has delivered marketable PV modules, whereas 
high efficiency cells are reserved for specialised applications in which cost is a 
secondary issue, such as space applications. 

1.3 Introduction to the Optimisation Program 

This work concerns the optimisation of a novel design of solar cell with the po- 
tential to reach high efficiencies through the use of novel physical processes. This 
design is the quantum well solar cell (QSWC). This is studied here particularly 
with regard to the efficiency with which the solar cell converts light into current, 
which is called the quantum efficiency or QE. 

Chapter |2] will give an introduction underlying semiconductor physics before 
defining important solar cell parameters. A short review of some standard and 
high efficiency solar cell designs will be investigated, including the fundamental 
principles behind the quantum well solar cell. 

Chapter [3] outlines the theoretical basis for QE calculations in general and 
describes the model we will use. This model requires knowledge of a large number 
of materials parameters, which are reviewed in chapter |H 

Chapter |5] describes previous work on the optimisation of the QWSC struc- 
ture, and illustrates additional theoretical results from the model for a number 
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of designs based on the QWSC A number of optimisations are presented. 

The practical aspects of growth together with a review of experimental char- 
acterisation techniques are reviewed in chapter O These include preliminary 
characterisation techniques which allow the quality of individual growth and pro- 
cessing runs to be assessed. Results for these preliminary techniques are given in 
chapter [71 

Chapter [8] presents experimental measurements of the QE for a wide range 
of samples grown by two different techniques. These samples cover the optimi- 
sations suggested by the theory. The experimental data are then compared with 
the modelling, and conclusions drawn regarding the viability of the techniques 
proposed in chapter [5l 



Chapter 2 

Introduction to Solar Cells 



In this chapter we introduce parameters such as short circuit current, open cir- 
cuit voltage and fill factor, and overall conversion efficiency, which are useful for 
characterising and optimising solar cells. The first sections summarise the basic 
semiconductor physics of semiconductor p - n junctions in equilibrium and under 
illumination. These allows us to derive an expression for the built in potential of 
a p - n junction. 

The behaviour of a p - n junction under different applied bias conditions is 
described as a basis for understanding the current voltage characteristic of a solar 
cell in the dark and under illumination. We describe the principles underlying 
the photovoltaic effect, and can then define the solar cell parameters mentioned 
above. 

Having defined these parameters, we can identify design requirements for 
efficient solar cells. We review a number of high efficiency solar cell designs in 

a range of materials. We introduce the quantum well solar cell and explain why 
this design is expected to yield higher fundamental efficiencies limits than are 
expected from single band-gap solar cells. 

2.1 Physics of Solar Cells 

In this section we introduce the fundamental physics of semiconductors required 
to understand and optimise solar cell designs, and to define all fundamental 
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physical concepts and parameters which will be used in subsequent sections. The 
discussion is mainly based in the treatment given in Sze [1]. 

The crystal potential of a crystalline solid restricts the energies which can be 
occupied by electrons into bands of electronic states. The highest occupied band 
in energy in a semiconductor is the valence band, and is completely filled by the 
outer valence electrons of the solid at zero Kelvin. 

The band immediately above it is called the conduction band, and is empty 
at zero Kelvin. These two bands are separated by a gap in energy where no 
electronic states exist. This band-gap is principal bandgap of the semiconductor 
Eg, and has important consequences for the electrical and optical properties of 
these materials. 

2.1.1 Intrinsic Semiconductors in Equilibrium 

The probability that an electron in a solid should occupy a given energy level is 
described by the Fermi-Dirac distribution function f{E) of the form 

= 1 + ,(E-E,)A.T (2-1) 

where Ep is the Fermi level. Equation 12.11 defines it as the energy at which 
there is a 50% chance of occupation by an electron at any temperature, fcs is 
Boltzmann's constant and T is the temperature. 

At room temperatures thermal excitation in an undoped semiconductor gen- 
erates a finite equilibrium density of electron-hole pairs known as the intrinsic 
carrier concentration n^. This quantity influences the properties of semiconductor 
devices, and solar cells in particular. 

The densities of holes Pq and of electrons no are equal in an intrinsic semicon- 
ductor because each electron leaves exactly one hole behind in the valence band. 
They are both equal to n,. 

no = po = rii (2.2) 
The semiconductor in this case can conduct current through the motion under 
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an applied bias of the electrons in the conduction band and of hole states in the 
valence band. 

We now consider the magnitude of this intrinsic carrier concentration by de- 
termining the number of electrons excited into the valence band at a temperature 
T. This can be found by integrating the product of the occupation probability 
f{E) and the number of available states N{E). We consider only direct band 
gap semiconductors for simplicity. 

In this case, the number of available states near the band-edges may be as- 
sumed to be parabolic. The relationship between energy and electron momentum 
is analogous to the free-electron case, with the electron mass mg replaced by an 
effective mass m*, to account for the crystal potential. It can be shown [4J that 
the energy-momentum relationship yields the following density of states near the 
conduction band edge: 

NAE) = '-^^^iE-E^y'^ (2.3) 

A similar relationship for the density of states near the valence band edge can 
be found by substituting the hole effective mass mass for m* and integrating 
in energy from the valence band edge to minus infinity: from the valence band 
edge to infinity 

N.iE) = ^'^lf'^\ K-EY/^ (2.4) 
For energies E — Ep > 3KbT, the Fermi-Dirac function can be approximated 

by 

fE ~ e-^^-^^)/'^^^ (2.5) 

Integrating the product of equations l2.3l and l2.5l over energies from the conduction 
band edge to infinity then gives the number of conduction electrons in the system: 



no 



where Nc is known as the effective density of states in the conduction band and 
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is given by 

The number of holes is found by assuming a similar form to equation 12.51 for 
Ep — E > 3KbT and integrating downwards in energy from the valence band 
edge. This gives 

where Ny is the valence band effective density of states, given by 



'2'KmlkBT' 



3/2 



Ny = 2 y j (2.9) 

Equations 12.61 and 12.81 together with equation 12.21 yield the intrinsic carrier 
concentration 

nl = Pono = NcNye-^'o/''^^ (2.10) 

The intrinsic carrier concentration in GaAs for example is approximately 1.79 x 10^ 
cm~^ and a rather higher 1.45 x 10^°cm~^ in Si at room temperature because of 
the smaller band gap of this material. 

We can now calculate the Fermi level in an intrinsic semiconductor. Equation 
12.21 shows that the ratio of equations 12.61 and 12.81 is equal to one Separating Ep 
gives 

The Fermi level is situated near the middle of the bandgap, with a small offset 
caused by the difference between the effective densities of states in the valence 
and conduction bands. 



2.1.2 Doped Semiconductors 

We now consider the case of a semiconductor which has been doped with impu- 
rity atoms, starting with the case of a n type semiconductor. Donor or n type 
impurities have valence electrons with energy levels located a few meV below the 
conduction band edge. At non-zero temperatures, some of these donor atoms 
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are ionised by lattice vibrations and contribute extra electrons to the conduction 
band edge. Charge neutrality is preserved, since an ionised donor exists for every 
extra conduction band electron. For most common semiconductor devices, the 
dopants are chosen so that their ionisation energy is of the order or less that 
SksT, where T is the operating temperature. In this case, the donors are effec- 
tively all ionised, and the extra conduction electron concentration is of the order 
of the doping density A^'^. 

The law of mass action states that the electron and hole concentrations p and 
n are related as follows: 



Useful doping densities in GaAs/AlxGai_xAs are typically in the range 10^^ - 
10^'' cm~^. We have seen previously that the is of the order of lO^cm"^ for 
GaAs at room temperature. Since n ~ Nd, it can be seen from equation 12.121 
that the hole concentration is typically many orders of magnitude smaller than 
the electron concentration, and is given approximately by 



The electrons in n-type material are called majority carriers, and the holes minor- 
ity carriers. We use the subscript nO to identify the minority hole concentration 
in n-type material in equilibrium. 

A p-type semiconductor is one which has been doped with electron acceptors 
with empty electron valence states a few meV above the valence band. Electrons 
from the valence band are thermally excited into these states at room temper- 
atures, ionising the acceptor and leaving a positive hole in the valence band. 
The majority carriers in this case are holes, and their concentration is approxi- 
mately equal to the concentration Na of acceptor atoms. The minority electron 
concentration is then given by 



where the subscript pO identifies the electron minority concentration in p-type 




(2.12) 




(2.13) 




(2.14) 
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material in equilibrium. 

We can calculate the Fermi level in a n doped semiconductor from equation 
12.6^ noting that n ~ N£,. This gives 



Equation 12 . 1 5 1 shows that as approaches the effective conduction band density 
of states, the Fermi level increases towards the conduction band edge. This 
reflects the increasing probability of electronic states in the conduction band 
being occupied by electrons. 

Similarly, the Fermi level in a p type semiconductor is 



The Fermi level for a p type semiconductor moves towards the valence band 
edge as electrons are removed by acceptors, thereby decreasing the probability of 
valence band states being occupied by electrons. 



In the preceding sections, we looked at the static behaviour of electron and hole 
populations for semiconductors in equilibrium. We now examine the behaviour 
of these populations when the system is moved away from equilibrium under 
illumination, when electrical power can be extracted via the photovoltaic effect. 

The photovoltaic effect converts incident light into electrical power. The 
process consists of two stages. The first is to create an excess minority carrier 
population by exciting electrons with a light flux of intensity -F(A), in photons per 
second per wavelength A. By excess minority carrier population we mean those 
minority carriers which exceed the minority carrier population in equilibrium. 
If the incident light has sufficient energy, an electron can absorb a photon and 
transfer from a valence band state to a conduction band state. The net effect 
is the creation of an electron-hole pair. The excess minority carriers form a non 




(2.15) 




(2.16) 
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equilibrium population at a higher energy than the ground state of the system. 

The second stage of the photovoltaic process is to extract the potential energy 
of this system. It is known as charge separation, and consists of spatially sepa- 
rating the electron-hole pairs. This prevents them from recombining and losing 
their energy as heat or light, and imposes motion in opposite directions which 
constitues the photocurrent. This photocurrent can then be extracted through 
an external circuit connected to the terminals of the solar cell. 

2.2.1 Light Absorption and Pair Generation 

Minority carriers are generated when an electron absorbs a photon and gains suf- 
ficient energy to be promoted to a conduction state, leaving a positively charged 
hole state in the valence band. This is illustrated in figure 12.11 

In p type material, electrons which absorb a photon with energy greater than 
the bandgap are promoted into minority carrier states high in the conduction 
band. They very rapidly lose this excess energy in collisions with the lattice, 
generating heat. This process, known as thermalisation, ceases when they reach 
the band edge. The minority carriers then remain in band edge states for a 
finite time known as the minority electron lifetime r„ before recombining with a 
majority hole in the valence band. The time taken to thermalise to the band-edge 
is much smaller than r„, and the carriers therefore build up a quasi-equilibrium 
thermal population at the conduction band edge. 

Minority holes in n type material which absorb more energy than the bandgap 
are promoted to holes in the valence band, and thermalise up towards the valence 
band edge in a similar manner, again losing the excess energy as heat. They 
then remain at the band edge for a minority hole recombination time Tp before 
recombining with a majority electron. 

Thermalisation is an important energy loss mechanism in solar cells. It en- 
sures that the maximum amount of power delivered by a minority carrier is 
proportional to the bandgap of the material, irrespective of the energy of the 
photon providing the excitation. 
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Figure 2.1: Minority carrier photo generation, thermalisation and radiative and non radiative 
recombination mechanisms 



The generation rate G at a particular point is given by the absorption coeffi- 
cient a and the hght flux F 

G = aF (2.17) 

Direct and indirect semiconductors show distinctly different optical and elec- 
trical behaviour. In direct material such as GaAs the lowest electronic transition 
is one which does not require a change of momentum. In this case the absorption 
as a function of photon energy E near the band edge is approximately given (H] 

by 

ad^rect{E) ^ {E - Egf'^f{E - Eg) (2.18) 

where f{E — Eg) is a step function with value 1 for E > Eg and zero for E < Eg. 

Indirect material such as silicon or AlxGai_.xAs with an aluminium content 
greater than about 45% has an indirect bandgap. Light absorption also requires 
phonon absorption or emission. As a result, the absorption coefficient near the 
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Figure 2.2: p and n type semiconductors in equilibrium a) before contact and b) after contact. 
The built in voltage Vu is equal to the difference between the Fermi levels in a), before contact. 
Also shown are ionised donor (+) and acceptor (-) atoms which give rise to Hi- 

band-edge is significantly smaller than in direct material and takes the form 



The electron hole pairs generated by light absorption are free to move through 
the crystal lattice until they either recombine, losing the absorbed energy as heat 
or light, or are collected and contribute to the photocurrent. 

2.2.2 Charge Separation 

Collection of photogenerated carriers requires a field to separate electrons from 
holes, and to impose ordered motion in the direction of the positive and negative 
terminals of the solar cell. This constitutes the photocurrent. 

The charge separation is achieved by creating a junction region in the solar 
cell with p type doped material on one side, and n type dopants on the other. 
This results in a potential drop across the structure. Electrons drift down the 
potential, whereas holes drift up in the opposite direction. The precise origin of 
this separation potential is explained in more detail below. 

The solar cell is divided into three regions. Charge separation occurs mainly 




(2.19) 
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in the junction region, which we will refer to as the space-charge region. In good 
material, the field separates the electrons and holes with efficiencies approaching 
100%. 

There is no net electric field in the charge neutral sections of the p and n 
regions adjoining the depletion layer, which we call the p and n neutral regions. 

Electrons and holes in these neutral regions move by diffusion. Spatial vari- 
ations in composition of doping may contribute a small drift component. This 
will be investigated in chapter |3l The diffusion process is much less efficient than 
the drift collection which dominates in the space-charge region. Minority carri- 
ers perform a random walk with a finite average lifetime and diffusion length. 
A minority carrier which does not reach the interface with the space-charge re- 
gion eventually recombines with a majority carrier and is lost. The collection of 
carriers generated in the neutral regions is described in more detail in chapter [3J 

Current generation in all regions of the cell therefore depends on the presence 
of a built-in potential, which arises when a junction is made between p and n 
type material. 

When p and n type semiconductors are placed in electrical contact, there is 
an imbalance between the large number of majority electrons on the n side and 
the thermal minority electron concentration on the p side. A similar but opposite 
concentration gradient exists for holes. These electron and hole gradients give 
rise to a diffusion current of electrons from the n doped side of the junction to 
the p, and of holes in the opposite direction. 

The electron current from the n side removes the electrostatic shielding of the 
donor atoms, which become positively ionised. Similarly, acceptor ions on the p 
side of the junction are negatively ionised. The resulting distribution resembles a 
parallel plate capacitor with positively charged donors on one plate and negatively 
charged acceptors on the other. The region populated by ionised impurities is 
known as the space-charge region, and is made up of the p and n depletion layers. 

The electric field between the ionised dopants opposes the diffusion current. 
In the absence of an external bias, the p - n junction in the dark reaches an 
equilibrium in which the electrical field cancels the diffusion potential. It can be 
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shown [4j that this condition is equivalent to a constant Fermi level throughout 
the system, as illustrated in figure 12.2b . 

The thickness of these depletion layers can be found from charge neutrality 
and assuming that the total potential dropped is equal to the difference in Fermi 
levels between the two materials before electrical contact. 

This difference in Fermi levels AEp is equal to —qVu where Vu is the voltage 
dropped across the depleted region and is known as the built in voltage. It is 
conveniently expressed in terms of the intrinsic carrier concentration rii and the 
dopant densities. Subtracting equation 12. 151 from equation 12. 161 and substituting 
nf for the product NcNy from equation 12.101 gives 



The built in voltage is negative, in the sense that it is opposite to the standard 
terminal connection terminology, which is that the n terminal is negative and the 
p positive. 

For heavily doped high-bandgap materials, the Fermi levels in the p and n 
neutral regions lie close to the conduction and valence band edges, and Vu is only 
a few percent smaller than the bandgap of the semiconductor. 

Since the space-charge region has essentially zero mobile charge density in 
low injection conditions, we assume that all the external bias is dropped across 
it, and that no net field exists in the neutral sections. The total potential drop 
across the space-charge region is therefore given by the sum of the applied bias 
Vext and Vu 



This potential drop can also be derived by integrating Poisson's equation across 
the space charge region of width W = x^p + x^^ . 




(2.20) 



V = Vext + Vu 





(2.22) 



where e is the permittivity of the material. The charge density p on both sides 
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of the junction is given by Na and Nr,. 

Since there is no electrical field in the neutral sections of the p and n layers, the 
sum of fixed positively and negatively charges must cancel. The charge density 
in the p depletion layer is negative and equal to —qN^, where q is the magnitude 
of the electronic charge. It is exactly opposite to +qND in the n depletion layer. 
If Xwj^ and Xw„ are the p and n layer depletion widths, this charge neutrality 
condition can be written 

Xn^.NA = x^„Nd (2.23) 

By solving equation 12.221 subject to 12.231 and 12.211 we find that the depletion 
widths (with Vbi < 0) are 



1 



2e{-Vbi - Vext) ^2 24) 



1 



where e is the permittivity of the material. We note that the depletion widths 
increase as the applied bias becomes large and negative, and tend to zero as the 
positive applied bias cancels the built in potential. 



2.2.3 Current-Voltage Characteristic of a p-n Junction in 
the Dark 

The current fiowing through a solar cell is the sum of the photocurrent and of a 
diode characteristic, known as the dark current. The dark current component is 
important in determining the open circuit voltage of the cell. A good understand- 
ing of the dark current is therefore important in optimising the overall efficiency 
of the solar cell design. 

Studies of dark current in p - n diodes were pioneered by Shockley em et al. 
[7]. They established that an exponentially increasing current fiows if the diode 
is driven into forward bias by making the p contact positive with respect to the 
n. This is because the applied bias reduces the electric field in the space charge 
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region away from its equilibrium value. The drift current no longer cancels the 
diffusion current, and a net current flows. 

A negative applied bias on ther other hand increases the field in the space 
charge region, inhibiting the diffusion of majority carriers from the p and n sides 
of the junction. The diffusion current therefore decreases and no longer cancels 
the drift current due to generation-recombination effects in the neutral regions. A 
small, voltage independent saturation current Jo is seen to flow until the external 
applied bias exceeds a critical break-down voltage, at which the structure cannot 
sustain potential difference across the terminals and a sudden increase in reverse 
current is observed. 

Shockley et al. |7j initially modelled the dark current of an ideal diode in 
terms of semiconductor majority and minority carrier transport equations. This 
considers only the diffusion currents of holes and electrons from the neutral sec- 
tions of the diode and neglects the generation-recombination current in the space 
charge region. This approach yields the Shockley ideal diode equation 

J.deaK^) = '/oi(e^''/'"^-l) (2.26) 

where the saturation current Jqi is the ideal saturation current. It can be ex- 
pressed in terms of the electron and hole transport parameters. It is a propor- 
tional to nf and is therefore given in Sze jS] as 

Joi ~ 6^''/'=^^ (2.27) 

A more advanced analysis takes the generation-recombination current in the space 
charge region into account. The discussion in Sze [1] based on analysis by Sah, 
Noyce and Shockley [5] shows that the electron and hole emission is dominated 
by recombination centres situated at energies in the middle of the bandgap. The 
generation-recombination current Jgr can then be expressed in terms of a recom- 
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Voltage 

Figure 2.3: Semilog schematic of the forward bias Shockley dark current characteristic show- 
ing the generation-recombination regime (n=2) and the diffusion regime (n=l) 

bination lifetime Xr, the total space-charge region width W and rzj: 

J = SEHl (e^v/2ksT _ ^\ (2.28) 

The sum of equations 12.261 and 12.281 gives the Shockley-Read-Hall diode dark 
current Jb 

= Joi(e''^/'=^^ - 1) + Jo2(e«^/''=«^ - 1) (2.29) 

Equation 12.291 neglects the thermal contribution to the dark current, which has 
a bandgap dependence of ~ ^-Eg/ksT ■\^^^ jg negligibly small in semiconductor 
structures suitable for solar cell applications. The expression is frequently sim- 
plified as follows when modelling data: 

Jd ~ Jo(e^^/"'^^ - 1) (2.30) 

where n is called the ideality factor. It can be seen from equations 12.261 and 
12.281 that n = 1 for an ideal diode, and n = 2 for p - n junctions with strong 
generation-recombination currents. Figure [231 shows a schematic of a diode dark 
current in forward bias. 
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Diodes made from materials with small, direct band gaps tend to follow the 
ideal diode equation adequately, because the relatively large value of rii ensures 
that the diffusion currents from the neutral regions dominate the dark current. 

Diodes made from high bandgap materials on the other hand frequently show 
two regimes. At high forward bias, the diffusion current dominates and a factor 
n = 1 is seen. For low forward bias, the reverse is true and n tends to 2. 

2.2.4 The p-n Junction in the Light 

We have mentioned before that the photocurrent from a solar cell is made up 
of two components which flow in the opposite way to a conventional current. 
Minority carriers in the space-charge region are rapidly separated by the built 
in field, and swept to the interface with the neutral regions where they become 
majority carriers. 

Minority carriers generated in the neutral regions must reach the space-charge 
region before they can contribute to the photocurrent, because no field exists in 
these regions. Transport is dominated by carrier diffusion. 

Both these cases are examined in detail in section [31 The photocurrent Jl is 
expressed as the sum of diffusion components from the p and n regions Jp and 
Jni and a drift component from the space charge region Jj 



A significant decrease in the photocurrent is observed if the solar cell is forced 
into forward bias of the order of the built in potential, because charge separation 
in the space-charge region is impaired. A small increase is seen in reverse bias as 
the space charge region broadens because of the enhanced internal field. However, 
in the operating range of high quality solar cells, the photocurrent Jl may be 
assumed to be constant with bias and equal to the short circuit value Jsc- 

The current-voltage characteristic of an illuminated solar cell JiV) is then 
approximately given by sum of and the dark current J/j of equation 12. 30^ 




n 



(2.31) 
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Figure 2.4: Light and dark currents of a solar cell, showing the maximum power rectangle 



noting that the photocurrent is negative 



A typical light current characteristic is given in figure 12.41 

We now define some of the more useful quantities for solar cell characterisation 
and optimisation. The open circuit voltage Voc is the forward bias at which the 
dark current exactly cancels the photocurrent, and no net current fiows. Setting 
J = in equation 12.321 and assuming that .Jd is dominated by diffusion currents 
at such voltages such that Jgr « Jideai, we obtain 



The maximum power which can be extracted from the solar cell is given by the 
largest rectangle which can be fitted under the curve, as shown in the schematic 
figure [231 This rectangle defines the maximum power point, and the associated 
optimum operational current and voltage Jp and Vp. The maximum power is 
then 



(2.32) 




(2.33) 
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It is useful to define a fill factor EF as follows 

FF = (2.35) 

^ocJsc 

The overall efficiency of the cell rj is given by the ratio of the the maximum power 
Pmax and incident power Pj„c- It is usually expressed as follows 

V I FF 

V = — p (2.36) 

2.3 The p-i-n Solar Cell 

In the last sections we have defined the principal characteristics of p - n junction 
solar cells. We saw in section 12.2.21 that carriers photogenerated in the space- 
charge region are collected with efficiencies approaching 100%. It is therefore 
clearly desirable to extend the width of this layer. Although equations 12.241 and 
12.251 show that this can be achieved by reducing the p and n doping densities, 
this technique has two unfortunate consequences. The first is the reduction of 
the Voc according to equation 12.331 The second is an increased resistivity in the 
neutral regions due to lower doping. 

We now introduce a structure which circumvents these problems. This is the 
p-i-n diode, which is commonly used for high efficiency photodetectors because 
the spectral response of the cell can be optimised as a function of the response 
time. It serves as a useful starting point for the design of an efficient quantum 
well solar cell, as we shall see below. 



2.3.1 Space Charge Region of a p-i-n Diode 

The p-i-n solar cell consists of an ordinary p - n diode, with a layer of nominally 
undoped material inserted between the p and n regions. This layer is known as 
the intrinsic region, of thickness Xj. The total space-charge region thickness is 
then given by 

Wi = Xnjp +Xi + Xu,^ (2.37) 
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Figure 2.5: Electric field variation in an ideal p-i-n diode with zero background doping at 
zero and reverse biasses 

Figure [231 shows the electric field variation in such a structure at reverse bias 
and at zero bias. The p and n depletion regions decrease in size as the field drops. 
The field in the i region will only drop to zero when the V^xt is equal and opposite 
to Hi- 

Although the intrinsic layer is nominally undoped, we always find a residual 
or background doping density Nbg, which may be either n or p type depending 
on the growth machine and growth conditions. It is of the order Nbg ~ lO^^cm"'^ 
in good material, but may be as high as 10^^ or even 10^''cm~^. 

The calculation of Wi is similar to the calculation outlined in section I2.2.2I 
The built-in potential is given by equation 12.201 and by Poisson's equation 12.221 
as before. The main difference is that the double integral is over the new space 
charge region thickness Wi, and must include any residual space-charge due to 
the backgound doping Nbg- Figure [22] shows the variation of field as a function 
of position in a p - i - n structure where the background doping is p type. Curve 
A shows the field for a large negative applied bias. In the case of curve B, the 
magnitude of the negative applied bias is lower, and the depletion regions are 
smaller. Curve C shows where the i layer is not completely depleted. 

For the p-i-n cell with a n type background doping shown in figure 12. 6[ 
the potential V dropped across Wt is again the sum of external and built in 
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Figure 2.6: Electric field variation in a residually n doped p-i-n diode for A high reverse 
bias B moderate reverse bias and C an undepleted i layer 



potentials, and is defined as before by equations 12.201 and 12.211 It can also be 
calculated by integrating Poissons equation across Wi, 
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(2.38) 



The charge neutrality condition of equation 12.231 must be re-written to take the 
background doping into account: 



Xw„ 

NA + XiNBG + X^„ND=0 



(2.39) 



Equations l2.38l and l2.39l define the depletion widths in terms a quadratic equation. 
For the p side depletion layer of the n type background doped cell we have been 
discussing, this quadratic equation takes the form 
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(2.40) 



Similar expressions describe the n depletion width x^^. If the intrinsic doping is 
n type, the sign of the term in Nbg of equation 12.381 must be reversed. 
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2.3.2 The Background Doping Problem in p-i-n Struc- 
tures 

In this section we consider the influence of the intrinsic region thickness Xi and 
the background doping Nbg on the performance of a p - i - n solar cell, according 
to theory described by Paxman [8]. We then introduce a method of calculating 
Nbg from the dependence of the photocurrent Jl on voltage. 

A calculation of depletion widths as a function of intrinsic region width Xi and 
Nbg raises a number of potential problems with this type of structure. The n and 
p depletion layer thicknesses decrease clS X{ IS made larger, thereby decreasing the 
overall space charge region width. Calculations however show that the decrease 
in n and p depletion layer thicknesses is negligible compared to the increase due 
to the intrinsic layer. 

Another consequence of increasing the intrinsic region width is a lower elec- 
tric field in the space charge region, which simply results from the fact that the 
constant built in potential is dropped across a wider thickness. This sets a prac- 
tical limit on the amount of intrinsic material which can be inserted between the 
doped layers. This is not usually a serious limitation in practical solar cells. 

Studies of the space charge region thickness dependence on background doping 
reveal more challenging technological problems, as we shall see in the next section. 
In this case, the ionised background charge may cause the electric field to drop to 
zero in the intrinsic region. Since Nbg is always much smaller that Na and Nd, 
Poisson's equation shows that the rate of reduction of the space charge region 
thickness with increasingly positive voltage becomes much greater. Furthermore, 
cell resistance may increase because of the presence of an undepleted and low 
doped thickness adjoining either the n or the p region. Both these effects will 
reduce the solar cell efficiency. 

We now consider the variation of photocurrent Jl for p-i-n cells with high 
and low values of Nbg- In the case of low Nbg, a nearly constant photocurrent 
is seen into forward bias. A residual voltage dependence does in principle exist, 
and is due to the change in depletion widths as a function of applied bias. The 
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photocurrent in forward bias may drop because of experimental considerations 
which will be considered in chapters M and [71 In theory, however, it drops at 
the point where the applied bias is equal and opposite to Vbi- This removes the 
charge separating potential, and the device no longer operates as a solar cell. 

In the case of a cell with a high background doping, the higher space charge 
density in the intrinsic region imposes a significant electric field gradient, as 
illustrated by the curves of figure [231 In this case, the width of the space charge 
region decreases slowly while the p and n depletion widths are greater than zero, 
and Jl is only weakly dependent on the applied bias. 

As the applied bias is increased, either the p or the n depletion layer vanishes. 
At this point, there is still a field in the space charge region, and charge separation 
still occurs. However, the space charge region diminishes much more quickly with 
increasing bias because of the lower doping density in the space charge region. A 
significant decrease in Jl with increasing voltage is suddenly seen. 

The corresponding shoulder in the curve of Jl versus applied bias provides us 
with a useful tool for estimating Nbg in poor quality samples. At the shoulder, 
the p-i-n cell is effectively behaving like a p~ncell with the p type depletion layer 
given simply by the intrinsic region width, and an effective p doping A^^ = N^q. 
We can therefore substitute for and N^q for Na in equation 12.241 giving 
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(2.4i: 



where Vgh is the measured voltage at the shoulder. The equivalent relation for a 
n type N^q = is given by equation 12.251 
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(2.42) 



Equations 12.421 and 12.411 can be used to estimate the maximum tolerable back- 
ground doping for a given p-i-n solar cell design. For a typical structure with 
a p type background doping, — lO^^cm^^ and an intrinsic layer thickness 
Xi ~ O.Sfim we find that the highest background doping density permitting the 
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cell to operate adequately to voltages of about a volt is N^q ~ 5 x lO^^cm ^ 

2.4 Limiting Efficiency of Solar Cells 

In the previous sections we have introduced the quantities describing the efficiency 
of a solar cell. Equation 12.361 provides a starting point for designing a high 
efficiency solar cell. 

High efficiencies require large values of Jsc, Vqc and EF. The series resistance 
of the cell is important in determining EF and is ideally small. Although this can 
be achieved by high doping levels, practical limits are set by the degradation of 
transport parameters in heavily doped samples. 

A high Jsc clearly depends on absorption of as much light as possible. Equa- 
tions 12.181 and 12.191 show that a low bandgap increases absorption by reducing 
the cut-off energy below which no photons are absorbed. 

The requirements for a high Vqc can be found from equation 12.331 A fun- 
damental requirement for a low Jd and hence a high Voc is a low recombina- 
tion. Recombination can be reduced to some extent by improved cell fabri- 
cation techniques and passivation which reduce the concentration and activity 
of non-radiative recombination centres. The fundamental and unavoidable ra- 
diative limit however guarantees a lower limit to the dark current. We shall 
see in subsequent sections that GaAs/AlxGai_xAs solar cells are dominated by 
Shockley-Hall-Read non-radiative recombination mentioned in section 12.2.31 

We have seen that the the saturation current contains an exponential bandgap 
dependence of the form Jqi ~ ^-^a/^BT ^j^g ideal diode. 

The dependence of Jgc on Eg is more complex. Experimental and theoretical 
studies (see chapters [3] and [8]) however show that it decreases with increasing 
bandgap. This can be partly explained by the change in absorption. An increase 
in bandgap moves the absorption edge to higher energies and decreases the ab- 
sorption coefficient above that threshold (see equations l2.18l and l2.19p . The direct 
result is a reduction in minority carrier generation and therefore a smaller light 
current. This decrease however is weaker in character and magnitude than the 
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Figure 2.7: Light currents of solar cells with low and high bandgaps under similar illumination, 
showing the compromise between voltage and current. Also shown is an optimum design with 
a manifestly greater maximum power rectangle. 

exponential dependence of the saturation current. The form of the Vqc therefore 
reduces to 



A compromise therefore exists between the high bandgap requirement for a high 
Voc and the low bandgap for a high Jgc, as illustrated in figure l2n 

Limiting efficiencies for solar cells as a function of bandgap have been cal- 
culated by Henry [9]. The figures quoted are for the AMl.Sd standard solar 
spectrum. This spectrum corresponds to direct sunlight at noon at a latitude of 
48°. Other standards include the AMO spectrum, which is sunlight without any 
atmospheric absorption. 

An optimum of 30.6% reported by Henry for materials with bandgaps in 
the region of 1.35eV, which is near to the GaAs bandgap (1.424eV) for direct 
semiconductors, and to Si (1.12eV) for indirect materials. Higher efficiencies 
can be achieved by concentrating the incident light, which increases both the 
short circuit current and the open circuit voltage and increases the maximum 
efficiency to about 37.2%. The concentrator approach is limited by cell heating 
under the concentrated light. The increased operating temperature increases the 



Voc ~ Eg ln( Jsc) 



(2.43) 
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dark current, thereby reducing the open circuit voltage. Cells for concentration 
must therefore be designed with this effect in mind, or cooled. 

Cells for concentrator applications must also be designed with a much lower 
series resistance, because of the higher current densities. Although cell resistances 
can be reduced by increasing the doping, this leads to a reduced minority carrier 
mobility because of the increased number of scattering centres. 

Cells with efficiencies in the high 20 's have been achieved in practice. Stan- 
dard tabulations of world record solar cell efficiencies under AMI. 5 ilumination 
published by Green [TI] quote the best GaAs solar cell as having an efficiency of 
~ 25%, and ~24% for Si. Laboratory cells are therefore close to achieving the 
fundamental efficiency limits, and increasingly complex designs are necessary to 
achieve higher efficiencies in single bandgap systems. 

Higher efficiency designs however, which are not subject to the same efficiency 
limits as single band-gap solar cells, have been suggested. These are reviewed in 
the next section. 

2.4.1 Higher Efficiency Designs 

We have seen that the two fundamental factors limiting the maximum efficiency 
of single bandgap solar cells are the trade-off between Vqc and Jsc, and the related 
problem of thermalisation losses. 

A solution to this trade-off is to combine several band gaps in a single solar cell. 
Different parts of the incident spectrum are absorbed in solar cell components 
with appropriate bandgaps. The component with the highest bandgap absorbs 
short wavelength light and delivers current at a high voltage. The cell with the 
lowest bandgap absorbs long wavelength light and operates at a lower voltage. 
The operational principles of the component cells remain those of single bandgap 
cells. The combination of the different bandgaps however reduces thermalisation 
losses, resulting in higher efficiency limits. This concept is known as the rainbow 
or cascade cell. Henry [9] predicts a maximum theoretical efficiency of around 
75%. 
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Unfortunately, this design suffers from serious technological problems. If the 
cells are not connected in series, the output voltage is that of the cell with the 
lowest bandgap. On the other hand, connection in series requires that the current 
outputof all the components be equal, because the output current of the cell will 
be that of the worst component. 

Finally, series connection of cells is technologically difficult. It can be achieved 
at high cost in complex electrical connections, or during growth by the incorpo- 
ration of tunnel junctions. 

A simplified approach is the tandem cell. This exploits the same principle 
as the rainbow cell but restricts the design to two materials, thereby simpli- 
fying the current matching and electrical connection problems. Henry [9] pre- 
dicts a fundamental efficiency for this system of 50%. A recent indium-gallium- 
phosphide/gallium arsenide tandem cell reported by Bertness [lOj has achieved 
a record efficiency of 29.5% in AM1.5g (global) illumination. 

2.5 Quantum Well Solar Cells 

In previous sections we have described the operating principle of the photovoltaic 
effect and given an overview of fundamental loss mechanisms. We have reviewed 
a number of higher efficiency designs which can minimise these losses and lead 
to higher efficiencies. 

We now present the main design studied in this thesis, which is the quantum 
well solar cell (QWSC), first suggested by Barnham [12J. The QWSC is a p - i 
- n cell with quantum wells grown in the intrinsic region. This cell is similar to 
tandem cells because it contains material with two different bandgaps. 

Figure 12.81 shows a schematic of the QWSC. The structure contains two 
bandgaps Eg and Ea which are the barrier bandgap and quantum well bandgap 
respectively. The material behaves like bulk material for electron and hole ener- 
gies greater than Eg and like a quantum well or superlattice states for energies 
below Eg. 

For incident photons with energies greater than Eg, the QWSC operates like 
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a conventional single bandgap p - i - n cell. 

Photons with energies between Eg and Ea however are absorbed in the quan- 
tum well states. Once excited into well states, they can recombine radiatively or 
non radiatively or escape from the well. The escape mechanisms are thermionic 
emission into a bulk state, tunnelling, or a combination of the two. Experimen- 
tal and theoretical studies of this system by Nelson [13] have shown that the 
escape probability from GaAs quantum wells embedded in AlxGai_xAs with a 
composition of 34% is about 100% at room temperatures. This very high escape 
probability is seen even when the field across the quantum well is low. Thermionic 
emission therefore is sufficient to ensure high carrier escape probabilities from the 
wells at room temperatures. 

Further work by Barnham et al. [H] has shown that the output voltage is 
greater than would be expected from a single junction cell with the quantum well 
bandgap Ea- This has been observed in three different materials systems. In the 
AlxGai_xAs system, this work showed that increasing barrier bandgaps Eg give 
increased voltage enhancements over the Vqc expected from a single bandgap cell 
with the quantum well effective bandgap Ea- 

The QWSC design therefore departs from the analysis used by Henry [9] in 
his analysis of the fundamental efficiency limit. Unlike the separate components 
of a single bandgap cell, the QWSC operates in a physically different way to 
single bandgap cells because of the minority carrier escape from the quantum 
wells. Since it is in principle capable of generating the same photocurrent as a 
single bandgap cell with a higher voltage, the design has a higher fundamental 
efficiency limit than a single bandgap cell. 

The fundamental efficiency limits of QWSCs are not yet clearly understood. 
The review article by Nelson gives an overview of the subject. 

2.6 The QWSC in the AlGaAs System 

Work by Paxman [8j has demonstrated efficient carrier escape from QWSC struc- 
tures in AlxGai_xAs with x ~ 30%. This resulted in a sample grown by Tom 
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Fi gure 2.8: Schematic QWSC structure, absorption, recombination and escape mechanisms 



Foxon at Phillips Laboratories, Redhill, which was 14.2% efficient in AMI. 5. 
This sample showed a good Voc of 1.07V, and a fill factor FF = 77%. More detail 
concerning optimisations carried out in this work will be given in chapter [51 

The work by Paxman however showed scope for improvement in the J^c of 
QWSCs in AlxGai_xAs. The advantages of this material include extensive char- 
acterisation and a material which is lattice matched for compositions ranging 
from GaAs to AlAs. The lattice spacing in GaAs and AlAs is only different by 
about 0.14% according to Adachi ^16j . 

The material however suffers from poor transport characteristics. As we shall 
see in chapter HJ minority carrier transport is poor, and the Jgc of single bandgap 
AlxGai_xAs cells decreases rapidly with increasing aluminium composition x. 

We conclude that QWSCs in the AlxGai_xAs system mainly suffer from poor 
minority carrier characteristics, which result in a low J^c- Work presented in fur- 
ther chapters will address this problem with theoretical and experimental analysis 
of design changes made to the basic QWSC design. 
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2.7 Conclusions 

This chapter has reviewed the operating principles of solar cells and defined the 
main parameters allowing characterisation of their performance. 

We have seen that single bandgap solar cells are limited to efficiencies in un- 
concentrated sunlight of about 30%. by fundamental loss mechanisms. Although 
designs such as the cascade cell and variations on this design are capable of 
reaching higher efficiencies, they are subject to complex technological problems 
because they combine two or more component solar cells in one design. Moreover, 
the approach represents a technological solution rather than one which introduces 
fundamental changes to solar cell operation. 

The QWSC solar cell however is governed by physically different efficiency 
limits, because it has a which is higher than that of a single bandgap cell 
with the quantum well effective bandgap. This principle has been demonstrated 
in three different materials systems. 

The AlxGai_xAs system shows a number of advantages for quantum well 
applications. Unfortunately, it suffers from poor minority carrier transport prop- 
erties, which degrade with increasing aluminium fraction. This results in a poor 
short circuit current, and poor efficiencies for solar cells fabricated from this 
material. 

We conclude that the voltage enhancement aspect of designing an efficient 
QWSC in the AlxGai_xAs cell has been demonstrated, but that the Jgc must 
be increased. The following chapters will present a theoretical and experimental 
analysis of this problem. 



Chapter 3 

Quantum efficiency and short 
circuit current models 

3.1 Introduction 

Optimisation of the short circuit current (Jsc) requires an understanding of the 
spectral response of the cell at zero applied bias. In this chapter we give a brief 
survey of previous theoretical work on such models, before going on to describe a 
model of the Jsc for GaAs/AlxGai_xAs QWSCs with position dependent trans- 
port and optical properties. Values of the materials parameters are discussed in 
chapter HI 

3.1.1 Review of Previous Work 

Modelling of AlxGai_xAs solar cells has largely been confined to theoretical stud- 
ies of cases which can be solved analytically. Hovel |69] gives a good review of 
solar cell physics, including discussions of a range of analytical and numerical 
models. 

The theoretical model by Konagai [20] considers a p AlxGai_xAs - p GaAs 
n GaAs compositionally graded p layers grown epitaxially on GaAs. He does 
not consider the position dependence of the transport parameters, and uses an 
analytical expression for the absorption coefficient. The effect of the grade is 
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expressed via the quasi-electric minority carrier field it produces in the graded 
region. The field is taken as the gradient of the bandgap due to the compositional 
change, and is discussed in more detail in chapter |H This work suggests that the 
grade is especially useful in reducing the minority carrier recombination at the 
surface of the cell. Efficiencies of the order of 20% are predicted under air mass 
zero (AMO). 

Hutchby [21] has studied a similar design, but for an n AlxGai_xAs - n GaAs 
p GaAs design. This work is again purely theoretical, and makes similar as- 
sumptions by assuming that the minority carrier transport is constant. The 
absorption coefficient is a semi-empirical fit to data. The quasi-electric field in 
this model is again taken as proportional to the band-gap gradient, and the mi- 
nority carrier transport parameters are assumed to be essentially independent of 
the aluminium composition. The theoretical efficiency in this case is 17.7% for 
an AMO spectrum. 

Sutherland [22] has presented a theoretical study which uses semi-empirical 
parameterisations to model the position dependence of the materials parameters. 
He finds good agreement with the earlier studies, predicting a peak efficiency of 
about 20% for a n AlxGai_xAs - n GaAs p GaAs design. 

Hamaker [23] has presented a theoretical model of a n - p and p - n AlxGai_x- 
As cell with graded emitter layers. This is the most accurate model published 
to date, and uses semi-empirical parametrisations combined with data for the 
materials parameters. An optimisation in terms of doping and compositional 
profiles is included. The predicted efficiencies are similar to those predicted 
by the previous studies. Very little difference however is reported between the 
different designs studied. 

Of the work mentioned above, the Hamaker ^24j model is closest to the one 
we shall present in this chapter. We note however that the models, with the 
exception of Hamaker p3] have used analytical approximations for the materials 
parameters and generally only consider the position dependence of the optical 
parameters. 

The models mentioned here are not compared with experimental data. In 
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the following sections we will present a model similar to that of Hamaker [2l] 
but which is equipped to deal with quantum wells and compositional grading 
throughout the structure. The model furthermore will be designed to repro- 
duce experimental data, and the relevant fitting parameters identified. The main 
methodology has been described by Connolly [25] . 



3.2 Preliminary definitions 

We start by introducing the external quantum efficiency or QE. It describes the 
spectral efficiency of the cell independently of the particular incident spectrum 
and is defined as the probability of an incident photon of wavelength A generating 
an electron-hole pair which contributes to the Jsc- It can be expressed as the 
ratio of the number of steady state current carriers Jsc/q collected under an 
illumination of wavelength A to the number of photons incident on the top surface 
of the cell F{X) where q is the magnitude of the electronic charge. 

QEW ^ (3,1) 

We also define the internal QEint, which is defined in the same way, but where 
the flux F is replaced by the flux actually transmitted into the solar cell. The 
two differ only by a factor dependent on the surface reflectivity. 

The model calculates the QE, which can then be used to predict the Jsc for 
a given cell design and illumination by multiplying the QE by the appropriate 
spectrum. 

The wide range of QWSC and p - i - n solar cell designs considered in the 
characterisation and optimisation program requires a general approach to each 
part of the calculation. These designs include cells with compositionally graded 
neutral and space charge regions. These will be referred to as graded cells. Cells 
without these compositional changes are referred to as ungraded cells. We also 
consider cells which operate as Fabry- Perot cavities because they feature a mirror 
which is deposited on the back surface of the cell. These will be referred to as 
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mirror backed cells. 

The first problem we consider in this chapter is the variation of the light 
intensity in the solar cell as a function of position, wavelength, and front and 
back surface reflectivities. The generation rate can then be calculated from the 
absorption coefficient and the light intensity at any particular position and wave- 
length. 

We then describe the photocurrent calculation in the two cases of depleted 
space-charge regions and neutral regions. In a general p - i - n cell, the space- 
charge region is made up of the p depletion region, the nominally undoped i 
region, and the n region depletion layer. Under the depletion approximation 
described below, the photocurrent calculation reduces to an integral of the gen- 
eration rate over the depletion region and wavelength. 

The second case deals with the photocurrent from the electrostatically neutral 
sections of the n and p layers. Solutions to the resulting transport equations are 
first discussed for ungraded cells where materials parameters are constant, and 
an analytical solution can be found. 

In the case of graded cells however an analytical solution cannot in general 
be found and the transport equations must be solved numerically. A standard 
difference method for this case is presented in the final part of the chapter. 

3.3 Calculation of light intensity in a QSWC 
3.3.1 Definitions 

In this section we derive the light intensity in a solar cell of total thickness d 
and in which the optical parameters may be position dependent. We define 
the optical parameters required to calculate the light intensity for cells where 
interference effects may be significant. We consider only the normal incidence 
case. A comprehensive discussion is given in Macleod ^35j. 

The first phenomenon experienced by a beam of light of intensity Xj which 
is incident on a solar cell is refiection. We quantify this by the front surface 
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Figure 3.1: Example generalised QWSC band diagramnic!, including position dependent 
bandgap in the p type layer and high aluminium fraction window layer 

reflectivity R which is the fraction of light intensity of wavelength A reflected by 
that surface. 

For graded material the absorption coefficient varies with distance. In this 
case the intensity at a particular point depends on the integral of the absorption 
which the light has suffered up to that point, minus the fraction of light reflected 
at the front surface. 



X{x, A) = X,(A)(1 - i?(A))e/o^"(^''^^'^^' 



(3.2) 



This expression allows us to calculate the light intensity for cases where hght 
only makes one pass. By this we mean that the amount of light internally reflected 
from the back surface of the cell is negligible. 

This is not the case for cells which are coated on the back surface with a 
mirror and which leave a significant amount of light unabsorbed on the first pass. 
In this case, successive reflected beams suffer phase changes with distance and 
upon reflection at the front and back surfaces. The superposition of beams of 
light with different phases leads to interference phenomena. To take these effects 
into account, we must derive the intensity from the complex amplitude of the 
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light. The electric field component of the light £^ is a convenient quantity to use 
and we refer to it as the amplitude for brevity. It is of the same form as X and 
is characterised by a complex decay constant k. The real part of k characterises 
the exponential decrease in electric field strength due to light absorption, and the 
imaginary part refiects the change in phase of the electric field with distance. S 
is then 

S{x,X) = ^o(A)e-/o^'=("''^)''"' (3.3) 

where Sq is the amplitude at position zero. 

The light intensity X is proportional the product of the electric field amplitude 
and its complex conjugate [35] : 

X(x, A) = evSS* (3.4) 



where S* denotes complex conjugate of 8, and e and v are respectively the per- 
mittivity and the speed of light at position x. 

Sq is related to the incident fiux Xj by substituting the value of X(0, A) from 
equation 13.21 into equation 13.41 Since the phase changes in the cell are all relative 
to the phase of the incident light, we choose to set this incident phase to zero. 
Sq is then real and of the form 



- (3,5) 

V ev 

The total amplitude S is found by summing the amplitude contributions from 
successive reflections £j, where i is the order of reflection. 

oo 

The dimensionless components Sj are a measure of the amplitude of the j^^ 
amplitude component relative to the amplitude Sq transmitted throught the front 
surface of the cell, and are calculated below. 

Equations 13.61 and 13.41 together deflne the intensity as a function of distance 
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and wavelength in terms of the ampUtude components Sq. 

3.3.2 Phase and reflectance 

Before calculating the components £j we must relate the reflectivity and absorp- 
tion coefficient to the corresponding quantities for the amplitude. Both of these 
can be derived from inspection of equation I3.4[ 

When a normally incident beam of light of intensity X is internally reflected, 
the intensity is reduced by a factor R. Consistency with equation 13.41 demands 
that the corresponding amplitude should be modified by a factor a/R. A factor 
e*'^ describes the phase change upon reflection at the back surface, where i = V—1 
and ip is the phase change. The complex reflectance is therefore the square root 
of the R multiplied by the ip. 

r = e'^y/R (3.7) 

The relationship between amplitude and intensity in equation 13.41 together with 
equations 13 . 21 and 13.31 shows us that the real part of k is simply half the absorption 
coefficient a. However, it is a complex quantity and may have an associated 
phase. 

The phase of an electromagnetic wave changes by a factor 27r over a distance of 
one wavelength. Furthermore, the wavelength of light in a medium with refractive 
index nis X/n. The phase change is therefore simply the ratio 27rm/A. The phase 
change over a distance x will be the integral of this phase over distance. 

The complex decay constant is then the sum of a/2 and the phase 

, , , , a 27r m , , 

Hx,X) = --^ (3.8) 

3.3.3 Generation Rate 

We can now explicitly write down the successive terms in the sum of amplitudes 
as a function of position x. The first is the simply the transmitted beam and 
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starts with the relative amplitude of 1. 



(3.9) 



The initial amplitude of the second beam is equal to the limit of £i at the back 
mirror, but with a reflectance factor. We introduce the back surface reflectance rb 
and it's associated phase change ipb which is related to the back surface internal 
reflectivity Rh as outlined above. The amplitude of the second beam is then 
dependent on an integral over —dx from the back surface d to the position x. 
This is the same as integrating forwards from x to d over +dx. Component £2 is 
then 

If r/ is the front surface internal reflectance, with associated phase ipf and reflec- 
tivity Rf, subsequent reflections take the form 



£4{x) 
£b{x) 



(3.11) 



Rewriting the even reflections in terms of an integral from to d minus an integral 
from to a; and factorising by the extra term in r^, we obtain two geometric series 
which differ only by a factor relating to the amplitude attenuation over one pass 
and the back surface reflectance. It has the limit 
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(3.12) 



The intensity in the solar cell can be written independently of the factor ev as 
follows 



I{x,X) =Xi{l-R) 



' ££* 
.£o£o 



(3.13) 



We emphasize here that this treatment consitutes a flrst approximation to the 
problem, since it does not consider the internal reflections taking place at the 
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quantum well interfaces. The reflectivity at these surfaces, however, is not great, 
since the change in refractive index between GaAs and AlxGai_xAs is not large for 
the barrier compositions considered. A brief calculation using optical parameters 
reviewed in chapter H] indicates that the reflectivity in the worst case of a 45% 
aluminium barrier is only 0.2%. On the other hand, a large number of wells may 
make this effect significant. Further discussion of this point is given with respect 
to experimental data in chapter ID 

The generation rate G is then readily given by the product of a and / as 
follows 



3.4 Photocurrent 

Having derived an expression for the generation rate we now investigate the 
mechanisms which transform these photo-generated carriers into a useful electric 
current. 

In general, the depleted or space-charge region of the QWSC includes the p 
depletion region of width x^^ (see figure 13. ip . the intrinsically doped i region 
(width Xi), and the n region depletion layer (width Xw„). We have seen in section 
12.31 that the i region may not be fully depleted either if the forward applied bias 
is larger than the magnitude of Vti, or if the background doping level is high. In 
this case, the electric field drops to zero within the intrinsic region. 

If, however, the background doping level is low, a field is maintained across 
the i region in the moderate forward bias regime and carrier transport is drift 
dominated. Under the depletion approximation described below, carrier recom- 
bination in the space-charge region is negligible and charge collection efficiencies 
are 100%. 

In the neutral window, p and n regions, transport is dominated by diffusion 
and recombination is not negligible. The photocurrent from these layers is found 
by solving minority carrier transport equations. Previous work by Braun [18] has 
shown that the high band-gap window layer does not contribute a measurable 




(3.14) 
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photocurrent and may be treated simply as a non-contributing surface absorber. 

In the following sections the photocurrent calculations for space charge and 
neutral regions are discussed for the general case of a p - i - n structure with 
position dependent parameters. We then examine the simplified case of an un- 
graded structure and the manner in which certain parameters can be determined 
from this simplified case. We then develop a numerical method to solve the more 
general graded case. 

3.5 Intrinsic region photocurrent 

Calculation of the photocurrent Jj from depletion layer is based on the deple- 
tion approximation. This assumes that minority carriers generated or injected 
into the i region are swept across it to the boundary with the depletion edge 
where they become majority carriers. That is, electrons are swept towards the n 
interface, whilst holes are swept to the p interface with no loss. This is a good 
approximation for cells in which a field is maintained across the i region. For cells 
with either compensated or low background doping, this condition is maintained 
into moderate forward bias. J, can then be found from the generation rate. 

The photocurrent contribution from the quantum wells is calculated assuming 
an escape efficiency of 100%. Once the carriers have escaped from the wells, they 
behave like other minority carriers. Previous work by Nelson [13] has shown this 
to be true for GaAs wells in AlxGai_xAs with an aluminium fraction of 30%. 

The photocurrent density Ji from the space charge region is then simply given 
by the integral of the generation rate 



Although in general an analytical form for Jj can only be found for ungraded 
solar cells, the calculation is identical for the graded case. The quantum well 
absorption coefficient is described by Paxman [19]. 




(3.15) 



3.6. p and n Region Contributions 



59 



3.6 p and n Region Contributions 

The photocurrent from neutral p and n regions is found from current and conti- 
nuity relations outlined below. Since the method in both cases is identical, the 
following discussion concentrates on the p for clarity. 

3.6.1 Drift and Diffusion Currents 

Minority carrier currents can be caused by drift or diffusion. The former results 
from the acceleration of carriers in an electric field, whereas the latter is the 
consequence of a carrier concentration gradient. In the low injection limit and 
for moderate electric fields, the net current flowing at depth x in the p layer jp{x) 
is the sum of these two currents 0]. 



where /i„ is the minority electron mobility, D„ the electron diffusivity and E 
the quasi-electric field. These parameters are position dependent in the general 
case. We note here that some workers use the terminology jp as the minority 
hole current in the n layer. We choose to define it in this way in order to keep 
the same subscript between the excess minority carrier concentration Up and the 
corresponding current jp. 

The total minority electron carrier concentration Up is defined as the sum of 
the equilibrium carrier concentration n^^ and the excess minority carrier concen- 
tration n as follows: 




(3.16) 



(3.17) 



No net current flows in equilibrium, and equation 13.161 gives 
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Subtracting 13.181 from I3.16| we obtain 



n{x) + qDn{x) 



dn 
dx 



(3.19) 



where we have used the Einstein relation /i„ = qDn/kBT. This is vahd for 
equihbrium [T7] where the diffusion current due to the electron and hole carrier 
concentration gradients can be equated with the drift current flowing in the 
opposite direction. We assume that in the low field diffusion dominated regime 
the mobility and diffusion coefficient are equal to their equilibrium values. This 
is confirmed by low field mobility measurements by Walukiewicz [TS] 

Equation 13.191 allows us to calculate the photocurrent in terms of the excess 
minority carrier concentration n only, because the intrinsic carrier concentration 
gives no net current. 

The total current contribution of the p layer is the current injected from the 
p layer into the i region. In the depletion approximation, n is zero at the edge of 
the space charge. In this case this is at the junction between the electrostatically 
neutral section of the p region and the p depletion layer. Equation 13.191 then 
becomes 



Calculating Jp therefore reduces to the problem of determining the gradient 
dn/dx at the depletion edge from the current and continuity equations. 

3.6.2 Current and Continuity 

A minority electron in the p layer can be generated, drift, diffuse or recombine 
with a hole. The rate of change of current as a function of position is given by 
the difference between generation and recombination rates: 




(3.20) 



dx 



qG{x, X) + q 



n 



T 
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In AlxGai_xAs the recombination time r is dominated by non-radiadive Shockley- 
Read-Hall recombination mentioned earlier in section I2.2.3I If we differentiate 
equation 13.191 and substitute for r using L„ = y/rDn, we obtain the following 
differential equation after dividing throughout by Dn 

(fn , qE D'dn , qED' qE' ^ G , , 

where D'^ = dDn/dx. In order to clarify subsequent discussions, we re-write 
equation 13.221 as follows 

d^n , , ^dn , ^ , ^ 

-T^ + b{x)— + c{x)n + g = {) (3.23) 
dx'^ dx 

where functions h and c and the generation function g are defined by comparison 
with equation 13.221 Given suitable boundary conditions, equation 13.231 defines 
the minority carrier density as a function of position. 

3.6.3 Transport Equations for n Layers 

For completeness we briefly indicate the differences between the transport equa- 
tions for p layer and n layers. For an n type semiconductor, the minority carriers 
are holes with a positive charge of q. The drift current takes the same form, but 
the diffusion current changes sign. Equation 13. 161 becomes 

d'Pn 

jn{x) = qfin{x)E{x)pn{x) - qDp{x)-^ (3.24) 

The differential equation 13.221 however does not change except the the hole dif- 
fusion length Lp substituted for L„, and Dp and p„ are substituted for Dn and 
Up respectively. 

Since the solution to the problem is identical in both cases, we shall not 
consider the n layer contribution explicitly in the following sections in order to 
avoid unecessary repetition. 
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3.6.4 General Analytical Solution 

Equation 13.231 is an inhomogeneous differential equation which has a standard 
solution described for example by Boas [26j. A study of this method however 
increases our understanding of the influence of the parameters governing the 
behaviour of the solution. In the dark, equation 13.231 reduces to 

^ + h{x)— + c{x)n = () (3.25) 
ax^ dx 

A general solution to this homogeneous second order linear differential equation 
exists of the form 

ndark = aiUi{x) + a2U2{x) (3.26) 

where oi and 02 are constants fixed by two boundary conditions and Ui, U2 are 
functions of position. 

Following the methods outlined in Boas [2^, we postulate that the solution 
to the inhomogeneous case 13.231 is of the form 

Us = A{x)Ui + B{x)U2 (3.27) 

where A and B are two arbitrary functions which can be found with standard 
solution methods, and take the form 

^(^) = u,ul-u[U2^^' ^3 

These two functions together with the dark solution form the particular integral 
solution 

The general solution is then 



n{x) = aiUi + a2U2 + U3 



(3.30) 
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The constants ai and 02 are found as described in the next section. In gen- 
eral, analytical solutions cannot be found and we resort to numerical methods of 
solving equation 13.231 described below. 

3.6.5 Boundary Conditions 

The first boundary condition is that the surface current flow is equal to a recom- 
bination current characterised by a recombination velocity Sn 

JpiO) = qSnuiO) (3.31) 

The surface recombination may be either due to interface recombination centres 
at an abrupt junction between layers with different aluminium compositions or 
may be a surface recombination associated with defects present at the interface 
between the semiconductor and the air (Hovel |69j). 

Real devices may depart from this behaviour because of increased recombina- 
tion at the front contacts [19j or damage to the front surface. Both these effects 
may lead to very high recombination velocities, and a possible reduction of the 
p layer thickness. This is discussed further in section 14.4.51 

Substituting equation 13.191 for the current and the general solution 13.301 for n 
into the previous expression we obtain 

ai (z^„f/{(0) + (^-5„)f/i(0)) 

+a2 {Dr^Um + - 5n)f/2(0)) (3.32) 

= ^3(0) {Sr^-f^-DMm) 

The second boundary condition states that the minority carrier concentration 
at the boundary with the depletion region is zero, according to the depletion 
approximation 

n.^^ = (3.33) 
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Susbituting for n from equation 13.301 we obtain 

aif/i(x^p) + a2t/2(x^J = -U-six^^) (3.34) 

The constants ai and 02 are then obtained from the two coupled equations 13.321 
andE31 

Two different cases requiring two methods of solving equation 13.231 are pre- 
sented in the following sections. The first covers the case of an ungraded solar 
cell, for which there is no electric field in the neutral regions and materials param- 
eters can be assumed to be constant. A hypothetical cell with constant materials 
parameters but with an effective electric field applying only to minority carriers 
is examined as a preliminary to the general compositionally varying case. 

The second case is that of compositionally graded QWSCs, for which all 
the parameters in equation 13.231 are position dependent. The same boundary 
conditions are applicable in both cases. 

3.7 Constant Transport Parameters 

In this section we consider the p layer QE of ungraded cells with and without 
effective fields. In this case, the parameters L„ and Dn are constant. The original 
work is described in the recent review by Nelson [15j. 

3.7.1 Zero Field 

For a layer with a zero field, a trial solution e^^ yields the functions 

Ui = e^/^" (3.35) 

Ui = e"^/^" (3.36) 

Constants ai and 02 are obtained from equations 13.321 and 13.341 thereby defining 
the minority carrier profile. The photocurrent from the p layer is then found 
from equation 13.201 
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Appendix |X] gives more detail, and shows that the photocurrent for these 
ungraded structures is only weakly dependent on the diffusion coefficient and the 
recombination velocity if 

(3.37) 

The appendix further shows that the QE in the limit of L„ < x^p near the band 
edge is well approximanted by 

QE ~ L^ae""^-'' (3.38) 

This result agrees with a brief treatment given by Ludowise [27], and reviewed 
in more detail by Hovel [69]. The AlxGai_xAs structures we will study generally 
obey this limit for acceptable p layer thicknesses. The photocurrent near the 
band edge is dominated by p and i contributions. Since, as we have seen, the 
i region photocurrent is well understood, L„ can be estimated by modelling the 
QE in this wavelength range. 

Appendix |X] further shows that the QE at short wavelengths remains sensi- 
tive to the recombination velocity, and tends to a limit determined primarily by 
the surface recombination velocity and the absorption coefficient. For diffusion 
lengths which are of the order or greater than the p layer thickness, a simple 
limit is found, depending solely on absorption coefficient and recombination ve- 
locity. Since the absorption coefficient varies slowly at such wavelengths, the 
overall QE becomes weakly dependent on wavelength. This behaviour is primar- 
ily observed in solar cells made from materials with high absorption coefficients 
at short wavelengths, such as GaAs. 

3.7.2 Non Zero Effective Fields 

A brief examination of the hypothetical case of an ungraded cell with electric 
fields present in the neutral regions is useful here. Detail is again provided in 
appendix lAl We now define the effective field £f due to the compositional grade. 
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in units of per metre, as follows: 

^f = ^ (3-39) 
A trial solution for cells with £f yields the solutions 

Ui = e^^i^ (3.40) 
U2 = e^^"" (3.41) 

with 



ki = +J£} + L-2 - Sf (3.42) 



k2 = -^£] + L-^-£f (3.43) 

Since Sj < 0, parameter ki tends to 2 \Sf\ for large fields, whereas ^2 tends to 
zero. The exponentially growing solution is therefore favoured, and the resulting 
minority carrier distribution shifted in a positive direction towards the depletion 
layer, increasing minority carrier collection. 

We now examine the effect of the field on the surface boundary condition. 
Appendix Rl shows that the condition of 13.371 is replaced by 

(3.44) 

In this limit, the QE is insensitive to the surface boundary condition and the 
surface recombination velocity in particular. This limit is readily achieved even 
for high surface recombination velocities. We shall see in section H] that it may 
fail however if the recombination velocity is unusually high. This may occur if 
recombination at the contacts is large, or if the surface quality is poor. 

We conclude however that the presence of a field renders the QE relatively 
insensitive to the surface recombination. 

We shall see in subsequent sections that the real picture is somewhat less clear 
in the case of a full numerical solution to equation 13.221 We shall see however in 



3.8. Graded QWSCs and Numerical Solutions 



67 



the modelling section that cells with electric fields are much less sensitive to the 
magnitude of the surface recombination velocity, as carriers in a well designed 
cell tend to diffuse in the direction of the greater effective diffusion length. 

3.8 Graded QWSCs and Numerical Solutions 

Analytical solutions to equation 13.301 generally cannot be found if the parameters 
b and c are position dependent. We now develop a numerical scheme to solve 
equation 13.231 in this case. 

We express the first and second differentials of n using standard centred dif- 
ferencing described by Potter |30j as follows: 

(3.45) 

where k is the position index and A the grid spacing of a one dimensional grid 
of length N. Substituting these expressions into equation 13.231 we obtain 

^k-i{^ - ^) + Mck -^) + ^k+i{^ + ^) = "^^ ^^-^^^ 

The surface boundary condition equation 13.311 can be written as 

The second boundary condition equation 13.331 becomes 

UN = (3.48) 
Equations I3.46[ 13.471 and 13.481 can be expressed as the following matrix equation 

MN = G (3.49) 



This equation is then solved by standard Gaussian elimination with backsubsti- 
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tution. A good review of this technique is given in Press [29] 



Ahhough more complex methods were investigated, this simple method was 
chosen because trials with a number of different numerical schemes shows that 
this method is preferable to methods which are nominally quicker and avoid 
explicitly calculating the inverse matrix. This is mainly because the matrix M 
need only be inverted once, but is typically used about six hundred times to 
cover the wavelength range we are considering. For the matrices we consider, 
which may have a dimension of a few hundreds, we observe that the calculations 
involved in solving equation 13.491 are the most time consuming. 

3.8.1 Surface Boundary Condition in Compositionally Graded 



The surface boundary condition is expressed by the topmost row of matrix M, and 
is identical to the boundary condition in the analytical case. We conclude that 
compositionally graded designs remain insensitive to the surface recombination 
velocity. 

More detailed conclusions regarding the general behaviour of the QE of such 
cells are best illustrated by numerical calculation, and are investigated further in 
chapter [5l 

3.8.2 Numerical Accuracy 

To determine the numerical accuracy of the scheme we express the minority 
carrier concentration as a Fourier series 



Cells 



oo 



(3.50) 



where kj is the wavevector of the j mode. The exact differential of an arbitrary 
mode e*'^^ is 

n' = aikn (3.51) 
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When we use our difference expression of 13.451 for tlie same differential we obtain 

~ (3 52) 

^(2^sin(A;A)) 

If we substitute the Taylor series expansion for the sine term this gives 

n' = |f(2A:A-(^ + 0(A:A)5) 
= aikn{l - ^ + 0{kA)^) 

where 0(/i;A)^ refers to fourth order terms. A similar treatment for the second 
differential gives the analytical result 

n" = -ak'^n (3.54) 

whereas the numerical expression is 



(kA) 



2 



n" = -ak\{l - + 0(A;A)5) (3.55) 

The error is given by the difference between the analytical expressions and the 
numerical result in both cases and depends on A^. This scheme is therefore 
second order accurate. 

We further note that the error also depends on the square of the wavevector 
k and will increase rapidly as k becomes comparable with A. A large value of k 
is equivalent to a minority carrier density which is varying rapidly in space. We 
must therefore be careful of trusting the results if the composition changes very 
rapidly over a small distance, and this design is not suitable for neutral regions 
with sharp discontinuities in any of the material parameters. This is unlikely to 
cause problems since sharp interfaces are generally avoided in the neutral regions 
because of the increased interface recombination they cause. 
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3.8.3 Comparisons of Analytical and Numerical Solutions 

In this section we brifely consider the stabihty of the numerical solution. We 
then continue to compare numerical and analytical calculation methods in order 
to test the numerical method. 

Stability of Boundary Value Problems 

Unlike initial value problems, boundary value calculations are not prone to nu- 
merical instability. This is because the solution is constrained at the two end 
points. 

However, rounding errors in the manipulation of large matrices can grow and 
swamp the real solution. However, this is rarely a problem in sparse matrices 
such as ours where the number of operations is relatively low. Furthermore, 
Potter [30] states that the use of double precision can circumvent this problem 
for medium size matrices where is of the order of some hundreds. 

Furthermore, Gaussian elimination with back-substitution is more stable than 
many more complex methods [29j . 

Zero Field 

In this section the results of the numerical model are compared to those of the 
analytical model in two cases where this is possible. 

Figure 13.21 shows a comparison of the internal p layer QE spectra predicted 
by analytical and numerical methods for an typical highly doped p layer which 
is 0.15/im thick and has a composition of 30% aluminium. The calculation uses 
a grid spacing of 400 points for the p layer. 

We note that the numerical method slightly underestimates the QE at short 
wavelengths. This is a result of increased numerical inaccuracy due to increased 
position dependence of both generation rate and minority carrier concentration. 
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Figure 3.2: Numerical and analytical internal QE plots for an ungraded 30% aluminium 
QWSC 




Figure 3.3: Numerical and analytical QE plots for an ungraded 30% aluminium QWSC with 
electric field 
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Cells With Fields 

No analytical method exists for a graded p layer, making precise evaluation of 
the performance of the numerical model difficult. As a first approximation, we 
compare results for the hypothetical case of ungraded cells with an effective field 
in the neutral regions. 

Figure 133] compares the analytical and numerical results for a QWSC identical 
to that of the previous subsection, but with an effective field of —lO^Vm^^ in 
the p layer. The results are again very similar, with the numerical model again 
underestimating the QE slightly at short wavelengths. 

The error in the numerical case is marginally greater than in the ungraded 
case, but still amounts only about one percent in the predicted QE at short 
wavelengths. It is again weakly sensitive to the grid spacing. 

Conclusions 

For the general case of an graded structure, we have no analytical model model 
for comparison. Analysis of the closest approximation however shows an error 
which is about a a fifth of the experimental uncertainty at short wavelengths. The 
model furthermore is relatively insensitive to the number of grid-points used. 

We conclude that the principal source of modelling error lies in the physi- 
cal assumptions underlying the model and in the values of physical parameters, 
rather than in the numerical method. This is investigated further by comparison 
with experimental data in chapter [HI 

3.9 Conclusion 

In this chapter we have described a method of calculating the light intensity for a 
solar cell with varying optical parameters and with or without multiple internal 
refiections for a normally incident spectrum. 

Having calculated the generation rate from the light intensity, we have looked 
at methods of calculating the resulting photocurrent. These methods fall into 
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two broad categories. The first deals with ungraded cells where we find that an 
analytical solution to the minority carrier concentration exists. The QE in this 
case is only weakly dependent on the diffusion coefficient for low recombination 
velocity. We have found that is it possible to estimate the diffusion length by 
modelling the QE in cells subject to this limit. 

Brief examination of the abstract case of an ungraded solar cell with effective 
fields in the neutral layers shows that such cells are significantly less sensitive to 
surface recombination. 

This useful property is also true for compositionally varying cells. We have 
described a numerical scheme for calculating the QE of such structures. The sec- 
ond case, which deals with graded cells, requires a numerical solution method and 
a good understanding of the dependence of transport parameters on aluminium 
content. 



Chapter 4 

Model Parameters and Program 
Design 

In this chapter we describe the way a QWSC is represented by the model. We then 
identify the main adjustable parameters used to fit experimental data. We give a 
review of the current state of knowledge of these and non- adjust able parameters 
used by the model. 

4.1 Program design 

The FORTRAN program SOLAL closely follows the progression outlined in sec- 
tion 131 The program design is modular and menu driven, and interfaces with the 
PAW graphics program. It calculates both analytical and numerical QE spectra. 

The main design parameters are summarised in table 14.11 The anti-reflection 
(AR) coat is represented by a wavelength dependent reflectivity, which may be 
arbitrarily specified by the user. Three standard reflectivities are provided. These 
are a free GaAs surface, and two standard reflectivities based on data from Eind- 
hoven and Sheffield. The latter two will be discussed in chapter [71 
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Surface 


GaAs reflectivity 

or SiN reflectivity 
or MgF:ZnS reflectivity 
or Measured reflectivity data 


Non-contributing 
surface layers 


Layer thicknesses 


Linear compositional profiles of the layers 


Number of layers 


p and n layers 


Layer thicknesses 


Arbitrary composition profile 


Constant doping 


Bulk i region 


Number of strips 


Thickness of strips 


Linear compositional profiles of strips 


MQW i region 


Total i region width 


Well width 


Barrier composition 


Number of wells 



Table 4.1: List of the main QWSC growth and design parameters. 



Sample structure 

The solar cells are represented as a series of AlxGai_xAs layers which may have 
linearly varying aluminium composition. These are split into surface layers which 
do not contribute to the QE, and contributing p, n, and i layers. Most of our 
designs have a single surface layer, which is the AlxGai_xAs window. 

The p layer can be subdivided into several strips, each with a linearly varying 
composition. The only recombination velocity included in the calculation is that 
between the AR coat or GaAs surface and the p layer. The doping however 
is assumed constant. Depletion widths are included in the calculation, and the 
effective aluminium fraction at the edge of the depletion width is readjusted in 
graded structures to match that at the edge of the calculated depletion width. 
The treatment for the n layer is identical, with the difference that the surface 
or interface recombination velocity is replaced by a back surface recombination 
velocity. 

The i layer may be of two types. The first is a conventional i layer, which may 
be composed of a number of linearly graded strips. A constant background doping 
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Ungraded p 


Sn \ Recombination parameter 

Dn ) "S = Sn/ Dn 


Diffusion length L„ 


Graded p 


Sn ] 
Dn{x = 0) J 


Recombination parameter 

■ S = Sn/Dn{x = Q) 


Diffusion length Ln{x) 


Diffusivity Dn{x) 



Table 4.2: Transport parameters for the p layer required by the model 



level is assumed throughout. The second type of i layer is one which contains 
quantum wells. This layer may not be compositionally graded except for the 
sharp interfaces between quantum wells and barriers. Details of the quantum 
well modelling parameters have been supplied by Paxman |49j. The parameters 
in this case are the total i region width, the barrier composition, the number of 
wells and the well thickness. 

Minority carrier transport in neutral layers 

The minority carrier transport parameters required to calculate the QE of the 
neutral p and n sections vary according to the type of structure being considered. 
They are summarised for graded and ungraded structures in table 14.21 

For ungraded structures, the ratio of surface recombination velocity and dif- 
fusivity are specified for each layer, together with the minority carrier diffusion 
length. 

Graded structures require a parametrisation of diffusion length and diffusivity 
as a function of position x. 

These parameters are reviewed in the following sections. The parametrisations 
used are reviewed in chapter [5l 

4.2 Design and Growth Parameters 

In this section we review the parameters summarised in tables WA\ and W?2\ These 
parameters fall into two categories, the first of which are functions of cell design 
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and growth. Examples are the dimensions of the cell layers, their composition 
and doping levels. Knowledge of these parameters is usually good because the 
QWSC fabrication process uses highly accurate epitaxial methods. 

The second category contains parameters which are a function of growth 
method, growth run and processing. Examples include surface recombination, 
transport parameters and optical properties. These parameters are largely taken 
from existing data in the literature. 

Tables H]T] and H]2] show the most important modelling parameters for a graded 
QWSC with an anti reflection (AR) coat. Dimensions such p,i and n layer thick- 
nesses, and doping levels are generally well known. 

The quantum well width is frequently found to depart from the design spec- 
ification by a few percent. However, this does not affect the overall electrical 
characteristics of the QWSC and simply restricts or extends the QW absorption 
edge in the infrared. Such deviations from the specifications are discussed further 
with respect to modelling in chapter |8l 

The background doping level Nbg in the intrinsic region, however, is fre- 
quently subject to large uncertainty. This is largely because residual impurity 
levels in the solid or gas sources are poorly known. In addition, the net back- 
ground doping level is a function of reciprocal compensation by donor and ac- 
ceptor impurities. The final background doping level may also be affected by 
unpredictable passivation effects during processing. 

As detailed in chapter |2l a large value of Nbg degrades the performance of 
the QWSC because the built-in electric field drops to zero in the intrinsic region. 

The QE measurements however are taken when the photocurrent is max- 
imised. In this case, the only effect of Nbg on the photocurrent is to influence 
the width of the depletion layers. A more detailed discussion is given with ex- 
perimental data in chapter [71 
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AR coat 


~70A SiN, measured reflectivity 


Window layer 


300A Alo.67Gao.33As 


p layer 


0.15/im Alo.4Gao.6As 
graded to Alo.2Gao sAs 
Carbon doping 1.34 x lO^'^cm"^ 


i region 


30 GaAs wells of width 87A 
60A barriers, composition Alo.2Gao.8As 
Well absorption ~1% per level per well 


n region 


0.6/im Alo.2Gao.8As 
Si doping 6 x lO^^cm"'^ 


Metallic back mirror 


Reflectivity ~95% 
across entire wavelength range 



Table 4.3: Typical design parameter set for an anti- reflection coated 20% Al 30 well QWSC 
with a graded p layer and a back mirror 



4.3 Optical Parameters 
4.3.1 Incident Flux 

The QE calculation is independent of the incident flux. Estimated values of the 
short circuit current given in subsequent chapters are given for the AMI. 5 solar 
flux which was described in section 12.41 



4.3.2 Absorption Coefficient 

Connolly [31] has shown that the QE is very sensitive to the absorption coefficient 
a, and has proposed an absorption model based on accurate GaAs tabulations 
published by Aspnes [32] • The absorption model uses a non linear wavelength 
dependent energy shift which depends on the aluminium composition x to fit the 
GaAs absorption to data for AlxGai_xAs published in [33]. This approach only 
gives the absorption of AlxGai_xAs for direct aluminium compositions. 

Subsequent work in [19] fits an exponential tail below the direct absorption 
threshold for indirect aluminium fractions to account for the weaker indirect ab- 
sorption. Since the indirect and direct absorption routines were initially slightly 
ill matched, this work introduces a further continuous transition between the two 
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parametrisations. 

Connolly [31] has also indicated that the absorption coefficient near the band 
edge is subject to significant uncertainty for aluminium fractions in the region of 
30%. This is seen both in the experimental data tabulations and in a range of 
models which are reviewed in the reference. 

We conclude that the absorption coefficient remains a source of systematic 
error near the band edge. Since, however, the absorption model uses the most 
recent tabulated absorption data available, it is difficult to reduce this source of 
uncertainty. 

4.3.3 Reflectivity 

Given the very small size of typical QWSC mesa-photodiode test devices, direct 
reflectivity measurements are difficult. In order to make such measurements 
possible, pieces of wafer were exposed to all processing steps at the same time 
as the material which was used to fabricate devices, with the exception of mesa- 
etching and metallisation. This is discussed in more detail in section [7] 

These pieces of AR coated wafer are not available for devices processed before 
1995. For such samples, an average over previous processing runs is used in the 
modelling. 

The GaAs reflectivity was obtained from published data [32]. Measurements 
on GaAs wafers has produced close agreement. 

Further reflectivity for the double-layer MgF:SiN coating deposited and char- 
acterised in Eindhoven Technical University was available for devices processed 
in Eindhoven. 

The reflectivity of the metallic back- mirror is given by Hass [34j- It is found 
to be very nearly constant across the wavelength range of interest for gold and 
aluminium. The front surface reflectivity was modelled as a fraction of the mea- 
sured surface reflectivity. The corresponding fraction is the prime parameter 
used in modelling the amplitude of Fabry-Perot peaks in mirror-backed devices, 
and varies from about 70% to 100% of the measured reflectivity. This factor is 
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described in more detail in chapter |8l 

The front surface phase change tpf for internal reflection in a dielectric material 
at an interface with air is given by Macleod [32] as zero. The phase change is 
TT if the interface is with a material with a higher refractive index. Work by 
Whitehead [36] on multiple quantum well optical modulators indicates that the 
phase change ipb at the metallic back mirror is variable and device dependent. We 
shall see in chapter [8] that different devices processed from the same wafer may 
have different modelled values of ipb- More recent discussions with the author [37] 
have suggested that this quantity can be used as a fitting parameter to adjust 
the frequency of the Fabry-Perot oscillations seen in the QE of mirror backed 
devices. 

4.3.4 Refractive index 

Reliable ellipsometric measurements of the refractive index of AlxGai_xAs as 
a function of energy have been tabulated by Aspnes [32] for eleven aluminium 
compositions for aluminium mole fractions ranging from to 1. For intermediate 
energies and compositions, a linear interpolation is used. 

4.4 Transport parameters 

We have seen that good agreement exists on the optical materials parameters of 
AlxGai_xAs, and that consistent values are available in the literature 

Most of the minority transport parameters however are sensitive to device 
history. Growth conditions between different runs on the same machine may 
differ significantly due to relatively small variations in the composition of the 
gas or solid sources used in the epitaxial growth. Furthermore, post-growth 
processing, involving device exposure to etches, plasmas and heat treatments, 
may also affect transport. Examples include accidental passivation, diffusion of 
impurities and a modification of the recombination velocity. Finally, a number 
of these parameters are dependent on the device structure. The following section 
explores current knowledge of these parameters in the literature. 
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4.4.1 Band gap Parametrisations and the Quasi-Electric 
Field 

We saw in chapter |2] that an externally applied bias in a p - i - n structure is 
entirely dropped across the intrinsic region for moderate applied bias where the 
majority carrier profile is not significantly affected by the external voltage. That 
is, no externally applied bias is dropped across the neutral sections of the n and 
p layers. 

The effective electric field acting on minority carriers in n and p type graded 
structures is examined by Sutherland [22], Hutchby [21] and is reviewed by 
Hamaker [21]. It is generally referred to as the quasi field. This electric field is 
sometimes taken to be proportional to the gradient of the electron affinity as in 
Sutherland [22]. Other workers however express it as the gradient of the bandgap. 
The review by Hamaker [21] explains the origin of the field and expresses it in 
terms of doping, bandgap and Fermi level gradients. In samples designed in the 
present work, the nominal doping is constant. We will not therefore consider 
fields due to doping gradients. 

Sassi [23] has estimated the field due to the Fermi level gradient as less than 
2% of the field due to the bandgap gradient. Sassi therefore neglects the gradient 
of the Fermi level. Taking the example of the minority electron quasi field in 
p type material, this gives 

^' = 

where Eg is in electronvolts and x the position. This assumption of low Fermi 
level gradient is especially true in heavily doped material such as we use in our 
solar cell designs. It is the most commonly found parametrisation for the quasi- 
field acting on minority carriers, and is used for example in work by Konagai [20] , 
Hutchby [21] and Sassi [23]. We shall therefore use this method in this study. 
The corresponding hole quasi field in electronvolts for n type material is 

„ dE„ , , , 
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We shall not consider graded n layers in the AlxGai_xAs system, and will revert 
to the nomenclature of chapter [31 which refers to E^. as E. 

An added complication is the transition from direct to indirect material in Alx- 
Gai-xAs. The review article by Pollack [39] put this transition at a composition 
in the region of 45% aluminium. 

Two expressions are given for Eg in the literature. We write these in terms 
of a composition dependent step function G for convenience, which has value 1 
below the direct-indirect transition, and above. The bandgap is derived from 
expressions by Casey and Panish |44j for the direct regime, and Hutchby [2T] for 
the indirect material. 

Eg = 6(1.424 + 1.247X) + (1 - e)(1.92 + 0.17X + O.OTX^) (4.3) 
4.4.2 Minority carrier lifetime 

Although the model outlined in section 13.41 does not explicitly require knowledge 
of the minority carrier lifetime r, this parameter is mentioned here for complete- 
ness. It is related to the minority carrier diffusivity and the diffusion length 



for electrons. A similar expression holds for holes. Minority carrier lifetimes 
in GaAs/AlxGai_xAs have been measured for a range of aluminium composi- 
tions. The most common methods are time resolved photo-luminescence (TRPL) 
(Ahrenkiel [H]) and zero-field time-of- flight (ZFTOF) (Zarem [12]) techniques. 
This data is summarised and discussed in the review article by Ahrenkiel [15] . Un- 
fortunately, little data is available for compositions above about 40% aluminium. 
The large scatter observed in experimental data indicates that the lifetime is 
dominated by defect-related Shockley-Hall-Read recombination processes. It is 
therefore sensitive to material quality and according to Ahrenkiel may vary by 
orders of magnitude compared with the published values in the review article. 



by 




(4.4) 
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4.4.3 Mobility and Diffusion Coefficient 

We have seen that the model does now exphcitly depend on the magnitude of 
the diffusivity. The functional form of this parameter however is important in 
determining the QE of graded structures. 

Values of the diffusion coefficients Dp and Dn for minority holes and electrons 
are not available in the literature, and, like the diffusion length values, must be 
obtained indirectly. 

The review article by Ahrenkiel |15] states that the diffusivity is most conve- 
niently found by extrapolation from the electron and hole majority Hall mobilities 
lie and fih in the low field, low injection limit where the Einstein relationship is 
valid. The two quantities are then approximately related as follows 

n ~ hBiEK (holes) 
P 1 ^ ^ (4.5) 

Dn ~ ^^^J^ (electrons) 

Low field hole and electron mobilities are supplied in the review by Walukiewicz 
|18] . A wide spread is observed, indicating that the mobility is material de- 
pendent. Both electron and hole mobilities decrease sharply from GaAs to 
Alo.5Gao.5As. For aluminium compositions greater than about 50%, both mo- 
bilities appear essentially constant. 

This behaviour is essentially reproduced by a parametrisation of the diffu- 
sivity as a function of composition and doping presented by Hamaker [25]. The 
Hamaker paper does not give details, but states that the parametrisation is based 
on experimental data. We use this as a guide to diffusivity values. More detail 
concerning this point is provided in chapter 13 

4.4.4 Diffusion length 

The diffusion lengths Ln and Lp in GaAs/AlxGai_xAs are sensitive to growth 
method, device structure and aluminium composition. As we shall see in chapter 
[5l they are significant QE fitting parameters. We review here some of the expected 
values in AlxGai_xAs. Specific values of this parameter will be referred to with 
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respect to modelling in chapter |8l 

Ahrenkiel recommends extracting L for both minority holes and electrons 
from the minority carrier lifetime and the diffusion coefficient discussed above 
through the relation 



for minority electrons. A similar expression holds for minority holes. These 
values can only provide a rough starting point for modelling, since the diffusion 
length is strongly device dependent. 

Ahrenkiel reports that the diffusion length dependence on aluminium fraction 
is similar to that of the mobility. It generally decreases with increasing aluminium 
fraction over the direct region, but is found to increase near the direct /indirect 
transition region at about 40% aluminium. It remains approximately constant 
for indirect material, at a value which is however only a fraction of the GaAs 
diffusion length. 

Hamaker [21] gives a parametrisation similar to that which was discussed for 
the diffusivity in the previous section. We shall compare values predicted by 
Hamaker with values estimated from the modelling of QE spectra in chapter |8l 



Griinbaum |46] has measured minority hole diffusion lengths in some of our 

AlxGai_xAs material using Electron Beam Induced Current techniques. This 
technique however was only applied to one sample relevant to this study. It is 
unfortunately currently incapable of the resolution required to measure the more 
significant minority electron diffusion length Moreover, the predicted Lp 
values are rather large. 

4.4.5 Recombination velocity 

Chapter [3] showed that surface or interface recombination velocities in our mate- 
rial cannot be separated from the diffusivity, which was reviewed earlier in this 
chapter. Modelled values of the surface recombination parameter S however in- 
dicate whether or not our values are in qualitative agreement with those reported 
in the literature. 




(4.6) 
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An unfortunate result of any heteroj unction in a semiconductor structure is 
the presence of interface recombination centres. These may result for a variety of 
reasons. For heterojunctions these include lattice mismatch and abrupt changes 
in growth conditions. 

Two more cases occur. The first is the junction between a semiconductor and 
air. Recombination velocities in this case are particularly high. Measurements 
presented by Fonash [19j quote values up to three orders of magnitude greater 
than for heterojunctions. The third case is that of a semiconductor and an anti- 
reflection coating, such as deposited SiN. 

Measurements are carried out by time resolved photoluminescence (TRPL) 
studies of double heteroj unction devices. The active region is grown with a lower 
bandgap, and its thickness varied. The variation of the PL signal decay time as 
a function of active layer thickness can then give an indirect measurement of the 
interface recombination velocity. This method has been used by 'tHooft and Van 
Opdorp [17j. Results are included in the review article by Timmons |40j . 

Devices which have suffered surface damage as a result of handling or pro- 
cessing may however deviate from this picture. Furthermore, Fonash [12] has 
indicated that recombination at the contacts increases the effective surface re- 
combination. This effect can be reduced by the use of appropriate contacts. 
Selective ohmic contacts, for example, present a potential barrier to minority 
carriers, and inhibit minority carrier recombination at the contacts. They are, 
however, ohmic to majority carriers [T9] . 

Theoretical studies however such as those reviewed in chapter E] usually use 
values which are one to two orders of magnitude greater than these. The largest 
reported value for interface recombination in AlxGai_xAs is quoted by Timmons 
PU] as 3 X lO^cm/s. Surface recombination velocities of 4 x lO^cm/s have been 
observed in Alo.08Gao.92 As material exposed to the air. 

Details concerning the recombination velocity and the way this parameter is 
combined with the diffusivity are given in chapter [5] which deals with modelling 
techniques. 
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4.4.6 Discussion 

In this chapter we have briefly described the method used to design the modelhng 
program, before discussing the parameters required to model the QE of real 
solar cells. We have seen that these parameters can be divided into two broad 
categories. The first are optical, compositional and geometrical parameters and 
are resonably well known. 

The second category involves the transport parameters. These are gener- 
ally dependent on sample growth and device processing, and are poorly defined. 
Values quoted in the literature are variable and show a large spread. However, 
experimental and parametrised values quoted in the literature provide bounds 
within which we can expect our values to lie. 

The treatment of the modelling techniques associated with using mutually 
consistent values requires simultaneous examination of the modelling results. 
These considerations are desribed in chapters |5] and [H 



Chapter 5 

QE Modelling and Optimisation 



The reliability of the analytical method has been established in previous work by 
modelling published data by Bernd Braun |18]. The first section of this chapter 
summarises the main results of this work. 

We go on to investigate a number of theoretical predictions made by the 
model for a number of promising structures. All results apply to the QWSC 
design shown in table 15.11 and to p - i - n control samples which are identical 
except for the lack of wells in the intrinsic region. 

The optimisation is broadly separated into two cases according to incident 
photon energy. The first is the case of photons with energies above the bulk 
bandgap Ei,2, which corresponds to aluminium fraction X{,2 (see table 15. 

The QE in this range is mainly determined by the minority carrier transport 
efficiency of the p layer. Graded structures are further characterised by the 
aluminium fraction Xbi at the top of the grade, and its associated bandgap Ebi. 

The separate optimisation of the quantum well QE applies to incident photons 
with energies between the quantum well bandgap Ea and the lower bulk bandgap 
Eb2 and its associated aluminium composition Xb2- This optimisation is mainly 
dependent on the absorptivity of the quantum wells and the light levels in the i 
region. 
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Layer 
type 
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i 
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i buffer 


400A 


i 


X = Xbi 


n layer 


0.6/im 


n 


X = Xbi 


GaAs substrate 






x = 



Table 5.1: Basic QWSC growth menu 



5.1 Ungraded Cells 

5.1.1 Previous Work 
Photocurrent and Voltage Enhancement 

Early work presented by Mark Paxman et al. in reference [IH] showed that 
the short circuit current of a QWSC is greatly enhanced over that of a control 
structure which is identical, but with the well material replaced by an equal 
thickness of barrier material. 

Further work showed that the Voc for the QWSC drops relative to that of the 
control sample, but remains greater than the Voc which would be expected from 
a cell with the quantum well effective bandgap Ea. 

One of the samples with fifty wells showed a Voc of 1.07V, which is greater 
than that of the worlds' best GaAs cell (reference [H]). The modelled light IV 
curves and QE spectra for this sample are illustrated in figure 15.11 The large 
increase in J^c more than compensates for the decrease in Voc leading to an overall 
efficiency enhancement in the QWSC relative to the p - i - n control. 

A set of identical ungraded samples were designed with the help of the opti- 
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Fi gurc 5.1: Modelled a) light IV curves and b) QE spectra for a 50 well QWSC and p - i - n 
control pair based on a successful device. 
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misation program we described in chapter [3] and grown by MBE. They feature 
different values of Xh2 but identical geometries. We concluded that the Voc in- 
creases with Xh2. Unfortunately, the material quality and QE decrease with 
increasing Xb2. We concluded that a value of Xb2 in the region of 30% is a rea- 
sonable compromise between the requirements of good material quality and high 

Voc. 

Anti Reflection (AR) Coating 

The reflectivity of the standard single layer SiN AR coat has been optimised in 
reference [19]. Further reductions in reflectivity across the wavelength range of 
interest are possible with more sophisticated multilayer AR coats. 

This work, however, is concerned mainly with the internal QE of QWSC 
structures, and has used the recommended SiN AR coat of thickness ~ 70nm 
throughout. 

Short Wavelength Optimisation 

Previous work has recommended a number of optimisations which may be carried 
out. This work concluded that thinning the p layer increases the minority carrier 
collection efficiency, particularly at short wavelengths. It also transmits more 
light to the i layer. Both effects enhance the QE. This method, however, is 
limited by the increase in series resistance as the p layer is thinned. This is 
particularly damaging for concentrator applications. 

The addition of a window layer has been found to reduce recombination at the 
front of the p layer. This is due both to the replacement of a surface recombination 
velocity by a lower interface recombination between the p layer and the window, 
and the reflection of minority carriers away from the surface by the bandgap 
discontinuity at the interface. 

Although this method successfully reduces the recombination velocity, some 
of this advantage is lost because even thin window layers absorb a significant 
fraction of incident light at short wavelengths. The previous work has shown 
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that the window layer is very inefficient at converting this absorbed radiation to 
photocurrent, and essentially contributes no photocurrent. 

Use of a very high aluminium fraction for the window layer can reduce the 
window absorptivity, but introduces new difficulties. This is due to an increased 
recombination velocity for high window aluminium fraction, caused both by an 
increased lattice mis-match at the interface and lower window material quality. 

The techniques above have been successfully demonstrated in reference |19] . 
However, we conclude that there is still much room for improvement at short 
wavelengths, which requires more complex structures. 

Long Wavelength Optimisation 

As shown indicated by Paxman |19], the long wavelength QE can be increased 
simply by introducing more wells in the i region. This can be done by reducing 
the barrier thickness. This approach however is expected to reduce the Voc of the 
cell. 

Increasing the width of the i layer is an alternative route. We saw however in 
section 12.3.21 that this leads to a reduced QE in forward bias and a degradation 
of the fill factor FF. Furthermore, a more fundamental limit is imposed by the 
increase in dark current for wider i layers. 

Previous work resulted in a QWSC with fifty wells, which nevertheless failed 
to absorb approximately half the incident light at energies below Ei,2. Failing 
material with lower background doping, more complex techniques are required 
to increase the QE in the long wavelength range. 

Conclusions 

We have seen that the short circuit current enhancement techniques presented 
by Paxman |19] are limited by a range of physical and materials considerations. 

Short wavelength efficiency enhancements in ungraded GaAs/AlxGai_xAs so- 
lar cells are limited by series resistance, window absorptivity and poor material 
quality at high aluminium fractions. 
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At longer wavelengths, the background doping sets a limit on the number of 
wells which may be inserted in the structure, and hence the total absorptivity of 
the wells. 

In both wavelength ranges therefore, there is potential for optimisation by 
improving p layer transport and by increasing the absorptivity of the i layer. 

5.2 Further QE Enhancements Above the Bar- 
rier Bandgap 

Chapter [3] shows that a variation of material composition with position can sig- 
nificantly alter the optical and transport properties of the QWSC This section 
investigates the effect of such changes on light intensity and generation rate as a 
function of position and incident light wavelength. The consequences for minor- 
ity transport are covered before reviewing the minority carrier density and QE 
optimisation of graded structures. 

5.2.1 Light Absorption and Minority Carrier Generation 
Light Intensity 

In chapters [3] and H] we saw that AlxGai_xAs is direct for aluminium fractions 
up to approximately 45%. Direct materials are characterised by high absorption 
coefficients. Unlike indirect materials such as silicon, direct materials can absorb 
most of the incident light spectrum in a relatively short depth. However, given 
the poor materials characteristics of AlxGai_xAs, it is desirable to increase light 
levels deeper in the cell so that minority carriers are generated preferentially near 
or in the depleted layers where they may be collected by the built in potential. 
This may be achieved by compositional grading. 

The relative light intensity for typical ungraded and graded p - i - n cells is 
shown for a range of wavelengths in figure 15.21 The intensity immediately starts 
to decrease in the ungraded cell for all wavelengths above the bandgap. At very 
short wavelengths, no light reaches the i region. 
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Figure 5.2: Relative light intensities in a a) ungraded p - i - n cell and b) a p - i - n cell with 
a p layer graded from Xbi = 67% to Xb2 = 30% 
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Figure 5.3: Fractional increase in light intensity in a p - i - n cell with a p layer graded from 
Xbi = 67% to Xb2 = 30% compared to an ungraded cell 
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In the graded sample for long wavelengths such as 600nm, absorption in the 
p layer only becomes significant at a depth of about 70nm, at which point the 
bandgap decreases below the corresponding absorption threshold. A constant 
plateau in light intensity is seen from zero to 70nm because of the lack of absorp- 
tion. 

Figure 15.31 gives the difference in light levels between ungraded and graded 
cells, expressed as a percentage of the light level in the ungraded cell. Higher 
light levels are seen for all wavelengths in the graded design. The gain in intensity 
increases with depth until the junction with the i layer at a depth of 150nm. For 
greater depths, no further increase is seen since the absorption coefficients of the 
i and n layer are identical for the two cell designs used in this example. 

The gain in light level is strongest at short wavelengths because the difference 
in the absorptivities of the graded and ungraded p layers is greatest for high 
energies. The calculations of intensity at 400nm represented in figure 15.21 show 
that ~ 95% of the incident light is absorbed in the first 40nm of an ungraded 
cell, whereas ~ 85% is absorbed in the same depth for the graded structure. 
Whilst the light intensity in the i and n layers appears much higher in the graded 
structure relative to the ungraded cell, the absolute light level remains very low. 
The absolute light intensity in the i and n layers in the short wavelength range 
therefore remains very low. 

Improvement in the short wavelength range will be due to an improved p 
layer response due to a better distribution of light, as we shall see in subsequent 
sections. 

For wavelengths between approximately 450nm and the bandedge, increases of 
between 10% and 40% are seen throughout the cell. This is particularly significant 
for the i layer, which converts absorbed light into photocurrent with nearly 100% 
efficiency. For energies very close to Ef,2, the intensity enhancement tends to zero 
since neither cell absorbs strongly. 
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Figure 5.4: Generation rates in an ungraded p - i - n solar cell for wavelengths above the 
bulk bandgap 
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Figure 5.5: Generation rates in a compositionally graded p - i - n solar for wavelengths above 
the bulk bandgap Eb2 
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Generation Rate 

Figures 15.41 and 15.51 show the generation rates for the structures of figure 15.31 as 
a function of depth and position. In the homogeneous cell, the generation at 
each wavelength decreases exponentially with increasing depth. The exponential 
decay however very nearly approximates a constant for wavelengths within about 
50nm of the bandedge. 

The graded case is considerably different. Considering the p layer first, we see 
that the generation rate at short wavelengths is decreased because of the lower 
absorptivity near the surface. 

At lower photon energies, the generation rate falls to zero as the high Al frac- 
tion material towards the front of the cell ceases to absorb light. The generation 
rate at 650nm for example only becomes significant at a depth of about 0.1/xm, 
which is over two thirds of the total p layer width. 

We conclude that the grade improves the generation profile both by shifting 
minority carrier generation rate in the p towards the i layer, and by increasing 
the generation rate in the latter. An increase in generation rate in the n region is 
also visible but is not expected to modify the n region QE considerably because 
the light levels in this layer remain low. The magnitude of the subsequent QE 
enhancements in all three regions of the solar cell is discussed below. 

5.2.2 Minority Carrier Transport 

In chapter [3] we saw that the QE of graded structures depends on the diffu- 
sion length and diffusivity of the neutral layers. Both of these parameters are 
composition dependent functions. 

Neither quantity can directly be measured in the QWSC structures. We are 
therefore obliged to make a number of assumptions regarding the functional form 
of both the diffusion length and the diffusivity. 
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n Layer Minority Carrier Transport 

In the structures we are considering, the n layer remains ungraded, and only 
constant values of Lp and Dp are required. Since the n region contribution in 
our case is small, we cannot estimate values of these parameters from QE mea- 
surements. We rely instead on the values reviewed in chapter |H Although this 
approach is inaccurate, the effect on the absolute QE is small because of the very 
low n region QE contribution in p - i - n AlxGai_xAs devices. 

p Layer Minority Carrier Diffusion Lengths 

The spatial variation of L„ and Dn is an important factor in determining the 
overall p layer QE. A number of assumptions must be made to compensate for 
the lack of direct measurement. 

Appendix |X] shows that the QE of ungraded p layers is directly proportional 
to Ln for low absorption and low Experimental measurements of the QE of 
ungraded devices with different values of Xb2 can therefore be used to build up a 
picture of the functional dependence of L„ on aluminium composition. In turn, 
this gives the spatial variation of L„ in compositionally graded structures. 

We make the further assumption that growth conditions are reproducible 
and the resulting table of L function of Xb2 can be generalised to graded 

structures. 

p Layer Diffusivity 

Chapter [3] showed that the QE of ungraded p layers depends on the ratio Sn/ Dn. 
Since the recombination velocity Sn is device dependent and cannot be measured 
directly in our structures, we cannot extract the functional form of Dn from QE 
spectra. We must therefore make a further assumption allowing us to parametrise 
the effect of the diffusivity in terms of a single parameter. 

Chapter [3] further showed that Dn affects the QE of graded structures in two 
ways. The first is similar to the ungraded case and involves the recombination 
velocity parameter S. It was also demonstrated however that S has a weaker 
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effect in graded structures because of tlie bandgap gradient effective field £f. 
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Figure 5.6: Effective electric field in per micron as a function of position for a p layer graded 
from Xbi = 67% to = 30% 

Dn also appears as a term D'^/Dn which we define as the diffusivity field X>/ 
in units of per metre: 



(5.1) 



This term expresses the current flow caused by minority carriers in regions of 
high mobility moving towards regions where minority carriers are less mobile. 
Although it is expressed here in terms of the differential with respect to distance 
VxDm it can be re-expressed in terms of the aluminium compositional gradient 



(5.2) 
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This form, which we will call the diffusivity field, has the advantage of being 
independent of the device geometry, and is dimensionless. The diffusivity in 
terms of these two forms is 



where X{x) is the aluminium composition at position x. The diffusivity field is 
related to Vf via the relation 



The diffusivity field defined this way has a value of about -4 if we use the Hamaker 
parametrisation of chapter HI for p type AlxGai_xAs doped at about 2x lO^^cm"'^. 

Since comparisons between samples with different geometries are more clearly 
apparent with this form, we shall use it in the modelling section. We continue to 
use the device dependent form Vf here, however, because it has the same units 
as the effective field Sf and is therefore better suited to this discussion. In terms 
of Pj, then, the total field St acting on the electron minority carriers is 



According to experimental data presented in reference [TS] (reviewed in chap- 
ter H]), Dn decreases rapidly with increasing aluminium fraction, before levelling 
off for indirect material. Figure \5M shows the resulting fields Vf,£f and £t where 
Sf is calculated as in chapter H] and Vf is found from the experimental data of 
reference [T8] . 

Since the diffusivity is higher in AlxGai_xAs with a low aluminium fraction, 
Vf is positive and, unfortunately, opposes the drift field £f due to the bandgap 
gradient. The figure shows that a significant decrease in St may be expected. 

The functional form of Dn therefore may have important consequences for the 
overall effect of a compositional grade on minority carrier transport efficiency. 



Dn{x) = Dn{0)exp{Vfx) 

= DniO)exp{V;^iX{x)-Xbi) 



(5.3) 




(5.4) 



St = Sf + Vf 



(5.5) 
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Since we cannot experimentally determine the compositional dependence of Vf, 
we must find a parametrisation which can reduce this poorly defined function to 
a single fitting parameter. 

An exponential parametrisation for D„ suggested by Hamaker [2l] was re- 
viewed in chapter HI The Hamaker parametrisation reproduces the functional 
form of the experimental data given in reference pB]. For the example shown in 
figure [5^ the Hamaker value is Vf ~ 10/im~^, which agrees qualitatively with 
the published experimental data. Because of material variability, however, we do 
not expect the Hamaker parametrisation to yield accurate values of Vf for our 
material. 

In the general case of an exponential parametrisation for Dn, is a position 
independent constant for linear compositional gradients. Moreover, it is indepen- 
dent of the magnitude of D„. This means that the magnitude of the diffusivity 
influences the QE of a graded structure only through the surface boundary pa- 
rameter S. Since, as we have mentioned ealier in this section, S has a weak effect 
in graded structures, we are left with a single principal fitting parameter which 
is the diffusivity field Df. 

We will therefore assume an exponential form for Dn along the lines suggested 
by Hamaker but will treat the diffusivity field as a position independent fitting 
parameter. 

Summary of the Transport Model 

Modelling of the QE from the neutral regions is therefore subject to a number 
of assumptions which we will review here. 

The first is that the n region contribution is small. We only treat ungraded 
n layers, and rely on published data for transport parameters in this case. 

For ungraded p layers, we have seen that parameters are well known except 
for Lrf, Dn and Sn- Because the diffusion length in our structures is smaller 
than the p layer width be estimated from the QE close to the bulk 

bandedge at energy Eb2. 

The parameters Sn and Dn in ungraded structures can be separated as a 
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single surface recombination parameter S (see chapter [3]) which is used as the 
sole fitting parameter in the short wavelength regime. 

For compositionally graded structures, we make the assumption that the epi- 
taxial material is of reproducible quality. The values of L„ derived from modelling 
the QE of ungraded structures with different values of Xb2 are used to build up a 
table of the compositional dependence of L„ . Under the aforementioned assump- 
tion, this provides the position dependence of L„ for the grandad structures. 

We finally make the assumption that the functional dependence of Dn on 
aluminium fraction is exponential. The diffusivity field I'/ is therefore a constant, 
which is independent of the magnitude of the parametrisation Dn{x). Since 
chapter [3] has shown that the QE of graded structures is only weakly dependent 
on the surface recombination parameter 5, we are left with a single dominant 
fitting parameter which is "Dj. 

The number of assumptions involved in this modelling prevents us from ex- 
pecting quantitative predictions of the QE of graded structures. We expect, 
however, to qualitatively describe the dependence QE enhancement factors due 
on the one hand to better light distribution and on the other, to improved mi- 
nority carrier transport in graded solar cell structures. 

5.2.3 Minority Carrier Concentration 

Figures 15.71 and 15.81 show the minority carrier concentrations for an ungraded 
p layer and a cell graded from 30% to 40%. The concentrations are given in 
arbitrary units but are subject to the same normalisation factor. 

We see that the minority carrier concentration in the graded case decreases in 
magnitude by a factor of about two, and is shifted towards the depletion region. 
The decrease in concentration is due partly to decreased absorption and partly 
to the modifications in minority carrier transport efficiency caused by the grade. 
The effective field in particular sweeps carriers to the junction. Increased bulk 
recombination due to the higher aluminium fraction also plays a role. 

Figure 15.91 shows the case of a p layer graded from 30% to 67%. In this 
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Figure 5.7: Minority carrier concentration in an ungraded p layer for a range of wavelengths. 
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Figure 5.8: Minority carrier concentration in a p layer graded from 30% to 40% Al for a 
range of wavelengths. 
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Figure 5.9: Minority carrier concentration in a p layer graded from 30% to 67% Al for a 
range of wavelengths. 
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case, the two separate field regions are clearly reflected in the minority carrier 
concentration. The lower direct transition region near the depletion layer has a 
lower minority carrier concentration in this case because the higher field sweeps 
carriers towards the depletion layer more efficiently than the lower field in the 
indirect region near the surface. 

The overall minority carrier concentration is again lower than the previous 
two cases we have looked at for the same reasons. The change is greater in this 
case because of the increased compositional gradient. 

5.2.4 QE and Short Circuit Current 

The quantum efficiency for the ungraded and graded cells described in the previ- 
ous sections is shown in figure 15.101 The breakdown of percentage increase in p, 
i, n and total QE is given in figure 15.111 We note that the p layer contribution 
is increased at short wavelengths but decreased near the bulk band edge. 

The improvement at very short wavelengths is the result of increased carrier 
collection efficiency in the p layer. This is due to improved minority carrier 
transport and the shifting of the minority carrier generation profile towards the 
i region. The deterioration near the band edge however is due to reduced light 
absorption in the p layer at low photon energies. 

Figure 15.12b shows a prediction of short circuit current as a function of the 
topmost aluminium fraction in the p layer. Such calculations are subject to large 
uncertainty however, particularly given the large distribution of possible diffu- 
sivity gradients. The calculation of figure 15.12b uses a gradient which is about 
half that of the published data of chapter |4] which as we shall see in chapter |8] is 
consistent with experimental QE spectra. Although the model cannot indicate 
the magnitude of the diffusivity, this deviation of a factor of two in gradient is 
well within the margins expected. 

With the uncertainties due to in mind, we note that Jgc initially rises 
sharply as the effects of surface recombination become less important. For higher 
grades, the gain in i region contribution is increasingly balanced by a reduction of 
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Figure 5.10: QE of a) a 30% p - i - n and b) a p - i - n with a p layer graded from Xh\ = 67% 
to Xhi = 30% 
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Figure 5.11: Percentage inerease in QE for p, i, n layers and total enhancement for a graded 
p - i - n device relative to an ungraded control. 
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the p layer current. In this example the QE reaches a clear optimum. Relatively 
small changes to the diffusivity can have significant effects on the behaviour of 
the QE at high aluminium fractions. 

The common feature of such theoretical studies is the increase of QE up to 
Xbi ~ 50% ± 5%. Thereafter, more optimistic parametrisations of V than that 
used in the figure show a QE which continues to increase monotonically as a 
function of X^i, albeit at a slower rate. 

We conclude that a successful grade in a cell with Xb2 ~ 30% requires a grade 
of approximately Xhi > 40%. A higher grade is not expected to significantly 
affect the overall short circuit current, and may in fact improve it, depending on 
the depth dependence of the diffusivity. This is partly due to the fact that two 
enhancement mechanisms are present, and that even devices with disappointing 
transport parameter profiles will benifit from enhanced i region QE because light 
intensity QE enhancement mechanism is independent of transport and is greater 
for higher grades. 

The choice of the higher aluminium fraction is not therefore critical from a 
modelling point of view. It is possible to favour one enhancement mechanism 
over another without seriously affecting the overall short circuit current. 

In real devices, however, this must be offset against the technological problems 
involved in growing and processing AlxGai„xAs layers which are graded to high 
aluminium fractions over the short distances involved. 

We note finally that in good material, a grade may in fact reduce the overall 
Jsc- Figure 15.121 b) shows the extreme example of a series of Jgc predictions 
starting with a good GaAs cell. In this case the QE of the cell is rapidly degraded 
as the efficient GaAs material at the front of the cell is replaced with AlxGai_xAs 
with increasing aluminium fraction. 

5.2.5 Windows 

Figure 1^.131 shows the QE of an ungraded p layer with and without a window for 
a recombination velocity of 500m/s, which is frequently observed in our samples. 



5.2. Further QE Enhancements Above the Barrier Bandgap 



111 



^ 17 
^16.5 

U 16 

4.5 


a) 
























































_ O O 


o 


o 


o o ° 










o 

o 


o 












o 
















o 


















— " 1 4 
13.5 
13 






o 






















o 














o 




































12.5 


















































12 

E - 

O 22 


























0.3 

b) 


0.4 


0.5 


0.6 


0.7 


0.8 

Xbi % 


o 





































o 




















<■ 20 

^ 18 






m. . 






















o 






































S 16 










....O 








































^ 14 
12 






























o 
















o 


o 






...Q O ' 












o 


o 


o 


O 








8 












1 


0.2 0.3 


4 


0.5 





6 0.7 

Xbi ^ 



Figure 5.12: Jsc oi a. graded p - i - n cell as a function of the higher bandgap X^i 
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If we assume the recombination velocity is unchanged, we see a significant 
increase in QE because less light is lost before the incident spectrum reaches the 
p layer. However, removal of the window layer is expected to increase surface 
recombination. The lowest curve on the graph shows the worst case scenario of a 
recombination which is large enough to force the minority carrier concentration to 
negligible levels at the front of the cell. We note that increased light transmission 
in this case cannot compensate the loss of carriers at the surface. 

The behaviour in graded structures however is different. Appendix 1X1 shows 
that the effect of surface recombination is reduced by the effective field. In good 
material with appropriate contacts (see chapter [H]), the surface recombination 
velocity is expected to be very small compared to the effective field term. 

Figure 15.141 shows calculations analogous to those of figure 15.131 but for a 
structure graded from 30% to 67% aluminium. In this case, both calculations 
assuming an unchanged surface recombination and a worst case infinite surface 
recombination after removal of the window predict an enhanced QE. 

Although the predicted gain in QE is not large, we conclude that window 
removal may be beneficial because the graded p layer fulfills the function which 
the window is designed to perform. 

5.3 Further QE Enhancement Below the Bar- 
rier Bandgap 

The quantum wells in a standard QWSC design typically absorb less than half 
the incident light. The remaining fraction is absorbed in the substrate and does 
not contribute to the photocurrent. 

The quantum well photocurrent can be increased by ensuring that the light is 
reflected back through the cell after the first pass. This is achieved by removing 
the GaAs substrate and coating the back of the n layer with a metallic back 
surface mirror, as described in chapter IH] 

This technique has consequences for all optical functions associated with the 
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Figure 5.13: Influence of window layer removal on the QE of an ungraded p layer, including 
surface recombination effects 
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Figure 5.14: Influence of window layer removal on the QE of a graded p layer QE, showing 
weak dependence on surface recombination velocity 
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Fi gure 5.15; Light flux in a mirror backed QWSC 

cell, which are examined in the following sections. 

The modelling results presented assume a measured wavelength dependent 
front surface reflectivity and a back surface reflectivity of 95%. 

5.3.1 Generation, Flux and QE in Mirrored Cells 

Figure l5 . 1 51 shows how the addition of a back mirror alters the light intensity as a 
function of position in a QWSC. The light levels at long wavelengths are signifi- 
cantly increased. The main overall consequence is a more uniform light intensity 
throughout the solar cell. The position and wavelength dependent Fabry-Perot 
oscillations in light intensity are due to interference phenomena. 

Figure ISTTBl shows the generation rate of a 30 well graded QWSC as a function 
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of position and wavelength. The extended quantum well generation is clearly 
visible for wavelengths between 700nm and 870nm. The detail of the generation 
in individual wells is not visible since, as we saw in section El the model considers 
the i region as a weighted average of quantum well and bulk material. 

Figure 15.171 gives the generation rate for an identical cell with a back surface 
mirror, again as a function of wavelength and position. Interference effects impose 
a strong position dependent oscillation in the quantum well generation which is 
absent in the non mirrored device. 

The magnitude of the generation enhancement is visible in figure [5?TSk . which 
shows the generation rates in both cases, integrated over wavelengths from 400nm 
to 900nm. Significant enhnancements are seen for the i and n regions. In the 
design used for this example, insufficient light makes a second pass through the 
light to significantly affect the generation rate in the p layer. 

Figure 15.18b shows the enhancment in the i layer, expressed as a percentage 
of the generation rate in the cell without a mirror. Since the QE of the i layer 
is proportional to the i layer generation rate, an enhancement of approximately 
30% in QE is expected for this example. 

Finally, figure 15.191 shows the QE enhancement expected in a 30 well un- 
graded QWSC with a back mirror, assuming 95% specular reflection from the 
back mirror and an internal front surface reflectivity given by the mean 70nm 
SiN AR coat. The breakdown in QE enhancement layer by layer is shown in 
figure 15.201 Contributions to the QE increase from the neutral p and n regions 
are expected to be negligible. 

Enhancement for photon energies above Ef,2 is minimal since most of the light 
at these energies is absorbed during the first light pass. 

The wavelength at which the Fabry-Perot oscillations in the QE spectrum 
appear are critically determined by the bulk bandgap Eh2. This provides a very 
useful characterisation tool for determnining the exact value of Eb2 in QWSC 
samples 

The overall predicted i region enhancement in the quantum well absorption 
range is of the order of 25% for such a device. This enhancement ratio decreases 
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Figure 5.16: Generation rate in a 30 quantum well QWSC 
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Figure 5.17: Generation rate in a mirror backed QWSC 
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Figure 5.18: Generation rate enhancement a) in all layers of a graded QWSC with and 
without a back mirror and b) percentage increase in the intrinsic region. 
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Figure 5.19: Quantum efficiency of a thirty well QWSC with a back mirror compared to an 
indentical cell without a mirror. 
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Figure 5.20: Absolute QE increase in a thirty well QWSC with a back mirror compared to 
an indentical cell without a mirror. 
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to about 20% for fifty well samples because of the greater light absorption on the 
first pass. 

5.4 Conclusions 

In the preceding sections we have reviewed previous work on QE and Jgc en- 
hancement in ungraded cells and explained the limitations of the optimisation in 
ungraded QWSCs. 

We have separately considered further enhancements above the barrier bandgap 
for graded QWSCs and below the barrier bandgap for any QWSC cell. 

Above the barrier bandgap we have seen that the spatial variation of optical 
and transport parameters increases the QE by reducing minority carrier genera- 
tion near the surface of the solar cell and by increasing light transmission to the 
highly efficient i region. 

The compositional grade results in improved carrier collection in a graded p 
layer because of the presence of an effective field. The consequent improvement 
in QE due to this mechanism is difficult to predict quantitatively because of the 
added complexity introduced by the spatial variation of materials parameters, 
and the assumptions made necessary by the lack of knowledge regarding the 
position dependence of transport parameters in particular. We have presented a 
modelling scheme which uses previously derived theoretical results in conjunction 
with published and measured experimental data to permit modelling of the QE 
of graded layers with a minimum set of free parameters. 

The model in this case can only be used as a qualitative guide to suggest 
optimal parameters for graded p layers, but may be improved by comparison 
with experimental data. 

Below the barrier bandgap the Jsc from the quantum wells can be greatly 
enhanced by coating the back of the QWSC with a mirror. The large fraction of 
light which is not absorbed upon a first pass is reflected back and forth in the 
cell, increasing light absorption in the quantum wells and hence the Jgc- 



Chapter 6 

Device fabrication and 
characterisation 

The QWSC design is based on a MQW photodiode detector described by White- 
head [36]. Each successive generation of QWSC design passes through design, 
growth, processing and characterisation stages. 

6.1 Sample Growth 

Successive growth runs are determined by a combination of the characterisation 
of earher samples and optimisations carried out using the short circuit current 
modeL Sample growth was carried out using two epitaxial growth techniques, 
which are briefly described below. 

6.1.1 Metalorganic Vapour Phase Epitaxy (MOVPE) 

MOVPE growth was carried out by John Roberts at the EPSRC III-V Semicon- 
ductor Facility at Sheffield University using a Quantax MOVPE reactor. Details 
are given in reference [70]. 

MOVPE growth consists in passing reagent chemicals in vapour form over a 
heated substrate. The precursor gases dissociate at the surface of the substrate, 
depositing atoms at the rate of about one monolayer per second. 
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Substrate temperatures in our case varied between 560C and 600C. For the 
AlxGai_xAs system, the reagents were trimethyl gaUium (TMG, Ga(CH3)3), 
trimethyl aluminium (TMA, A1(CH3)3) and arsine (ASH3). 

The p-type dopant in early studies was zinc in the form of dimethyl zinc 
(DMZ, Zn(CH3)2). This was later switched to carbon doping in the form of CCI4 
to reduce p-type dopant diffusion, n-type doping was achieved using disilane 
(Si2H6). These reagents are transported to the growth chamber in a hydrogen 
carrier gas at atmospheric pressure. 

The MOVPE growth technique has a number of attractive features. The most 
interesting from the perspective of solar cell production is the high productiv- 
ity, given that commercial MOVPE reactors are currently capable of growth on 
many wafers at the same time. Moreover, the gas sources used in growth can 
be replenished without shutting down the reactor. MOVPE growth is therefore 
significantly cheaper than some other epitaxial growth methods. 

Compositional control is, in principle, lower in MOVPE than in MBE and 
suffers from a lack of in-situ measurement techniques such as reflection high 
energy electron diffraction. Accordingly, growth mechanisms in MOVPE are less 
well understood than in MBE in particular. Nevertheless, near monolayer control 
is routinely achieved in good MOVPE material. 

However, most of the MOVPE material we have characterised suffers from 
electrical problems due to high background doping levels. As we saw in chapter 
[21 a high background doping adversely affects the open circuit voltage of the 
QWSC Early material showed tolerable net dopant densities. We shall see in 
subsequent sections that a lower silicon background doping density, for example, 
may in fact lead to a higher background doping density due to the loss of a 
compensation effect for the carbon background. 

6.1.2 Molecular Beam Epitaxy (MBE) 

Molecular beam epitaxy was performed by Christine Roberts in the III-V Inter- 
disciplinary Centre for Semiconductor Materials in Imperial College. 
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MBE growth consists in firing coUimated beams of molecules at a substrate 
lieated to between 480C and 630C, wliicli is situated in an ultra liigli vacuum 
growtfi cliamber. Collimated molecular beams are generated in fieated crucibles, 
wfiicli contain tlie precursors in solid or liquid form. Evaporation produces the 
molecular beam. The molecules dissociate at the surface of the substrate, de- 
positing the metal atoms. 

Compositional control is achieved by changing the temperature of the solid 
sources. Very sharp compositional changes of the order of a monolayer are pos- 
sible simply by placing a shutter in front of any given source. The sample is 
generally rotated during growth to ensure uniform growth across the wafer. 

Dopants in this study are Be for p-type and Si for n-type AlxGai_xAs. As we 
shall see in subsequent sections, the Be dopant is prone to significant diffusion at 
the doping levels and growth temperatures used in this study. 

This growth method is relatively well understood, since an array of char- 
acterisation tools may be placed in the growth chamber. A review of these 
characterisation tools has been written by Foxon |38] . 

The method is complementary to the MOVPE technique. It is inherently more 
costly, but can achieve higher purity due partly to the use of UHV techniques, 
cryogenic paneling, the wide range of in-situ characterisation techniques available, 
and the relatively well understood growth mechanisms. 

Although the background doping levels in MBE material were generally low, 
decreasing material quality with increasing aluminium fraction has consistently 
been observed. 

6.2 Processing 

Processing was carried out by Malcolm Pate at the III-V Semiconductor facility 
at Sheffield. This last stage before obtaining finished devices is critical, since the 
devices are exposed to potentially damaging environments. 

The most common device configuration is the mesa photodiode. The early 
studies described in reference ^ used both anti-refiection coated (AR) and non 
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AR coated (NAR) versions of this structure in order to save processing time on 
poor material. This study concentrates on the AR coated devices. 

The principal stages of the processing routine are tried and tested, and involve 
relatively few processing steps. Reference [61] gives a summary. The principle 
innovation is the addition of a back surface mirror. We now outline the main 
steps involved in turning a piece of QWSC wafer into the finished device. 

6.2.1 Device Geometry 

The photodiode devices are lOOOyum circular mesa structures defined by pho- 
tolithography followed by an etch. Gold ring contacts are then laid down on the 
mesa, defining the 600 fim diameter optical window. 

NAR devices retain the protective ~ 40nm thick GaAs capping layer which 
is grown on the front surface of the cell to prevent AlxGai_xAs oxidation and 
surface damage in handling. AR devices are subjected to two further steps. The 
first of these is the removal of the GaAs capping layer with a selective etch which 
stops short at the AlxGai_xAs window, which acts as an etch-stop. A reliable 
etch stop requires an aluminium fraction of at least 30% Al. 

Following the removal of the cap, a SiN AR coat of thickness ~ 70A is de- 
posited onto the optical window. This step involves a brief exposure of the un- 
protected device to air, followed by device exposure to a silane/nitrogen plasma 
deposition machine at about 300C. SiN is deposited at a rate of approximately 
50nm per minute. 

Where possible, a piece of unprocessed QWSC wafer is AR coated at the same 
time as the devices, since these are too small to allow accurate measurement of 
the reflectivity. Such pieces of wafer were received for about half the samples 
processed. 

Design and processing of large area devices follows similar steps and is de- 
scribed in [66] . 
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6.2.2 Mirror Backed Devices 

The addition of a back mirror to the QWSC devices was first tried with mesa 
photodiode devices. The initial steps are identical to those listed in the previous 
section. In addition to these, however, the light absorbing GaAs substrate is 
removed only under the optical window with a selective etch which uses the n 
layer as an etch-stop. 

This leaves a ~ 300/xm hole of diameter ~ 600//m. The gold metal consti- 
tuting the mirror is then deposited in this hole, in direct contact with the back 
surface of the n layer. 

6.2.3 Effects of Processing 

Exposure of QWSC wafers to processing leads to unavoidable consequences in 
device performance. Effects which are directly related to processing include the 
wavelength optimisation of the AR coat and the electrical characteristics of the 
solar cells. The former may be improved by tuning the thickness of the AR coat 
to suit a particular cell, whereas the electrical characteristics such as series and 
parallel resistance may be optimised by the choice of appropriate alloys and mask 
sets at the contacting stage. 

Other processing effects are accidental and relate to the internal quantum effi- 
ciency of the devices. The removal of the cap has been known to cause problems 
in some MBE material for reasons as yet unknown, leading to a poor window 
surface. Malcolm Pate has reported difficulties with the etch stop in some cases. 
This may be due to poor material quality at the surface, or an unexpectedly low 
aluminium mole fraction at the front surface. 

Another source of damage is oxidation following the removal of the protective 
cap and subsequent exposure of the topmost AlxGai_xAs layers of the cell to 
air. Since these AlxGai_xAs layers are particularly susceptible to oxidation, this 
exposure must be as brief as possible to avoid oxygen diffusion through the thin 
AlxGai_xAs window to the p layer beneath. 

A general decrease in device quality is generally seen as the number of pro- 
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cessing steps is increased. This is particularly visible in mirrored samples which 
show variable device quality. Surface contamination is frequently observed in 
these samples. 

On the other hand, we shall see in subsequent sections that there is evidence 
to suggest that the AR coating stage has a passivating effect, leading to improved 
electrical characteristics in forward bias for initially poor devices. High quality 
samples on the other hand are generally degraded by the AR coating process. 

We conclude that the processing stage unfortunately introduces a variability 
in the characteristics of the finished QWSC devices. This variability is difficult 
to ascribe to any given cause, even with the limited number of steps involved. 

6.3 Characterisation 

Experimental data were obtained using equipment designed and built by Mark 
Paxman and described in reference [Hj. The two routine characterisation tools 
are photocurrent (PC) spectroscopy, monochromatic current- voltage (MIV), and 
measurements of the anti- reflection coating (AR measurements). 

Other characterisation techniques carried out by collaborators include trans- 
mission electron microscopy (TEM) and secondary ion mass spectroscopy (SIMS). 

6.3.1 Photoconductivity Experiment 

QE, MIV and AR measurements were all carried out on the photoconductivity 
rig which has been extensively covered in references [S], |SZ] and [52] • 

Illumination is provided by a tungsten-halogen lamp. The light is passed 
through a grating monochromator which is capabable of wavelength resolution 
of the order of 0.3nm (reference [8]). Two power supply units (PSUs) enable the 
device to be driven to reverse or forward bias. 

A pinhole placed at the exit slit of the monochromator transmits light to a 
collimator. The image of the pinhole is then projected on to the surface of the 
device. The diameter of the pinhole must be made smaller than that of the optical 
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window in order to ensure that no light is lost by reflection from the contacts. 
The light power incident on the device is usually of the order of lOOnW. 

The light beam is modulated by a chopper wheel at ~ 180Hz. A lock- in 
amplifier is then used to measure the resulting photocurrent. The use of lock- in 
techniques ensures that the DC background due to ambient light and the dark 
current is removed from the measurement. 

The incident power is measured by a calibrated silicon photodiode placed at 
the same position as the QWSC This additional calibration measurement can 
then be used to compute the QE as a function of wavelength. 

The experiment is controlled by a computer running a Turbo Pascal program 
written by Mark Paxman. It should be noted that measurements at short wave- 
lengths are prone to large errors because of the low reponse of the Si photodiode 
in this wavelength range. The error over most of the wavelength range has been 
quoted as ~ ±5% in reference [8]. 

Measurements performed on this rig are AC and are designed to eliminate DC 
signals due to the dark currents of both the solar cell devices and the detector. 
The resulting measurement gives the photocurrent, as opposed to the DC light 
current which includes a dark current component. 

6.3.2 Monochromatic current- volt age characteristics 

Monochromatic IV (MIV) measurements consist of measuring photocurrent as a 
function of applied bias at a fixed wavelength. The measurement is normalised 
to 100% at the voltage independent saturation level reached in reverse bias. This 
technique provides information on the quality of both growth and processing. 

Measurements are normally carried out above the bulk band-edge in order 
to prevent distortion of the MIV with voltage because of the modification of 
quantum well absorption with Stark shift. Below the bulk band-edge, full QE 
measurements (see section 16.3.31 below) as a function of bias must be taken in 
order to identify a wavelength range in which the QE is relatively insensitive to 
the Stark shift. 
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Figure 6.1: Typical MIV curves for good and moderate solar cell devices. 



Typical MIV curves for good and moderate devices are shown in figure 16.11 
Tlie "square" curve, sfiowing a MIV which starts to fall off at zero bias was 
measured for a 30 well MBE grown ungraded QWSC with Xb2 — 20%. The 
second curve showing a MIV which starts to decrease at -4V represents data for 
an identical sample grown in the same series with Xf,2 — 30%. 

Good devices such as the 20% sample maintain a constant photocurrent into 
forward bias. At a moderate forward bias, the signal falls abruptly. A cut-off 
occurring in forward bias may be due to the failure of the experimental technique, 
since the lock-in amplifier begins to detect the increasingly large AC variation in 
the dark current caused by the chopped light. Details of these problems have been 
given by Paxman [S] and Barnes [57]. In reverse bias, the MIV remains constant 
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until the breakdown voltage, not shown here because the MIV measurement 
does not measure the dark current. At this point, the dark current increases very 
rapidly. Detection of the photocurrent against this background becomes difficult, 
and the resulting data noisy. 

Devices with significant background doping levels such as the 30% device show 
a photocurrent which falls off at voltages well below the Voc- This current drop 
is due not to the increasing influence of the dark current, but to the decrease in 
field across the intrinsic region, and subsequent failure of the p - i - n structure 
to operate as a solar cell. 

With increasing reverse bias, the PC of such cells reaches a maximum. We 
saw in section 12.3.21 that the voltage at which the drop in photocurrent occurs 
can be used to estimate the background doping. The method however becomes 
inaccurate for devices with doping levels in excess of ~ lO^^cm"^ and in samples 
which do not conform to the abrupt p - i - n doping proflle assumed in the model. 

Device quality 

Many devices of lower quality depart signiflcantly from the picture described. 
This may be due either to growth or processing problems. 

Samples with doping proflles which do not conform to the sharp p - i - n 
doping proflle assumed in the previous discussion will behave in a more complex 
manner. If for example the voltage is varied from negative to positive bias, the 
p and n side depletion widths will vary in a manner which reflects the doping 
proflles on both sides. An estimate of this proflle may only be obtained if the 
one is much greater than the other. 

In practice, this method cannot be exploited further because such variations 
may also result from device history. Difficulties in processing in particular may 
lead to electrical problems. This may be manifest as a very noisy MIV mea- 
surement or as the absence of a plateau in reverse bias allowing the curve to be 
normalised to 100%. 

We conclude that the MIV measurement is useful in identifying promising 
devices and furthermore indicates whether the devices conform to the idealised 



6.3. Characterisation 



132 



structure assumed by the modelling programs. 

6.3.3 Quantum efficiency measurements 

Measurements quoted in chapter [7] are always taken at voltages for which the i 
region of the cell is fully depleted, in order to maintain the field across the i region 
and fulfill the assumptions in the model. Although the PC spectrum in the wells 
moves with applied bias due to the Stark shift, this dependence is quoted by 
reference [50] as being of the order of meV for typical biasses of up to a few volts 
and shifts the quantum well absorption edge by only a few nm for the moderate 
biasses used in the measurements. 

The MIV technique of the previous subsection provides a useful guide to the 
voltage at which the QE measurement should be taken. 

6.3.4 Reflectivity measurement 

The last routine characterisation tool is the reflectivity experiment. A slight 
modification of the experimental setup described above is required. As mentioned 
in section 16.21 pieces of wafer are generally AR coated at the same time as the 
photodiode devices. Direct reflectivity measurements are possible on large area 
devices. 

The reflectivity measurement consists of placing the AR coated sample at 
the position normally occupied by a device. Light reflected from the surface 
is collected and measured by the Si photodiode at as low an angle as possible. 
The sample is then removed, and the incident power measured directly by the Si 
photodiode, the optical path length being kept constant. The reflectivity is then 
given by the ratio of the two signals. Beam angles are typically below 5 degrees. 
Measurements were repeated on different dates and a maximum error of about 
5% in reflectivity at long wavelengths was found. 

A good knowledge of the reflectivity is required for accurate modelling of QE 
spectra. For some photodiode devices however, pieces of AR coated wafer were 
unavailable. In this case, we are forced to use an average of previous measure- 
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ments for the AR coating method used. A significant error is introduced in this 
case, which may be quahtatively reduced by examining the devices under a hght. 
Typical AR coats are a characteristic blue colour. 

Reflectivity measurements give additional indications of processing or growth 
problems. This is typified by a large variation in reflectivity over the AR coated 
sample. 

6.3.5 External characterisation 

In addition to the experiments described above, more sophisticated characterisa- 
tion techniques were applied to a selection of structures by collaborators. 

Transmission electron miscroscopy (TEM) 

TEM measurements were carried out by Professor Grunbaum in the Materials 
Department of Oxford University. 

Transmission electron microscopy consists of flring a collimated electron beam 
through a thinned cross section of a sample. The electron beam is diffracted by 
the crystal lattice, and the diffraction pattern recorded on a fluorescent screen or 
on photographic fllm. Details of this technique are given in the review articles 
by Holt ^5 and Hutchinson [61]. 

Resolution of the order of O.lSnm is possible on state of the art microscopes. 
Compositional contrast is possible because of the different electron beam scat- 
tering factors for AlxGai_xAs with different values of x. TEM micrographs can 
therefore give unprecedented direct measurement of the physical dimensions of 
QWSC samples. Measurements were carried out on selected samples to verify 
the dimensions of the quantum wells in particular, and of the larger scale p layer 
thickness. 

Secondary ion mass spectroscopy (SIMS) 

SIMS measurements were carried out by Federico Caccavale of Padua University. 
The technique is described in detail by Clegg [65]. This measurement involves 
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firing a collimated beam of ions at tlie sampie. Tfie liigfi energy incident beam 
removes atoms from tlie first few monolayers of tlie sample. The flux of these 
sputtered atoms is then measured by a mass spectrometer. 

The incident beam progressively digs a pit in the surface of the sample, pro- 
viding compositional information as a function of depth. The technique must be 
tuned by choice of incident atomic beam and mass spectrometer setting to detect 
the spatial variation of different elements. The focus in the two measurements 
carried out for this project was the variation of doping density and aluminium 
composition with depth. 

Depth resolution is dependent on factors including the type of atom used in 
the incident beam and the beam energy. Use of oblique incident atomic beams 
and non-normal incidence to prevent mixing at the surface of the sample can 
achieve resolutions of the order of a few nm. 

6.4 Conclusion 

In this chapter we have described growth, fabrication and characterisation of 
QWSC devices. We have seen that the behaviour of the finished devices may 
depart significantly from the theoretical picture given in previous chapters. Both 
growth and processing may contribute to this deviation. 

The most readily identifiable problem is that of background doping. This is 
revealed by MIV measurements, which are also used to suggest an appropriate 
bias at which to measure the QE. 

Although QE measurements may be taken at zero bias in good material, it 
is frequently necessary to force the device into reverse bias to ensure that the 
approximations in the model arc consistent with the experiment. 

Reflectivity measurements are essential to a correct analysis of the QE of a 
device. In some cases such measurements had to be extrapolated from an average 
over previous AR coating runs. 

Apart from these routine experimental tools, we have reviewed the SIMS and 
TEM characterisation techniques which were carried out by collaborators. These 
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additional characterisation methods may be used to verify or adjust the geometry 
and composition of QWSC samples. 



Chapter 7 



Preliminary Characterisation 

In this chapter we present results of the MIV and reflectivity techniques reviewed 
in chapter [6l These are the preliminary characterisation tools used to make 
an estimate of the device quality. The reflectivity furthermore is required for 
subsequent analysis of quantum efficiency data. 

The samples in this chapter are referred to by their growth name, and relevant 
design characteristics are given in the text. 

7.1 Monochromatic IV Measurements 

A good device is characterised by a MIV curve which remains constant into 
forward bias, before showing an abrupt cut-off for reasons given in chapter |6l 
If this cut-off occurs at forward biasses greater than about 0.6V, attribution is 
difficult because of experimental considerations. 

If the cut-off occurs at reverse bias and the acceptor, net background and 
donor doping levels A^^, Nbg and No layers are abrupt, the drop in photocurrent 
can be related to the net background doping level A''^^, as described in chapter 

m 

Calculations of Nbg for a range of cut-off voltages are shown in table 17.11 
The variation of predicted Nbg is only weakly dependent on the cut-off voltage 
Vs- The technique can therefore only provide a rough estimate of Nbg-, because 
the values of Vs generally vary only within a few volts. A further limitation is 
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Vs (V) 





-2 


-4 


-6 


-8 


-10 


Nbg (lO^^cm-^) 


1.0 


2.0 


3.0 


4.0 


5.8 


5.9 



Table 7.1: Predicted value of net background doping Nbg for a range of MIV cut-off voltages 
Vs, for a standard 30 well QWSC sample with Xb2 = 30% 

introduced by the poorly understood electrical influences on the MIV character- 
istic, and by the fact that the abrupt doping assumption is frequently erroneous, 
as shown by SIMS measurements in chapter |8l The position dependence of the 
p layer doping in particular has been investigated by SIMS in both MOVPE and 
MBE samples. This position dependence introduces a further voltage depen- 
dence which makes the precise attribution of the MIV voltage shoulder difficult. 
Furthermore, the shoulder is generally imperfectly defined in poor samples, and 
may be subject to errors of several volts. 

Jenny Nelson has developed a model, in collaboration with Alex Zachariou, 
which estimates the background doping by modelling the QE at zero volts, and 
using Nbg as the main fitting parameter. Varying Nbg in this model alters the 
effective p, i and n layer thicknesses. The p layer in this case is modelled as two 
undepleted layers with doping levels Na and Nbg respectively. This approach 
allows a more quantitative estimation of Nbg, but nevertheless still assumes 
square doping profiles. 

The model presented in this thesis does not use this approach. The MIV data 
in the following discussion is principally used as an indication of material quality. 
In particular, we will use it to qualitatively assess those samples in which Nbg 
is large enough to degrade device operation. We further use the MIV data as 
indicative of the reverse bias at which to take QE measurements. 

MIV data for QWSC devices grown for this project is generally poorer than 
the corresponding data for the p - i - n controls. We now examine this trend as a 
function of growth method and describe the effect of changing growth conditions 
on this problem. 

In the final part of this section we look at the effect of the AR coating process 
on the MIV characteristics of the samples, and compare these results with data 
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Figure 7.1: MIV data for an ungraded MOVPE grown 30 well QWSC and p - i - n control 
pair with Xi,2 — 33% aluminium before AR coating, showing deterioration of MIV in QWSC 
devices 

from a passivation study carried out at Eindhoven University of Technology as a 
result of these observations. 



7.1.1 Trends in QWSC and p - i - n Device MIV Data 
MOVPE Devices 

Figure O shows MIV data for the MOVPE grown 30 well QWSC and p - i 
- n control pair QT468a and QT468b. The nominal aluminium composition is 
Xb2 = 33%. Neither device in this example is AR coated in order to minimize 
the variability introduced by processing. Also shown is data for a 50 well QWSC 
grown previously on a different MOVPE reactor. 

The MIV data for the 30 well QWSC device shows a voltage independent 
plateau in reverse bias. A cut-off is seen at about -2V. Assuming that this 
current decrease is due to a loss of field in the i region as explained in chapter 
El we can estimate the net background doping as Nbg — 2 x lO^^cm"^. The net 
background doping for our MOVPE material is p type residual carbon [56] . 
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The curve for the p - i - n device is constant up to the point at which the 
measurement becomes unrehable for reasons given in chapter O The MIV cut-off 
in this case cannot therefore rehably be attributed to a loss of field in the i region. 
The cut-off point at about 0.8V however allows us to estimate the upper limit of 
the net background doping as Nbg < 5 x lO^^cm"^. 

The 50 well device MV370 of figure 17.11 was grown on the Sheffield MRIOO 
machine, which was subsequently removed from the AlxGai_xAs growth program. 
The high quality MIV of this 50 well sample shows that QWSC devices with a 
large number of wells and a low net background doping can be grown by MOVPE. 
Data for other high quality MOVPE grown QWSC samples of this type have been 
published |19]. The data in this case again only allows us to set an upper limit 
estimate on the net background doping, which is Nbg < 2 x lO^^cm"^. 

John Roberts [5S] has suggested that the systematic difference in MIV seen 
between p - i - n and MQW samples may be caused by a lack of self compensation 
between oxygen and carbon impurities in samples with multiple quantum wells 
due to reduced oxygen impurity levels in the quantum wells. This leads to a 
higher than expected p type carbon background. Moreover, the difference in 
net background doping levels on either side of a quantum well barrier interface 
is expected to lead to poor interface quality and high interface recombination 
velocities. 

John Roberts has further suggested that AlxGai_xAs material will have poor 
minority carrier transport properties despite the low net background doping lev- 
els, because of the high impurity levels. We shall see in subsequent sections that 
the AlxGai_xAs minority carrier transport properties decrease with increasing 
aluminium content, even though the net background doping levels observed in p 
- i - n structures such as that of figure 17.11 appear low. 

Very poor MIV characteristics have been recorded for some MOVPE samples. 
Figure 17.21 shows the MIV data for optimised sample QT794a. This 50 well 
sample has no window, and a thin p layer nominally graded from X^i = 40% to 
Xh2 = 30%. The MIV data in this case begins to drop off at -9.5V, suggesting 
Nbg — 3 X lO^^cm"'^. SIMS data taken by Federico Caccavale were aimed at 
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Figure 7.2: MIV measurements for optimised device QT794 for NAR and AR coated devices 

verifing sample geometry, composition and intentional doping profiles. The p type 
carbon background doping sensitivity was ~ lO^^cm^^ and was unfortunately too 
insensitive to verify this value of Nbg- 

Further discussions concerning this sample in the QE modelling section of 
this chapter show that the minority carrier efficiency in this sample is severely 
reduced, as suggested by John Roberts. 

MBE Devices 

QWSC and p - i - n samples grown by MBE show a similar trend to MOVPE 
samples discussed in the previous section. Figure 17.31 shows the MIV data for 
three ungraded 30 well QWSC samples grown and processed in the same run. 
The nominal composition of these samples was Xb2 = 20% for U2027, X52 = 30% 
for U2029 and Xh2 = 40% for U2031. Similar data for their respective controls 
U2028, U2030 and U2032 which have matching composition but no wells is shown 
in figure 17.41 



The data for the 20% p - i - n U2028 has been smoothed to eliminate noise 
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Figure 7.3: MIV data for a series of 30 well QWSC ungraded devices grown by MBE with 
nominal aluminium fractions Xs2=20%, 30% and 40% 
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Fi gure 7.4: MIV data for a scries of p - i - n ungraded devices grown by MBE with nominal 
aluminium fractions X_b2=20%, 30% and 40% 
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Figure 7.5: MIV data for the 50 well graded sample U6041 and its p - i - n control U6045, 
indicating low net background doping 

spikes. We note that the MIV data at 20% are nearly ideal in both QWSC and 
p - i - n cases with estimated limits on the net background doping of Nbg < 
7 X lO^^cm^^. The same limit is found for the 30% and 40% aluminium p - i - n 
samples. 

The QWSC MIV data however become increasingly poor as the aluminium 
content is increased, whereas the p - i - n control MIV data remain good. Es- 
timates of the net background doping yield values of Nbg — 2 x lO^^cm"'^ for 
the 30% and 40% QWSC structures. Although the two curves are significantly 
different, the MIV cut-off method of determining Nbg is not sufficiently sensitive 
to differentiate between the net background doping levels in the two samples. A 
model such as Jenny Nelson's model which was mentioned in the introduction to 
this section would however be more sensitive to Nbg in this case. 

Finally, figure 17.51 shows MIV data for NAR devices processed from p - i - n 
and 50 well QWSC structures grown in the latter part of this project. Although 
the nominal lower aluminium fraction Xi,2 for these samples was 30%, modelling 
in both cases indicates a much lower Xi,2 — 24%. Similar behaviour is again 
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apparent, with the p - i - n sample performing better than the QWSC 

In this case, however, a phenomenological difference with the 30 well samples 
of figure 17.31 can be seen. Instead of a steady decrease from the cut-off voltage 
seen in reverse bias for the 30 well samples, the current for the 50 well device 
of figure 17.51 decreases slowly from about -3V, but is still at about 95% of the 
reverse bias maximum at zero bias. Furthermore, SIMS data presented in chapter 
[8] indicates diffusion of the Be p type dopant from the p layer into the i layer. 
This leads to an exponentially decreasing Be doping profile in the i layer. This 
means that the depletion width in the p layer will begin to decrease more rapidly 
with an applied bias varying in the positive direction, as the depletion edge scans 
through regions of the p layer with exponentially decreasing doping. 

We therefore ascribe the slow decrease in MIV from 100% at -3V to about 
95% at OV to the inhomogeneity in the doping profile, and the fact that it does 
not match our assumptions regarding square doping profiles. 

The voltage cut-off at about 0.6V is then too low to allow us to reliably 
suggest an absolute value of Nbg- However, we can again set an upper limit 
which is Nbg < 2 x 10^^. This agrees with the upper limit set by the SIMS 
measurements. Net background doping levels of this magnitude are sufficiently 
low to allow efficient solar cell operation for 50 well QWSC designs. 

Previous studies presented in Mark Paxman's thesis [S] include measurements 
of MBE material grown by Tom Foxon at Phillips Research Laboratories in 1990. 
This is sample G951, which was a 50 well QWSC with an aluminium fraction 
of Xh2 = 33%, and was shown to have a net background doping level of Nbg — 
3 X lO^^cm-^. 

As in the case of MOVPE material, we conclude that the background doping 
problem is not unsurmountable. 
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Figure 7.6: MIV measurements for the identical pair of QWSC devices QT468a and QT528. 
The latter was grown with a cleaner Ga source and shows an improved MIV characteristic. 

7.1.2 Growth Optimisation Studies 
MOVPE Ga Precursor Optimisation 

An attempt was made to reduce the background impurity concentration in MOVPE 
samples in collaboration with John Roberts of the III-V Semiconductor Growth 
Facility in Sheffield. This involved comparing device quality before and after an 
improved TMG source was used for the MOVPE growth. 

Figure 17.61 shows MIV measurements for QT468a and the repeat structure 
QT528. The latter is identical to QT468a but was grown with the cleaner TMG 
source. The MIV data for both NAR and AR devices is better for the material 
grown with the cleaner TMG source. This supports the principle that poor MIV 
data is indicative of high impurity levels in the growth precursors. 

We note however that the voltage at which the MIV drops is very similar 
in both cases. The simple theory outlined in chapter |2] would therefore predict 
the same net background doping level in both cases of Nbg — 2 x lO^^cm"^. 
This further indicates that the main source of background impurities is the TMA 
source. 
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Figure 7.7: MIV data for identical MOVPE single quantum well samples grown at 560C, 
580C and 600C, showing little variation. 

MOVPE Growth Temperature Optimisation 

MOVPE growth temperature may significantly affect growth morphology and 
impurity incorporation during growth. In particular, carbon incorporation de- 
creases at lower temperatures. This is offset by an increased defect density [56] . 

A series of identical single quantum well (SQW) samples were grown at three 
different growth temperatures in order to identify an optimum growth tempera- 
ture. The temperatures were 560C for QT574b, 580C for QT874a and 600C for 
QT874c. The MIV characteristics for these samples are given in figure 17.71 We 
again consider NAR devices in order to minimize the device variability introduced 
by the AR coating process. The p layer in all samples was grown at 700C. 

Although the data for 600C appears higher than the other two at OV, the 
difference is only of the order of 3%, which is comparable or smaller than the 
experimental uncertainty. The MIV measurement in this case cannot distinguish 
between the three samples. Clearer conclusions can be drawn from PL spectra 
taken by J.S. Roberts piB]. The relative strengths of PL signals measured from 
the wafers relative to the 560C sample were 3.5 for 580C and 1.35 for 600C. A 
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temperature of 580C was therefore proposed as the optimum by John Roberts. 
He suggests that this temperature represents a good compromise between in- 
creased background carbon incorporation at higher growth temperatures and 
greater defect densities at lower temperatures. 

7.1.3 Passivation effects 

Comparison of MIV Data of AR Coated and As-Grown Devices 

Figure ES shows MIV data for the p - i - n sample QT468b for samples with and 
without an AR coat. The initially good MIV characteristic is degraded by the AR 
coat deposition. The decrease in MIV in reverse bias however shows a constant 
gradient. Linear behaviour of this type may be due to a change in the effective 
doping profile brought about by the AR coating process. Other influences such as 
electrical problems related to the extra processing step cannot however be ruled 
out. 

Figure [73] shows similar data for QT640a, which is a QWSC sample with the 
same dimensions as QT468b but which has a graded p layer. The initially poor 
MIV characteristic is improved by the AR coating process. The improvement 
however does not significantly change the voltage point at which the MIV starts 
to fall. The net background for both NAR and AR devices is therefore estimated 
at Nbg ~ 2 X 10^^cm~3. 

The degradation of good samples and improvement of poor ones upon the 
deposition of the AR coat has frequently been observed in the devices processed 
for this project. 

Such behaviour was also observed in samples grown at Eindhoven Technical 
University. Peter Ragay et al. [S3] of the III-V Solar Cell group at Eindhoven 
have carried out passivation studies which show similar behaviour in QWSC 
samples. 

MIV data for the 50 well QWSC sample W276 is given in figure rTTUl The 
wafer was grown by MBE at Eindhoven Technical University. The two curves in 
the graph show the device characteristic before and after exposure to a hydrogen 
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Figure 7.8: Degradation of MIV for a good p - i - n device QT468b 
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Figure 7.9: Improvement in MIV for poor QWSC device QT640a 
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passivating plasma. Hydrogen ions diffuse through the crystal lattice and bond 
to impurities and defects. These may include materials defects, minority carrier 
traps, and ionised donor and acceptor impurities in the crystal. The mechanisms 
are reviewed for GaAs in [51] and for AlxGai_xAs in [52j. 

Although the MIV data for the passivated device is much improved, electrical 
characterisation by Guido Haarpainter indicated that the treatment significantly 
increased the series resistance of the device. This is ascribed to the degradation 
of the intentional p and n layer doping due to unintentional passivation of the 
donor and acceptor dopant atoms. 

The analogies between passivated and AR devices suggest that the hydrogen 
present in the silane plasma may have a passivating effect on AR devices. The 
exposure time of one to two minutes is however very much smaller than the typical 
times of 30 minutes quoted by Ragay [53] , and does not benefit from the driving 
field which is usually employed to assist ionic diffusion. This however must be 
set against the fact that the samples we are considering are very thin. The i layer 
in the samples quoted above is approximately 0.15/im from the surface of the 
sample, which is smaller than the thicknesses of several microns usually studied 
in the literature. Furthermore, our samples show a weaker change in MIV than 
the Eindhoven samples. Confirmation of this hypothesis would require further 
characterisation. 

Passivation effects in indium phosphide solar cells have been reported by 
Chatterjee and Ringel [51]. This study observed an increased Voc in passivated 
p"*" - n samples. Deep level transient spectroscopy confirmed that deep levels 
were passivated. Consequences for the series resistance and overall efficiency of 
the cells are not given in this reference. 

We conclude that passivation may be of use in improving the dark currents 
of solar cells by passivating deep levels, and may also decrease the Nbg- The 
expected degradation of the intentional doping however presents a significant 
problem. It is expected to increase the series resistance and hence reduce the 
efficiency of the solar cell. 



7.1. Monochromatic IV Measurements 



149 





















- 




— 
































- 


















- 










\ 

\ 




















\ 










AR 








\ 

\ 

\ 










AR + 


pass 


ivati 


on 


\: 

! 
:l 

i. 

























-3 -2.5 -2 -1.5 -1 -0.5 0.5 1 i.5 

Voltage (V) 



Figure 7.10: MIV data for the 50 well QWSC sample W276 before and after exposure to a 
passivating hydrogen plasma 

7.1.4 Conclusions 

High net background doping levels have proved to be a significant problem in 
both MOVPE and MBE material with values of A^'^^ qualitatively estimated 
as greater than lO^^cm"^ observed in both systems. The systematic difference 
between QWSC and p - i - n devices has been ascribed to improved compensa- 
tion properties in p - i - n structures. The background doping problem is more 
systematic in MOVPE grown material. 

Efforts to reduce the background doping levels in MOVPE material have so 
far met with little success and the high levels of Nbg estimated for this growth 
method for QWSC devices continue to pose problems. 

Previous MOVPE and more recent MBE material has shown that both sys- 
tems are capable of net background doping densities of the order of lO^^cm^^. 

The AR coating process has been shown to significantly alter the MIV char- 
acteristics of solar cell devices. In general, the AR coat degrades the MIV data of 
good devices but improves that of poor devices. Samples which fall into the latter 
category show similarities with a passivated sample from Eindhoven University 
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Figure 7.11: Reflectivities for QWSC and p - i - n devices for the U2027-U2032 MBE series 

of Technology. This is hkely to be due to unintentional passivation of intentional 
and unintentional impurities in our samples by the hydrogen in the silane plasma 
used for the AR coating processing step. It is difficult to make use of this effect in 
practice since the passivation process unavoidably alters the intentional doping 
levels. This has been indirectly observed as an increased series resistance in the 
passivated Eindhoven device which degrades solar cell efficiency. 



7.1.5 Reflectivity Measurements 

We now briefly consider spread and reproducibility of reflectivity measurements. 
These are important in determining the internal QE of devices, and cannot always 
be measured directly. 

Figure 17.11b shows the spread in reflectivity for a range of samples processed 
in the same AR coating run. Of the six samples shown in the flgure, only sample 
U2029 deviates signiflcantly. 

These data form the basis of the average SiN reflectivity curve which was 
discussed in chapter [61 This mean reflectivity is shown flgure 17.11b . and is used 
to model the reflectivity of those samples for which a dedicated reflectivity mea- 
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Figure 7.12: Reflectivities of a variety of samples with and without windows, showing that 
the removal of the window does not noticeably affect the reflectivity 

surement sample is not provided. 

Also shown in figure ITTTk is a standard refiectivity for the Eindhoven MgF:ZnS 
double layer AR coat deposited at Eindhoven. We note that it is optimised at 
longer wavelengths, but is small over a wider wavelength range than the SiN 
single layer coat. 

The curves in figure 17.121 however show refiectivity data for samples grown 
towards the end of the project, just before the AR coating kit was upgraded. 
A large spread is seen for the samples AR coated at this time. This is consis- 
tent with discussions with Malcolm Pate of the III-V Semiconductor Facility at 
Sheffield University regarding the problems experienced with the original depo- 
sition equipment prior to the upgrade. 

It has been suggested [49j that the thin window may be used in conjunction 
with the SiN layer to produce a double layer AR coat. The optimum suggested 
window thickness was 30nm for the 70nm SiN AR coat. 

The data of 17. 121 includes two samples with and without windows, but with a 
nominally indentical SiN AR coat. The window thickness in this case was 20nm 
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instead of 30nm. Both samples were coated in the same processing run. We note 
that the data for the large area device processed from the U6041 wafer, which 
has no window, is very similar to the mean reflectivity of figure 17.111 

The reflectivity of the fifty well sample U7049 is again very similar in magni- 
tude and wavelength dependence, despite the fact that it has a window. The two 
other reflectivity curves shown in figure 17.121 are very high at short wavelength, 
indicating that the SiN coat is thicker than expected. The reflectivity for these 
samples is effectively optimised at longer wavelengths. 

We note furthermore than if we were to normalize the reflectivities of these 
four samples to the same wavelength, the windowless sample U7050 would have 
the lowest reflectivity. 

The window thickness in this example is 50% different from the optimum. 
The large differences observed in these data however lead us to conclude that the 
reflectivity of our samples is primarily determined by factors other than the pres- 
ence or absence of a window. The reflectivity overall is reasonably reproducible, 
with the exception of samples processed towards of the project, just before the 
AR coating kit was overhauled by Malcolm Pate. 

7.1.6 Conclusions 

In this chapter we have considered results and implications of the MIV and 
reflectivity characterisation techniques. The first section demonstrated the use of 
the MIV measurement to estimate sample quality in MOVPE and MBE samples. 
The technique has been shown to qualitatively demonstrate the magnitude of the 
background doping Nbg- Quantitative analysis however requires further analysis 
of the kind developed by Jenny Nelson and Alexander Zachariou [52] . 

We established that QWSC structures generally suffer from higher levels of 
Nbg than their p - i - n controls. This has been assigned to the loss of a com- 
pensation mechanism in QWSC structures by discussion with John Roberts. We 
further saw that previously grown QWSC samples have demonstrated that this 
problem is not insurmountable. Some recently grown MBE samples have con- 
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firmed this, although they have not reached the standard laid down by a sample 
(G951) grown at Phillips Research Laboratories in 1991. Overall, we conclude 
that MOVPE material is more variable than MBE material, and is subject to 
higher values of Nbg and a wider departure from nominal design parameters. 

Consequences of the AR coating process for device performance have been 
considered. An interesting obervation is that the MIV characteristics of poor 
NAR samples tend to improve in corresponding AR devices. This has shown 
interesting correlation with work by Peter Ragay |53] on hydrogen passivation. 

Reflectivity measurements were considered and shown to be reasonably con- 
sistent. Samples processed towards the end of the project however showed more 
variability. An average reflectivity was presented in order to enable modelling of 
samples for which the reflectivity was not directly measurable. 



Chapter 8 

Experimental and Modelled QE 
Spectra 

The first section of tliis cliapter gives an overview of tlie techniques used to model 
selected types of devices and outhnes the information which may be extracted 
from the modelling, together with limitations in each case. 

We then present experimental and modelled QE spectra for a set of samples 
with unusual i region compositional profiles. We then examine a set of ungraded 
samples grown by MOVPE and MBE, in order to estimate the dependence of 
the minority carrier transport on aluminium fraction. A set of graded devices 
are characterised, using the information derived from the previous sections. The 
final section presents QE spectra for QWSC devices equipped with back surface 
mirrors. 

The samples are referred to by their growth name, and relevant design char- 
acteristics are given in the text. Supplementary information is included in the 
graphical output of the modelling program. 

8.1 Principles of Data Fitting 

In this section we present experimental QE data and the corresponding modelling 
results for a range of ungraded samples grown by MBE. We start by reviewing 
theoretical results from chapter |3] relevant to modelling experimental data, and 
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establishing the modelhng techniques which will be used throughout this section. 

8.1.1 Sample geometry 

Sample descriptions are generally accurate, because of the epitaxial methods used 
in sample growth. Compositions however frequently need to be adjusted. This is 
usually possible only if a sharp cut-off in QE is present, for example in p - i - n 
devices. Aluminium fractions in QWSC samples can be estimated by fitting the 
shape of the QE above the base aluminium fraction and by looking at the 
position of the excitonic peaks in the quantum well QE. 

Adjustments are occasionally necessary to the window thickness, and are jus- 
tified by the presence or absence of a characteristic dip in the QE for photon 
energies greater than the window bandgap. 

An additional problem in MBE material is dopant diffusion. The sample 
U6041 which will be described in some detail below was characterised by SIMS. 
This revealed that the effective p thickness was approximately 40nm larger than 
expected, which for this sample represents an increase of nearly 50%. This casts 
some doubt over modelling other MBE samples with nominal layer thicknesses. 
In practice, modelling has shown that the extra thickness is not important in 
most MBE samples because they have p layers which are thicker than U6041. 
Furthermore, these samples do not show the symptoms of this problem. These 
are an overestimated i region QE at wavelengths of about 600nm, or an overes- 
timation of the quantum well QE resulting from unintentional doping of wells 
near the p layer. 

8.1.2 Samples with efficient minority carrier transport 

Chapter [3] shows that the QE is dominated by different parameters in different 
wavelength ranges. In particular, the QE at photon wavelengths corresponding 
to the band edge energy Eii2 is dominated by the i region contribution, and by the 
absorption and diffusion lengths in the p and n electrostatically neutral layers. 
Section 18.21 will show that the i region QE of even quite complex samples 
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may be modelled accurately, as long as the i region geometry and composition 
are close to specifications. We will further establish that the model is useful as 
a diagnostic tool. It can be used to adjust compositions, in the cases where the 
latter determine a sharp cut-off in QE at a certain wavelength. 

In the most general case, fitting of the model to the experimental QE is sub- 
ject to three main parameters, which are the recombination velocity parameter 
for the front of the p layer S = Sn/Dn, and the p and n diffusion lengths L„ 
and Lp. The recombination velocity at the back of the n layer is not a significant 
parameter. This is because the n layer is negligible at short wavelengths, and 
small overall. Furthermore, the n layer thickness of ~ 0.6/im is much greater 
than Lp. 

Appendix |X] shows that it is impossible to draw any clear analytical conclu- 
sions without making assumptions about the relative magnitudes of one or the 
other of these parameters. Experience with modelling however shows that differ- 
ent parameters dominate in different wavelength ranges. The diffusion lengths 
Ln and Lp are most significant near the band edge energy Eii2, whereas S is 
important at short wavelengths. Fitting the theory to the data in this case is 
an iterative process involving repeated adjustment of the three parameters until 
a best fit is found for all wavelength ranges. Samples falling into this category 
are generally those with good material quality and thin i layers, such as the low 
aluminium fraction 30 well samples. 

We will present results from samples grown at Eindhoven Technical University 
which also conform to this picture. In this case, however, the QE at short 
wavelengths is good, suggesting a very low value of S. The fitting then reduces to 
a two parameter fit, with the short wavelength QE being exclusively determined 
by Ln, and Lp being estimated subsequently by fitting the band edge QE. 

8.1.3 Samples with simplified modelling 

More exact modelling is possible in samples conforming to the simplified analyt- 
ical expressions for the QE derived in Appendix |Xl 
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Samples with small n region contributions in particular allow the QE near the 
band edge to be fitted mainly in terms of L„. A further coincidental simplification 
arises from the fact that the n region contribution is generally small in samples 
where L„ is less than the effective p layer thickness x^^. This is due to the fact 
that both p and n contributions to the QE are decreased in poor material. In 
this case, results from chapter [3] and Appendix [A] have shown that the QE for 
photon energies near Eg2 is very weakly dependent on S. The overall QE in this 
case can essentially be modelled independently at short and long wavelengths, 
in terms of S and L„ only. Samples falling into this category include those with 
base aluminium fractions Xi,2 greater than about 30% and samples with thick i 
regions. 

8.1.4 Samples with high surface recombination parame- 
ters 

At short wavelengths, the QE generally depends on both the recombination pa- 
rameter S and the electron minority diffusion length L„. As we saw above, the 
discussion need only consider the value of S in the p layer. As mentioned in 
chapter [5] we cannot separately estimate Sn and D„. 

A consequence of the analytical model of chapter [3] is that values of S com- 
parable to the term 1/L„ are required to affect the QE significantly at short 
wavelengths. Since diffusion lengths in the majority of our samples are only a 
fraction of a micron, it follows that surface recombination only becomes signifi- 
cant if iS > 10^m~^. Values quoted in chapter H] however set an upper limit on S 
of about lO^m"^. The larger values required to fit some samples indicate either 
unusually low mobilities together with high recombination velocities, or with a 
p layer width which is greater than the specifications. More fundamentally, this 
behaviour may be due to a failure of the recombination velocity to adequately 
represent minority carrier transport at the front of the cell. 

The recombination velocity represents the high defect concentration near the 
surface as a sheet of charge, which acts as a carrier sink. In reality, these defects 
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are distributed over a finite depth, and may be caused not only by crystal defects 
such as dangling bonds, but by surface oxidation which may extend some distance 
into the sample. If this damage is extended over a depth which is not negligible 
compared to the thickness of the p layer, the latter will be effectively shortened, 
and very high surface recombination velocities required to reproduce experimental 
QE data. 

The following discussion attempts to model samples with recombination ve- 
locities and diffusivities which are consistent with published data. For samples 
with thin p layers and low diffusion lengths, this in practice means values of S 
which are too low to have a significant effect. 

When such an approach is not possible, we shall assume that the surface 
recombination picture does not accurately describe the minority transport at 
the front of the cell, and substitute a recombination parameter which is large 
enough to saturate this mechanism. In practice, this means values such that 
S > 1/Ln, which are effectively infinite insofar as they are large enough to reduce 
the minority carrier concentration at the front of the cell to negligible levels. 

8.1.5 Graded samples 

A number of important differences appear when modelling graded samples. The 
first is that the p region QE for wavelengths corresponding to energy Eb2 is 
generally small, particularly for grades with a large compositional gradient. The 
QE near the band edge in graded samples is therefore principally determined 
by the i and n region contributions. Since we will not consider samples with 
compositionally graded n layers, this simplifies the fit near band edge Ei,2. 

Fitting the overall QE of graded p samples is necessarily limited by the as- 
sumptions we are forced to make regarding minority carrier transport in com- 
positionally graded material. In particular, the assumption that the estimated 
dependence of L„ on aluminium fraction and position can be extrapolated from 
the values used to model samples with different Xh2 is clearly liable to lead to 
large errors. Figure 18.11 shows two QE calculations. Curve A has no diffusion 
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Figure 8.1: QE for two linear parametrisations of Ln as a function of composition. For 
curve A), L„ is kept constant at O.OSfim. Curve B represents the worst case scenario of 
Ln{0) = 0.08/Ltm and Ln{xwp) = 0.01/im 

length gradient, and curve B has a higher than expected Ln gradient. The calcu- 
lation is for a value ofVx = —4. We note that even in this extreme case, the QE 
is relatively insensitive to the variation in L„. The QE in this case is therefore 
determined primarily by the difTusivity field Vx defined in section 15.2.21 

The theoretical considerations in chapter|3]and empirical tests using the model 
have shown that the recombination parameter S has little effect on the QE of 
the p layer at short wavelengths. Fitting the model to experimental QE data for 
graded samples is greatly simplified, since Vx is the determining parameter at 
short wavelengths. 

Although the diffusivity field Vx is in principle determined solely by the 
diffusivity Dn, it will in practice be influenced by other parameters. This is a 
consequence of reducing the modelling to a two parameter flt which is dominated 
hj Vx- Inaccuracies in the assumptions we have made regarding the diffusion 
length, for example, are subsumed into the modelled value of Vx- 

Modelling graded structures therefore follows similar guidelines to modelling 
the QE of ungraded structures, with the difference that the L„ is eliminated as 
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a parameter and the band edge response is determined by the i and n region 
contributions. The short wavelength QE will be modelled in terms of and 
to a lesser extent the recombination parameter S. This simplified approach is 
made necessary by the lack of detailed knowledge of the position dependence of 
the minority carrier transport parameters. 

8.1.6 Modelling of the quantum well QE 

The quantum well QE has been modelled in previous work by Jenny Nelson, 
which has been reviewed in |19] . For photon energies corresponding to electronic 
levels near the top of the well, the model for this wavelength range underestimates 
the QE to some extent because of superlattice states which are not included, and 
because of the change from a two dimensional to a three dimension density of 
states, which the model assumes is abrupt. 

With the exception of this correction, the model has been shown to reliably 
reproduce the quantum well QE, including excitonic features, and is sufficiently 
accurate for our purposes. The model is expected to be accurate for the first 
quantum well continuum, which is deep enough in energy to approximate the 
theoretical model accurately. 

The parameters for the quantum well QE are kept fixed, since they are well 
understood and have been reviewed extensively \49\ . 

8.1.7 Mirror backed samples 

Modelling of mirror backed samples introduces two new parameters. The first 
is the phase change ipb experienced by light internally refiected from the back 
surface mirror. This is a variable quantity which is poorly known [37] . 

The parameter ipi, can however be accurately determined from the position 
of the wavelength dependent Fabry-Perot oscillations which are observed in mir- 
rored samples. 

In principle, no other parameters are required, since the phase change ipf upon 
internal reflection at the front surface semiconductor-air interface is radians [35] , 
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the back surface reflectivity is well known from the reflectivity of the relevant 
metal, and the front reflectivity is generally measured. 

In practice, we shall see that some samples show clear signs of non specular 
reflections from the back surface, and show an enhancement in QE which is 
greater than that expected for light reflected normally from the back surface. 
This case is further characterised by lower Fabry- Perot amplitudes than expected. 
In some samples, they may be absent altogether. 

Finally, although the front surface reflectivity is measured, the model fre- 
quently predicts larger Fabry-Perot peaks than are in practice observed. The 
modelling in this case can be improved by the introduction of a front reflectivity 
factor J-". This reduces the modelled amount of light reflected back into the cell 
upon the third pass, and is a factor of the front surface reflectivity -R/ of section 
13.3.31 A value of J-" smaller than 1 reduces the amount of light reflected back 
into the cell from the front surface. This reduces the amplitude of the Fabry 
Perot peaks, but has little or no effect on the overall enhancement factor. The 
origin of this factor is not clear, but may be due to the fact that the light is not 
incident normally, or to an unintentional roughnening at the back surface mirror. 
Furthermore, we saw in section [3.3.3l that the internal reflections occurring at the 
interfaces between the quantum wells and the barriers are not included in the 
model. A brief calculation showed however that the effect is expected to be small. 
More discussion of this point is given in section 18.61 Finally, the uncertainty in 
reflectivity also plays a role, since measured reflectivities were not available for 
the mirrored samples. 

8.2 Graded i Region Samples 

A series of samples grown by MBE in Eindhoven University of Technology was 
studied. All structures have identical, ungraded window, p and n layers but have 
different i region specifications. These samples were initially grown to compare 
the efficiency of QWSC structures, p - i - n controls and p - i - n heterostructures 
with graded i layers, as reported by Ragay [53j. These samples are also useful 
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Figure 8.2: g^; of 32 well QWSC Eindhoven structure W405 

because they allow us to test the accuracy of the QE model in reproducing 
experimental data for devices with complex i region designs. 

MIV data for these samples were uniform and reproducible, indicating good 
device quality and a low value of Nbg < 5 x lO^^cm"^. 

The samples consist of a QWSC device, its p - i - n control, a GaAs cell and 
two samples with linearly graded sections in the i layer. Although the nominal 
p layer thickness was 0.5/im, Peter Ragay has confirmed that the real thickness 
is 0.2yum for all devices. Furthermore, the nominal window layer of 0.04/im 
was estimated as at most 0.02/im by the modelling, on the basis of the high QE 
value at short wavengths. A thicker window would produce a characteristic sharp 
decrease in QE for energies above the window bandgap. The QE measurements 
for these Eindhoven samples were taken by Ernest Tsui. 

Figures \K2\ and [831 show the experimental data and best fit from the model for 
a QWSC structure W405 and its AlxGai_xAs control p - i - n W406. The i layer 
of W405 contains 32 wells of width 90A separated by 60A barriers. The i layer 
of the control W406 consists of 0.5/xm of AlxGai_xAs with a nominal aluminium 
fraction of Xb2 = 20%. The model however indicates that the true aluminium 
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Figure 8.3: QE of p - i - n Eindhoven sample W406 

content is about Xb2 ~ 24% because of the sharp cut-off in QE observed at about 
730nm. 

We first consider the QE for photon energies above the neutral layer bandgap 
Ef,2. The interface recombination parameter S is estimated as too small to have 
any effect on the QE at short wavelengths, and is set to zero. 

In this case of neghgible recombination parameter 5, the QE of the neutral 
layers depends only on L„, Lp and the reflectivity. The modelling assumes a 
standard reflectivity supplied by Peter Ragay, and assumes that the layer specifl- 
cations are correct, with the exception of the window layer thickness as mentioned 
above. 

Modelled diffusion lengths are of the order or greater than the neutral layer 
widths at about 0.3/xm. The hole diffusion length is slightly smaller than the 
electron diffusion length L„ as expected. These are close to Hamaker's values, 
which are L„ ~ 0.39/xm and Lp — 0.20/im respectively. 

We conclude that the minority carrier transport in the neutral layers of these 
devices is efficient, and the MBE material of high quality. 

Below bandgap Eb2, the quantum well QE in sample W405 is accurately re- 
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Figure 8.4: QE of GaAs structure W407 



produced. The main source of error in this case is the reflectivity. The accuracy 
of the fit in this wavelength range suggests that the Eindhoven AR coating pro- 
cess is reproducible. Furthermore, it serves to confirm the assumption of carrier 
collection efficiencies close to 100% in the space charge region. 

Figure 18.41 shows similar data for sample W407. The device specification is 
identical to W406 except that the AlxGai_xAs is replaced by GaAs leading to 
an ungraded GaAs cell. The modelling in this case indicates minority carrier 
transport which does not differ significantly from the previous two samples. In 
general, significantly higher diffusion lengths are expected in GaAs. 

The schematic of figure 18.51 shows compositional profile of graded i region 
sample W408. In the depletion approximation, the QE below the bandgap Xb2 
is dependent only on the refiectivity and the integrated absorption. Figure 18.61 
shows the experimental and modelled QE for this sample. We note that the 
short wavelength response for photon energies above Eii2 requires slightly smaller 
diffusion lengths than samples W405-W407. The aluminium fraction Xh2 is again 
estimated at about 24%. 

The model predicts the i region contribution remarkably accurately, confirm- 
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Figure 8.5: Schematic compositional profile of graded i region sample W408 
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Figure 8.6: QE of graded i region sample W408 
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Figure 8.7: Schematic compositional profile of graded i region sample W409 

ing the reliability of the absorption coefficient parametrisation reviewed in chapter 
Hlfor aluminium fractions between and, in this case, 24%. 

Figure 18.71 shows the schematic compositional profile of the last Eindhoven 
sample W409. This sample is similar to W408 except that the i region grades are 
sharper, and a 0.2/im GaAs layer is grown in the centre of the i region. Modelled 
minority carrier transport in the neutral layers is similar to W408. The i region 
QE, for photon energies below Eb2, is again satisfactorily reproduced. 

Figure [H3] shows a schematic compositional profile of the GRINCH laser struc- 
ture NRC1038. This structure has been investigated by colleagues as a promising 
structure for a high Voc QWSC structure since studies have shown that the dark 
current in such structures is very low. The measured and modelled QE data for 
this structure are given in figure 18.101 

The p layer of this structure is very thick. It consists of a heavily doped graded 
layer starting at Xbi = 25% at the surface and rising to a nominal Xb2 = 70% 
at a depth of 120nm. This is followed by 1.484/im of AlxGai_xAs of nominal 
composition Xb2. Only a negligible amount of photons with energy greater than 
the higher p layer bandgap Xi,2 reach the intrinsic region. The p layer QE is 
negligible because of the very high aluminium fraction in this layer. The associ- 
ated low minority electron diffusion length in this case prevents photogenerated 
minority electrons from reaching the depletion layer. A sharp cutoff is therefore 
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Figure 8.8: QE of graded i region sample W409 
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Figure 8.10: QE of GRINCH structure NRC1038 

observed in the QE at energy £'^2 corresponding to a photon wavelength of about 
500nm. The main p layer composition can be estimated from this cutoff, and is 
of the order of Xh2 ~ 81%. 

The i region consists of seven different sections, as illustrated in the schematic. 
These include both ungraded sections with composition and grades from 
to Xbi. The centre of the i layer is GaAs. The Ino.2Gao.8As well at the centre 
of the i layer is modelled as an equivalent thickness of GaAs since the model is 
not equipped with an InGaAs quantum well absorption routine. The aluminium 
fraction X^i ~ 25% is assumed accurate, since we have no sharp threshold from 
which to make a rehable estimate. The QE cutoff at the nominal GaAs band edge 
occurs at a slightly higher photon wavelength than expected, and corresponds to 
AlxGai_xAs with a composition of 2% Al. 

Overall, the QE of this GRINCH laser structure is well reproduced by the 
model. Since the QE of the i layer is critically determined by the amount of 
light transmitted through the p layer, this fit confirms that the absorption of the 
graded section at the front of the p layer has been accurately modelled, with the 
proviso that the lower aluminium fraction Xbi is assumed to be 25% as specified 
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in the growth menu. 
8.2.1 Conclusions 

The series of Eindhoven samples presents a number of design features which have 
allowed us to test the reliability of the theoretical model for complex structures 
grown on a different MBE system. In particular, the graded i region samples 
allow us to test the accuracy of the depletion approximation and of the inte- 
grated absorption. The QE in both long and short wavelength ranges has been 
accurately modelled with typical values of the transport parameters. 

These structures also demonstrate high diffusion lengths for both electrons 
and holes, albeit at relatively low aluminium fractions. Finally, the interface 
recombination parameter S is negligibly small for these samples, despite the 
rather large compositional change from ~ 20% to ~ 67% between the p layer 
and the window. 

The GRINCH laser structure NRC1038 has been used to test the model in the 
case of a complex structure where the composition changes rapidly and where the 
QE is essentially determined solely by the integral of the absorption coefficient 
over the cell. 

8.3 Ungraded QWSC and p - i - n samples 
8.3.1 MOVPE samples 

All MOVPE samples were grown with nominal aluminium fractions of about 
Xb2 = 30%. They conform to the same basic design, which consists of an 80% 
aluminium window layer of thickness 0.02/im, a 0.15/im p layer, carbon doped 
at 1.34 X lO^^cm^^. The n layer consists of 0.6/im Si doped at 6 x 10^''cm~'^. 

No direct reflectivity measurements were available for the samples presented 
in this section. We therefore use the mean SiN data mentioned earlier in section 

Em 

The composition of samples frequently departs from the specifications by a 
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Figure 8.11: QE modelling and data for nominally 30% aluminium 30 well sample QT468a 




Figure 8.12: QE modelling and data for QT468b, the p - i - n control to QT468a 
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Figure 8.13: QE modelling and data QT528, which is a repeat of QT468a with a cleaner Ga 
source 

few percent. Experimental QE and modelling for the 30 well QWSC QT468a 
and its control sample QT468b are shown in figures 18.111 and 18.121 The nominal 
aluminium fraction in both cases is = 33%. Modelling however suggests the 
rather high value 42% for QT468a and 36% for QT468b. The value of 42% for 
QT468a is, as we shall see in section l8T6] supported by modelling the Fabry-Perot 
peaks of mirrored samples. 

The fitting for QT468a suggests a p layer diffusion length of L„ = 0.075/im. 
The n layer QE is negligible. The modelling uses a negligible surface recombina- 
tion parameter, which is set to 0. 

Sample QT468b shows a shortfall in QE near the band edge. This is not 
ascribed to a large n layer QE, given that a value of Lp much larger than L„ 
would be required to raise the theoretical QE in this wavelength range to the 
experimental level. The shortfall may be due either to an increased effective i 
region thickness or to the error introduced by the use of the mean SiN reflectivity. 
Finally, it may be due to the uncertainty in the absorption coefficient which was 
mentioned in chapter HJ The error in modelling however is not serious, in view 
of the ±5% error on the experimental data. 
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Finally, figure 18.131 shows modelled and measured QE for sample QT528, 
which is a repeat of QT468a but with a cleaner Ga source. The composition is 
modelled as Xf,2 = 40%. The modelled diffusion length is significantly shorter 
than QT468a, at Ln = 0.04/im. For material with these very low diffusion 
lengths, minority carrier loss at the front of the cell is dominated by bulk re- 
combination. The modelling becomes insensitive to the surface recombination 
parameter S, which cannot therefore be ascertained with any degree of confi- 
dence. It is therefore set to zero. 

The model significantly underestimates the experimental QE at wavelengths 
between about 550nm and 650nm. This may again be due to the assumption of 
a mean SiN reflectivity. Possible departures from sample specification might also 
account for the shortfall in theoretical QE but cannot be established whithout 
further characterisation, for example SIMS. 

The modelling however suggests that minority carrier transport is much lower 
than that observed in the Eindhoven samples. The diffusion lengths are slightly 
lower than the values suggested by Hamaker, and vary significantl between growth 
runs. The surface recombination parameter in all three samples was negligible. 

8.4 MBE samples 

The samples U2027 to U2032 and U4036 mentioned in sect ion [?. 1 . 1 1 were modelled 
with a view to establishing how the minority carrier transport depends on the 
base aluminium fraction Xb2 for MBE material. The samples all consist of a 67% 
aluminium window of width 0.02/im, a 0.15/im p layer Be doped at 2 x lO^^cm"^, 
a 0.481/im i layer, and finally a 0.6/im n layer Si doped at 6 x lO^^cm"^. 

The samples took the form of ungraded QWSC devices and their p - i - n 
controls. All samples were grown in the same growth run with the exception of 
the 20% p - i - n sample U2028 which was initially grown with a GaAs i region 
instead of the intended AlxGai_xAs. This device was repeated in a subsequent 
growth run as sample U4036. 

Experimental QE and modelling data for the AR coated devices processed 
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from these samples are presented in graphs 18.141 to I8.20[ 
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Figure 8.14: QE modelling and data for nominally 20% aluminium 30 well sample U2027 

8.4.1 20% samples 

Samples U2027 and U2028 are accurately modelled with aluminium fractions of 
Xh2 = 20%. U4036 however has a composition close to 23%. 

Although the modelling inaccuracies for these samples are of the order of 
the experimental error, the fitting procedure for these samples presents some 
problems concerning choice of parameters. The modelled QE peaks at a shorter 
wavelength than the data. This could in principle be modelled more closely by 
increasing Lp to boost the n layer near the band edge, at about 750nm, and 
reducing L„. The short wavelength reponse is poor, and sets a limit on the 
maximum value of L„. Since, in general, holes are less mobile than electrons, the 
electron diffusion length should be greater than the hole diffusion length. The 
Hamaker parametrisation for example suggests that L„ is about twice as large as 
Lp for all aluminium fractions. 

We cannot in principle rule out that Lp may be greater than for example 
due to particular growth or geometrical considerations. In the absence of more 
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Figure 8.15: QE modelling and data for U2028, which consists of 30% p and n layers with a 
GaAs i region 
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Figure 8.16: QE modelling and data U4036, which is a p - i - n repeat of U2028 with a base 
aluminium fraction of Xb2 = 23% 
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detailed information, however, we choose to set a hmit on the maximum value of 
Lp, which is that it should be less than or equal to 

The short wavelength response then poses a problem because a value of L„ 
high enough to fit the band edge causes the model to overestimate the QE at 
short wavelengths, even though we use a quasi-infinite recombination parameter. 

The fits presented in figures 18.141 to 18.161 are the result of a compromise be- 
tween these two constraints. The resulting diffusion lengths of Lp and L„ ~ 
O.lOyum are smaller than those given by the Hamakers parametrisation, which 
predicts L„ ~ 0.39/im and Lp ~ 0.20yum. 

Finally, the quantum wells in the QWSC sample are modelled as 70A, which 
is significantly smaller than the nominal 87A. Furthermore, we note that the 
model overestimates the QE contribution in the first quantum well continuum, 
at about 800nm. The misfit in this wavelength range is not serious, however, and 
is of the order of the experimental error of ±5% in QE. 

This may be interpreted as indicating a loss mechanism in the quantum wells, 
possibly due to the increased recombination at the interface between quantum 
wells and the barriers. Alternatively, dopant diffusion, combined with a diffusion 
length which is closer to the Hamaker value would reduce the QE at short wave- 
lengths and increase it near the band edge Ef,2. More characterisation is required 
to establish whether or not this is the case. 

8.4.2 30% samples 

Figures IHTTl and KTEl show QE data and modelling for U2029 and U2030. The 
true aluminium fractions are estimated at Xb2 = 33% and Xb2 = 36% respec- 
tively. The n layer contribution of the QWSC device appears negligible, and no 
reliable estimate of Lp is possible. The recombination parameter S is modelled 
as small but finite. L„ is then the single significant fitting parameter, and is 
modelled as about L„ ~ 0.07/im. This compares favourably with the Hamaker 
value of 0.053/im. 

The quantum well QE is again overestimated by the model. This may be 
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Figure 8.17: QE modelling and data for nominally 30% aluminium 30 well sample U2029 
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Figure 8.18: QE modelling and data U2030, the p - i - n control to U2029 
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ascribed to the same reasons as were given for U2027 in the previous section. 

For the p - i - n samples U2030, the n layer QE is again negligible. The 
surface recombination parameter in this case appears effectively infinite, although 
the predicted diffusion length is identical to the QWSC sample. Problems in the 
fit are similar to those seen for sample U4036, but are again comparable to the 
experimental error. 

The estimate of L„ in this case is more reliable than in the case of the 20% 
samples because the fitting procedure is simpler and the resulting match between 
theory and experiment better. 

8.4.3 40% samples 

The 40% samples U2031 and U2032 are mutually more consistent, and are shown 
in figures [8. 191 and [8. 201 Both have a modelled aluminium fraction of 47% instead 
of the nominal 40%. The n region QE is again vanishingly small. S in both cases 
appears too low to model, and is set to zero. 

The diffusion lengths in this case are modelled as L„ ~ 0.040/im for the 
QWSC and L„ ~ 0.055/im for the p - i - n sample. Hamaker in this case suggests 
0.02//m. 

The quantum wells in this case fit very closely in the first continuum. This 
however contradicts the possible explanation for an overestimated well QE in 
the 20% and 30%, since the interface recombination in a 40% sample should in 
principle be higher. 

8.4.4 Diffusion length tabulation 

Figure 18.211 shows the modelled value of L„ as a function of aluminium fraction 
Xb2 for QWSC and p - i - n MBE grown ungraded structures of the previous 
section. Also shown is the Hamaker parametrisation [24] for the relevant doping 
level. 

The modelled values show a weaker dependence on aluminium fraction than 
the Hamaker parametrisation. The modelled values decline slowly with increasing 
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Figure 8.19: QE modelling and data for nominally 40% aluminium 30 well sample U2031 
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Figure 8.20: QE modelling and data U2032, the p - i - n control to U2031 
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Figure 8.21: Modelled minority electron diffusion length _L„ for p - i - n and QWSC samples, 
compared with Hamaker's parametrisation for the MBE doping level doping levels 

Xb2. These values define the piecewise hnear parametrisation which is used to 
define the position dependence of L„ in the graded structures of the section 18.51 
We conclude however that the transport parameters derived from the QE 
spectra systematically indicate poorer minority carrier transport efficiency than 
has been reported in the highest quality AlxGai_xAs. Bachrach[43j. for example, 
has reported diffusion lengths of the order of 2.3/im in p doped AlxGai_xAs with 
a composition of x = 36%, although at a lower doping level of A^^ ~ lO^^cm"^. 
The review by Ahrenkiel |15] shows that diffusion lengths greater than those 
given by the Hamaker parametrisation are possible. 

8.4.5 Conclusion 

In this section we have modelled a selection of samples from the two available 
growth systems. The resulting p layer diffusion length has in all cases proved 
smaller than the p layer width. Although the MOVPE material seems to have 
more variable minority carrier transport properties, the values obtained are close 
to those of the MBE material. A small mismatch between theory and experiment 
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has been observed in a number of samples. This may be due the uncertainty in 
the absorption coefficient mentioned in chapter HJ although other influences such 
as the reflectivity, depletion widths and i region width cannot be ruled out. 

The resulting tabulation of the electron minority carrier diffusion length L„ 
is more weakly dependent on the aluminium fraction than suggested by the 
parametrisation due to Hamaker |24j. The results further indicate that the mi- 
nority carrier transport efficiency in both MOVPE and MBE material is relatively 
poor compared to values in the literature. 

The surface recombination parameter S is low in most of the samples we have 
studied. Notable exceptions are the MBE grown samples with low aluminium 
composition and one p - i - n sample from the same growth system. These samples 
used a quasi-infinite S. Supplementary characterisation is required to establish 
whether this is due to innacurate modelling of the minority carrier transport near 
the surface of the cell, or to geometrical and compositional factors. 

8.5 Compositionally Graded p layer samples 

In this section we present data and theoretical modelling of the QE of graded 
p layer samples grown by MOVPE and MBE. We will quote estimates of the 
short circuit current density J^c calculated from the experimental QE data for 
an illumination of AMI. 5. These figures do not include contact shading. 

In the case of MBE samples, we will use the diffusion length parametrisation 
based on the modelling of ungraded samples. In the case of MOVPE samples, 
we do not have as wide a range of samples with different aluminium fractions. 
The previous section however has demonstrated that the diffusion length in both 
growth systems is similar. We shall therefore use the same table of L„ as a 
function of aluminium fraction for MOVPE material. 

8.5.1 MOVPE Samples 

The graded p design was initially tested by graded p - i - n device QT601c and 
its control structure QT601a, which is identical but has no grade. The graded p 
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Figure 8.22: Experimental and modelled QE of ungraded MOVPE p - i - n QT601a 
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Figure 8.23: Experimental and modelled QE of graded MOVPE p - i - n QT601c 
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in sample QT601c is 0.15/im thick and is graded from Xbi = 67% to Xb2 = 33%, 
and has a window of thickness 0.02/im with a composition of 80% aluminium. 
John Roberts however has indicated that growth problems cast some doubt over 
the real value of the aluminium fraction X^i at the top of the grade. The QE 
data and modelling for these samples is shown in graphs 18.221 and 18.231 

The fit for both samples requires a quasi-infinite recombination parameter 
and Ln and Lp ~ 0.06yum. The fitting problems here are again reminiscent of 
samples such as the 20% ungraded devices of the previous section, but are more 
significant. A larger value of Lp would again improve the QE fit, but cannot be 
justified without more information. 

The QE of graded sample QT601c suffers from a similar problem near the 
band edge. The short wavelength QE however is accurately modelled with Vx = 
— 1.5, versus the Hamaker value of -4.2. The grade in this case enhances the 
short circuit current in AMI. 5 from Jgc = 8.06mA/cm~^ for QT601a to Jsc = 
10.35mA/cm~^. This represents an overall increase in projected Jsc of 28%. 

A set of three samples was subsequently grown to investigate the performance 
of a grade in MOVPE QWSC designs. The first of these is QT640a, which 
is similar to QT468a but for a grade in the p layer with nominal composition 
Xhi = 67% and Xh2 = 33%. Sample QT640b is identical to QT640a and has the 
same grading scheme, but has a wider i region of 0.78/im containing 50 wells. 
Finally, QT641 is a p - i - n control to QT640b but does not contain a grade. 
This last p - i - n sample was grown to compare material quality between the two 
growth runs. 

The experimental QE and modelling for QT641 is shown in figure 18.261 The 
experimental QE is very poor at short wavelengths, indicating poor minority 
carrier transport. The modelled p layer QE is dominated by bulk recombination 
losses, and the parameter S is not significant in the fitting. It is therefore set to 
zero. Finally, the n layer QE is negligible. 

The two QWSC samples bear out the conclusion that minority carrier trans- 
port efficiency for these samples was lower than for the previously grown samples 
QT601a and QT601c. Although the grades in both QWSC samples are of iden- 
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Figure 8.24: QE of 30 well QWSC QT640a, with p layer grading ranging from 30% to 67% 
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Figure 8.25: QE of 50 well QWSC QT640b with p layer grading ranging from 30% to 67% 
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Figure 8.26: QE of QT641, the p - i - n control to 50 well QWSC graded 30-67% 

tical design to that of QT601c, they perform more poorly. 

The modelhng and data for these samples is shown in figures 18.241 and 18.251 
The modelling assumes a mean SiN reflectivity. The diffusivity field for QT640a 
was = —3. The p layer QE of sample QT640b in particular is seriously 
degraded, and is modelled with = —5.5. The rather large value of for 
QT640b is expected to be at least partly due to the assumption that the L„ 
tabulation based on the MBE samples of the previous section is over-optimistic 
for these samples. 

Figure 18.271 compares the QE of QT640a with QT468a, which has the same 
design but for the absence of the grade. The QE enhancement at short wave- 
lengths is too small to be significant. We can state however that the grade has 
at the very least maintained the p layer QE in poorer material. This does not 
take into account the different reflectivities in the two samples however, which 
was better in the ungraded sample. 

At intermediate wavelengths, a significant QE enhancement is observed be- 
cause of increased light transmission to the i layer. 

Comparison of short circuit current densities gives Jsc = 14.94mA/cm~^ and 
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Figure 8.27: Comparison of thirty well sample QT468a and similar structure QT640a which 
has a graded p layer 

Jsc = 14.02mA/cm^^ respectively. This represents an enhancement of 7% in Jsc 
for the graded sample, despite the inferior material carrier transport in sample 
QT640a. We conclude that these samples illustrate the principle of a grade in 
QWSC samples, despite the fact that the ungraded benchmark sample benefits 
from superior minority transport properties. 

Finally, we present experimental and theoretical data for sample QT794a 
which cannot be explained by the model. This sample has been mentioned in 
section [7?T] and was shown to have an high Nbg- 

This sample was designed with a reduced p layer thickness of 0.1/im. The 
grading scheme was aimed at maximising the p layer response by choosing a 
moderate grade ranging in composition from X^i = 40% to Xf,2 = 30%. This 
is at the lower end of the optimised theoretical grading predictions presented 
in section I5.2.4[ and is intended to ensure good minority carrier transport by 
reducing the amount of aluminium towards the front of the p layer. We shall 
present results in the MBE section below which demonstrate successful grades 
over this compositional range. Furthermore, the results of section 15.2.51 were 
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Figure 8.28: MOVPE version QT794a of the optimised 50 well graded QWSC design 
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Figure 8.29: SIMS measurements of C doping and aluminium composition in sample QT794a 
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exploited and the window in this sample removed. 

The results are shown in figure IH2S1 The modelling benefits from a measured 
reflectivity. The MIV data indicated that the QE measurement had to be per- 
formed at at least -lOV in order to ensure full depletion of the device. We note 
that the experimental QE between about 500nm and 620nm is approximately 
equal to the i region QE alone, suggesting the the p layer QE is essentially zero 
in this range. It is impossible to model the p layer contribution consistently under 
this constraint, particularly since the short wavelength response is not negligi- 
ble, indicating that at least some minority carriers from the front of the cell are 
collected. 

SIMS and TEM measurements were performed on this sample and its MBE 
grown companion U6041, which will be presented in the following section. The 
TEM measurements confirmed that the dimensions of this sample were correct. 
The compositional SIMS measurements had to be calibrated by inferring the 
aluminium content from QE spectra. The resulting aluminium profile is shown 
in figure 18.291 The depth calibration used by Federico Caccavale in the SIMS 
results was independently verified by the TEM measurements. The mean sum 
of the well and barrier thickness was estimated as (15.3nm ± 0.5nm), which is 
experimentally indistinguishable from the nominal value of 14.7nm set down in 
the sample description. 

8.5.2 MBE 

For MBE material, it was decided to continue the study of QE as a function 
of base aluminium fraction Xb2 for QWSC graded structures. The samples are, 
in the first instance, based on those studied in section 18.41 They are all Be 
doped at a level of A^^ = 1.34 x lO^^cm"^ in the p layer and Si doped with 
Nd = 6.00 X 10^''cm~^ for the n layer. The main design characteristics of the 
samples studied here are summarised in table 18.11 No refiectivity measurements 
could be performed on the samples U4034-U4035 due to the absence of refiectivity 
measurement samples in the processing run. 
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Sample 


U4032 


U4033 


U4034 


U4035 


U6041 


AMiidow (/////) 


0.03 


0.03 


0035 


0.035 





Xbi (%) 


40 


40 


67 


67 


40 


Xk2 (%) 


20 


20 


30 


40 


30 




0.81 


0.481 


0.481 


0.481 


0.81 


Wells 


50 


30 


30 


30 


50 



Table 8.1: Principle design parameters of MBE grown graded p QWSC samples 




Figure 8.30: Comparisons of QE for QWSC samples with and without grades, for a variety 
of grading schemes. Also shown is comparison of similar graded QWSC structures grown by 
MBE and MOVPE 
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Figure 8.31: MBE 50 well QWSC U4032, graded 20-40% 

The four graphs of figure 18.301 compare the experimental QE of ungraded 
QWSC structures U2027, U2029 and U2031 with their graded counterparts. En- 
hancements in QE are observed in all cases. The increases in predicted Jsc are 
21%, 24% and 11% for the samples with X^a = 20%, 30% and 40% respectively. 
The last figure is lower because of the degraded QE in the wells. 

Figure 18.301 also shows a comparison between the MOVPE sample QT640a 
reviewed in the previous section and U4034, which has similar specifications. 
Although U4034 has a thicker window than QT640a, it nevertheless has a signif- 
icantly better short wavelength QE. Comparison of modelling results for these 
two samples indicates a p layer QE about three times better for the MBE sample. 
This confirms the conclusion reached in the previous section regarding the poor 
material quality in the graded MOVPE samples. 

Experimental and modelled QE results for these samples are given in graphs 
18.331 to 18.341 Modelling is accurate for U4032 and U4034, but shows significant 
deviations for U4033 and U4035. 

Sample U4032 requires a diffusivity field Vx =-3 and shows a slight shortfall 
in the wells, although this is not significantly greater than the experimental error. 
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Figure 8.32: MBE 30 well QWSC U4033, graded 20-40% 




Figure 8.33: MBE 30 well QWSC U4034 graded 30-67% 
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Figure 8.34: MBE 30 well QWSC U4035 graded 40-67% 

Sample U4033 presents a serious shortfall in modelled QE near band edge Eb2. 
This has previously been observed in several MBE and MOVPE samples such as 
U2027, U4036 and QT601c. As in the fitting carried out for these samples, the 
value of Lp required to fit the QE more accurately in this wavelength range is 
too great to provide a realistic explanation without direct measurement of this 
parameter. We have therefore used the value of Lp found to fit previous MBE 
samples with aluminium fractions of 20%, which is Lp ~ O.l/xm. The diffusivity 
field in this case is Vx=-2.5. 

Sample U4034 (figure I8.33P is modelled accurately with a diffusivity field 
T>x =0. In theory, this indicates a constant diffusivity in the p layer. In practice, 
however, this is probably due to the assumptions we have made regarding the 
dependence of L„ on aluminium fraction and a reflectivity which is given by the 
average SiN data. 

Sample U4035 (figure I8.34p is accurately modelled in the short wavelength 
regime with =0. This is ascribed to the same reasons given for U4034. For this 
sample, however, there is a clear misfit in i region contribution. Devices for this 
sample were generally poor, and the QE data noisy. The correct attribution of the 
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Figure 8.35: Comparison of MOVPE and MBE versions of an optimised 50 well design 
QT794a and U6041 with graded p and no window 

i region shortfall requires reprocessing, and could be assisted by supplementary 
characterisation such as SIMS measurements. The short wavelength response 
however suggests that the p layer contribution is being accurately modelled. 

Finally, figure [8735] shows optimised sample U6041. The concepts underlying 
the design of this sample are identical to MOVPE grown QT794a and were out- 
lined in section 18.5.11 The figure also shows the QE of the MOVPE device. As 
in the case of QT640a and U2029, the MBE device again performs better. 

Initial modelling of the experimental data with nominal parameter values 
showed similar problems to those experienced for QT794a. In particular, the 
modelled i region QE was greater than the experimental QE between about 
550nm and 600nm. This suggests departure from design characteristics. Fur- 
thermore, the model significantly overestimates the QE at short wavelengths, 
indicating an impaired p layer QE. 

Extensive characterisation was carrier out on this sample. Calibration of the 
SIMS depth scale was verified by TEM, as for sample QT794a in previous sec- 
tion 18.5.11 TEM data indicated that mean sum thickness of barrier and well 
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was (14.1nm ± 0.5nm), with the error estimated from the resolution of the mi- 
crograph. The specification was 14. Trim, which is not significantly greater. The 
total thickness of the multiple quantum well system confirmed the depth cali- 
bration of the SIMS data presented in figure 18.361 The graphs show only the 
most relevant data, which are the aluminium fraction and Be doping for the first 
0.6/im of the sample, including the 0.04/im GaAs capping layer. Lower resolution 
measurements covering the entire sample were also supplied. 

Federico Caccavale has confirmed that the aluminium SIMS signal in counts 
per second will be linearly dependent on the aluminium composition for the range 
of aluminium fractions in these samples. Modelling of a mirrored device processed 
from the U6041 wafer will be presented in section 18.61 This has allowed us to 
fix the lower aluminium fraction at Xf,2 = 24%. The linear relationship between 
SIMS signal and aluminium fraction, together with the value of Xf,2, enable us 
calibrate the SIMS mass spectrometer signal. 

The resulting aluminium compositional profile is shown in the lower part of 
the figure. We note first that the aluminium composition peaks at only 29%. 
Christine Roberts who grew this sample has confirmed that grades over thin 
layers are technologically difficult. This is a consequence of the growth scheme 
for graded structures, which treats the grade as a series of ungraded strips grown 
close together. In practice, a smooth grade is obtained by ensuring that the 
aluminium source is never in equilibrium, and has a constantly increasing output 
flux. 

If however the nominally ungraded strips are too thin, the aluminium source 
has too little time to ramp up to the higher temperature required for the ma- 
terial with a higher aluminium fraction. The result is a lower final aluminium 
composition than expected, as seen in this sample. The grading technique can be 
improved by characterisation such as the SIMS data presented here, which can 
help in determing an optimum strip width resulting in both a smooth grade and 
an accurate target aluminium fraction X^i. 

Furthermore, the unexpectedly slow decrease in aluminium at the front of the 
p layer has resulted in a positive bandgap gradient with increasing distance at 
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the front of the p layer. The origin of this feature is not clearly understood. 

This rounding off of the aluminium profile towards the top of the sample 
results in a thin layer with a bandgap gradient pointing towards the front of the 
cell. The resulting field will tend to drive minority carriers generated towards 
the front of the cell away from the depletion region, thereby degrading the short 
wavelength QE. This unforeseen effect is a result of removing the window layer. 
A similar effect in a sample with a window would not cause this problem, because 
the step in the conduction band at the interface between the window and the p 
layer acts clS 8b minority carrier reflector. 

Finally, we note that the Be dopant profile extends significantly further than 
the nominal O.l/xm. Illegems ^8] studied Be diffusion in MBE GaAs and Alx- 
Gai_xAs. He concluded that this dopant was promising insofar as it allows high 
doping levels (~ lO^^cm"^) whilst maintaining good surface morphology in the 
MBE material. Be diffusion coefficients are quoted in this work. Be diffusion 
has more recently been reported as a problem. Lee [59] has studied this problem 
more recently in AlxGai_xAs as a function of aluminium composition and anneal 
temperature. 

Using the correlated TEM and SIMS data, we predict an effective p layer 
thickness of 0.14/im. This fairly conservative value, is based on the fact that the 
quantum well signal fits theory accurately, suggesting that the first wells on the 
p side of the i layer are not significantly doped. We note finally that the residual 
background doping in this sample is n type silicon, and were estimated in chapter 
Has being below lO^^cm"'^ by both SIMS and the MIV technique. 

The resulting QE modelling is shown in figure 18.371 Theory reproduces the 
experiment closely with a diffusivity field of -4.5. The surface recombination 
parameter is set to an effectively infinite value to account for the drop in alu- 
minium fraction over the first few nm of the sample. The aluminium composition 
used in this case is based on the SIMS scan, and includes the ungraded 0.04/xm 
unintentionally doped layer between the p and i layers. 

The SIMS and TEM data also indicated that the n layer thickness is 0.72/im 
thick instead of the nominal 0.6/im. 
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Figure 8.36: SIMS measurements of Be doping and aluminium composition in sample U6041 



The dopant diffusion has proved crucial in understanding the QE of this sam- 
ple. The extension of the p layer and consequential reduction in i layer thickness 
explains why modelling using nominal sample specifications overestimated the i 
region QE. Similarly, the extra thickness in the p layer and lack of grade over 
nearly a third of this layer explains the reduced p layer QE. 

The final efficiency of this sample was estimated by Paul Griffin at (14.2% ± 
0.5%) under tungsten illumination corrected to AMI. 5 light levels. This is compa- 
rable with the efficiency reached by the previous best AlxGai_xAs sample G951, 
which was mentioned in section 17.1.11 There is no control sample to compare the 
estimated short circuit current enhancement for this sample. 

Despite the departure from device specifications, we have shown that the 
model has adequately reproduced data given corrected device specifications. This 
increases confidence in the ability of the model as a diagnostic tool. 
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Figure 8.37: MBE 50 well QWSC without a window graded from Xbi29% to Xb2 = 24%. 

8.5.3 Conclusions 

The study of a wide range of grades has demonstrated their capacity to increase 
the Q£; of p - i - n and QWSC samples grown both by MBE and by MOVPE. 

The case of sample U6041 has served to increase our confidence in the model 
both as a diagnostic tool and as a basis for further optimisations. 

A large amount of detailed information about each sample is required. Since 
the advanced experimental characterisation cannot be carried out as a matter of 
course, we are frequently forced to make assumptions about composition. The 
modelling has furthermore had to use transport data extrapolated from ungraded 
samples. Finally, reflectivities are not always available. 

As a result, theoretically well defined parameters such as the the diffusivity 
field and the recombination parameter which are used in the modelling are in 
practice ill-defined, since modelled values contain systematic errors due to the 
assumptions we have mentioned. 

Despite these difficulties, the accuracy of the compositionally graded Alx- 
Gai_xAs model has been demonstrated for a wide range of samples. In addition, 
the model can reproduce changes in QE due to a number of design changes to 
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Figure 8.38: QE of mirrored sample QT468a. Also shown for reference is the QE of an AR 
coated device without a mirror 

the original structure. These include compositional grades in neutral layers and 
depleted layers from three different growth systems. Enhancements in QE for 
samples with graded p layers have been shown in a wide range of samples. 

The diagnostic capabilities of the model have furthermore been confirmed 
by comparison with quantitative characterisation techniques. To our knowledge, 
this is the first model of compositionally graded AlxGai_xAs cells to be used to 
model experimental data. 

8.6 QE of Mirror Backed Samples 

In this section we investigate the QE enhancement in the quantum well wave- 
length range. The technique was tested on the 30 well sample QT468a. Theoret- 



ical modelling and experiment for this sample is shown in figure 18.381 together 

with the QE of an AR coated device without a mirror for comparison. The in- 
crease in predicted Jgc for this device is 27%. Main fitting parameters are given 
on the graphs. 

The wavelength dependence of the Fabry-Perot peaks is accurately modelled. 
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Figure 8.39: QE of a mirrored QT468a device demonstrating non-specular reflection from 
the back surface and QE enhancement greater than two passes. 

Variations in model accuracy over the quantum well wavelength range are visible, 
and are due to the underlying inaccuracy of the quantum well absorption model. 

Amplitudes of Fabry- Perot peaks in the experiment are variable, and generally 
lower than those predicted by the model. The broadening mechanism responsible 
for this is not well understood, but may be partly due to roughness at the surface 
of the mirror. 

Figure 18.391 illustrates a mirrored device processed from the QT468a wafer at 
the same time as that of figure 18.381 We note a complete absence of interference 
peaks in this case. Moreover, the QE enhancement is more than can be explained 
by a double pass. We interpret this as indicating that the light is not being 
specularly reflected from the back surface. The randomising effect of reflection 
from a rough mirror causes the relative phases of internally reflected light beams 
to be randomized. Since the mean path length is also increased, this results in 
both increased light trapping in the cell and the disappearance of the Fabry- Perot 
peaks. Both devices processed from the QT468a wafer are modelled with a zero 
back surface phase change. 
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Figure 8.40: Experimental QE enhancement in two mirrored QT229 devices, compared with 
a AR device without a mirror 

Figure 18.401 shows experimental data for QT229, which is a 50 well QWSC 
with a very thick p layer of nominal thickness 0.5/im. The three QE curves relate 
to an AR coated device and two AR coated and mirrored devices. We note that 
the positions of the Fabry-Perot peaks for the two mirrored samples are signifi- 
cantly out of phase. The modelling for these two samples is given in graphs 18.411 
and 18.421 The modelled phase change differs by a factor of nearly tt. This result 
is not unexpected in the light of discussions with Mark Whitehead [37] which 
indicated that the phase change is device dependent at metal- semiconductor 
interfaces. 

We note that the frequency of the theoretical peaks in wavelength is smaller 
than that of the experiment. This suggests that the model is underestimating 
the width of the optical cavity defined by the mirror and the front surface. The 
accurate fits obtained in other samples support the conclusion that the cell is 
wider than expected. 

Experimental and theoretical QE for a mirrored device processed from the 
U6041 wafer is shown in figure [H331 Only two devices of this type were processed, 
one of which was not successful. 
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Figure 8.41: QE of a QT229 mirrored device showing a back surface phase change ijjb = 
radian 




Figure 8.42: QE of a a QT229 mirrored device showing a back surface phase change ipb = —3 
radian 
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Figure 8.43: QE of a mirrored device processed frora the U6041 wafer 

Despite this, an accurate fit is obtained using the same parameters as in 
section 18.5.21 together with a back surface phase change of -2.8 radians and a 
front reflectivity factor J-" = 0.4, such as was estimated for the QT229 modelhng 
of figure 

No measured reflectivity was available for this sample. We have assumed 
a poor reflectivity, basing this on the fact that the mirrored device has a sig- 
nificantly poorer response at short wavelength than the device without a mirror 
discussed in section [8.5.21 This was clearly visible on the devices, which showed a 
green tint characteristic of high reflectivity at short wavelengths. This device was 
processed at the time as the sample whose reflectivity was used in the modelling. 

As we explained earlier, this sample was studied by SIMS and TEM. The 
wavelength dependence of the Fabry-Perot oscillations is reproduced. 

This sample again illustrates the diagnostic capacilities of the model, by show- 
ing that accurate modelling of the experimental data is possible if the structure 
is of good quality and the composition and geometry is well known. 
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8.7 Optimisations 

Optimisation of Available Material 

The model can now be used to optimise the QWSC design we have been consider- 
ing. Two approaches to this exist. The first of these is to identify design changes 
made separately to samples that have been characterised, and to combine these 
in a single optimised sample. 

This relatively conservative approach is illustrated in figure [H311 This design, 
which we will call sample A, is based on a combination of the Eindhoven AR coat, 
an 80% window layer, and the grade demonstrated in sample U2029 together with 
the back surface mirror seen in 50 well QWSC devices processed from QT229 and 
U6041. Also shown is a prediction for the correponding control without a grade 
or a mirror. 

The modelled short circuit currents in this case are 22.3mAcm~^ for A and 
and 17.0mAcm^^ for the control. This corresponds to an enhancement of 31% 
with the new design modifications when compared to the control shown in the 
figure. 

The second and less reliable approach is an extrapolation to design parameters 
which have not been observed. Figure 18.451 shows a similar comparison between 
an optimised sample B and its control. The design in this case has an Eindhoven 
AR coat, no window, and a thinner O.l^m p layer. In the case of the optimised 
sample, this p layer is again graded between X^i = 67% and Xi,2 = 30%. 

The modelled short circuit current for these predictions are Jgc = 23.6mAcm~^ 
for the optimised design and Jsc = 19.7mAcm~^ for the control, or an enhance- 
ment of 20% overall compared to the control without a grade or mirror. 

These results indicate that the more conservative approach can nearly match 
the current output of the more adventurous design. Both samples compare 
favourably with the previous best MBE sample G951 which was mentioned in 
sections 17.1.11 and I8.5.2[ and which has a theoretical J^^ = 18.2mAcm~^, com- 
pared to the experimental value of (17.5±0.5)mAcm^^ quoted in Barnham et al. 
[50] . This estimate is an independent measurement carried out at the Frauenhofer 



8. 7. Optimisations 



203 



institute in Germany. 

A prediction of the efficiency requires a model of the dark current of the celL 
This is the subject of ongoing work by collaborators. However, we can make a 
rough estimate by assuming that the new samples have the same J-'J-' = 78% and 
Voc = 1.071^ which were observed experimentally in G951. This sets a lower limit 
to the projected efficiency of optimisation A since the higher current levels in this 
cell will increase the Voc, and the similarity between G951 and design A should 
produce the same JT". With these assumptions, the final estimated efficiency for 
this sample is then 18.6%. 

The corresponding efficiency for sample B under the same assumptions is 
19.7%. These results compare well with a good conventional ungraded p - n 
cell reported by Amano et al. [73j, which had an aluminium composition of 
30% and was 14.6% efficient under AMI. 5 illumination. This cell was the result 
of extensive growth optimisation and had large diffusion lengths Lp and L„, a 
high Voc = 1.281^, a ffil factor FF = 70.7% and a short circuit current Jsc = 
16.2mA/ cm~'^. 
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Figure 8.44: Theoretical prediction for optimised sample A combining the Eindhoven reflec- 
tivity, a 0.15/im grade with Xbi — 67% and Xb2 = 30%, a 0.02/zm window with 80% aluminium, 
a high surface recombination and a back surface mirror. 
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Figure 8.45: Theoretical prediction for optimised sample B combining the Eindhoven reflec- 
tivity, a thin O.lO/im grade with Xbi = 67% and Xf,2 — 30%, no window, a negligible surface 
recombination and a back surface mirror. 

Material with More Efficient Transport 

Further efficiency enhancements are possible in material with improved minority 
transport parameters. We saw in section 15.2.41 that a graded p layer may re- 
duce the efficiency of samples with good material quality. Section 18.1.51 further 
showed that the QE of graded structures in particular is sensitive to the position 
dependence of the diffusivity. Without a systematic study of the compositional 
dependence of the minority carrier transport in such material, therefore, we can- 
not reliably make predictions of the QE and Jsc enhancements expected in this 
system. 

We therefore limit ourselves to short circuit current predictions of the best 
possible ungraded AlxGai_xAs QWSC cell, using the highest reported literature 
values for the transport parameters. We will compare these projected efficiencies 
with the results of the previous optimisation section. 

The paper by Ludowise [27] mentioned in sect ion [3 . 1 . 1 1 rep ort s very high diffu- 
sion lengths in a series of solar cell samples grown by MOVPE. These values were 
derived from internal QE spectra, using a method similar to the one described 
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in this thesis. 

An optimisation based on the Ludowise diffusion length of L„ ~ lAfim at 30% 
aluminium composition yields a value of Jgc = 25.6mAcm~^. This corresponds 
to an efficiency between 21% and 24% for low and high values of FF and Voc in 
AMI. 5. These results apply to a 50 well sample with a back surface mirror, a 
0.02fim window of composition 90% aluminium, no compositional grade and an 
optimum p layer thickness of 0.15/im. 

The Ludowise diffusion lengths, however, were estimated from a set of n-p 
algaas samples. The structure of the samples is not given in detail. The front n 
layer, however, is very thin. This indicates that the diffusion lengths are estimated 
from the QE near the AlxGai_xAs band-edge. The paper does not make clear 
how the n and p contributions to the QE in this regime were separated, since the 
methods we have used in this work do not apply to the Ludowise n-p structure. 
Loo [28] however has reported hole diffusion lengths of 2/im at an aluminium 
composition of 30% and a minority carrier lifetime of 13ns, which is close to 
GaAs. The Loo paper combined with data from Ludowise indicates that high 
minority carrier transport parameters are achievable in AlxGai_xAs. 

More readily achievable minority carrier lifetimes are reviewed in the review 
by Ahrenkiel |15]. These qualitatively agree with the values we have derived for 
the Eindhoven samples studied in section 18. 2[ and reproduce the functional form 
of the Hamaker parametrisation. Using a value of L„ ~ OAfim and a structure 
similar to that described in the previous paragraph, we find an optimum short 
circuit current of the order of Jgc = 24.2mAcm~^. This corresponds to lower and 
higher limits to the estimated final efficiency of 20% to 23%. 

The second optimisation is more likely to be achievable, since it is based 
on reproducible material quality. It should be noted that the final efficiency 
estimates are subject to large errors because of the lack of a reliable model for 
the Voc and the FF. 

We note that the estimated efficiencies for optimised QWSC structures with- 
out grades but with a back surface mirror are only about 5% to 10% greater than 
those of the optimised graded structures presented in the previous section for 
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the AlxGai_xAs material with significantly poorer minority transport efficiency 
studied in this thesis. 



8.8 Conclusions to chapter [8 

This chapter has demonstrated the accuracy of the theoretical model for the QE 
of mirrored and compositionally graded QWSCs by comparison with experiment 
for a wide range of AlxGai_xAs samples. 

Results from a set of samples grown and processed in Eindhoven University 
of Technology showed that the theoretical model is capable of reproducing ex- 
perimental data for QWSC and p - i - n structures with graded i layers. These 
samples further demonstrated accurate ffis in samples with good material. 

MOVPE and MBE structures with and without quantum wells and with un- 
graded neutral layers were characterised. These are used to establish the depen- 
dence of the minority carrier diffusion length L„ on aluminium fraction. These 
samples have been successfully modelled, and values of L„ for both growth sys- 
tems extracted. MOVPE values are more variable than MBE ones but remain 
only fractionally different to MBE ones. The fitting is most accurate for in- 
termediate aluminium fractions of Xi,2 — 30%. The model has some difficulty 
reproducing experimental QE data near the band edge Ef)2 for material with 
low aluminium fractions. At higher aluminium fractions, the well QE could not 
be modelled in one MBE sample. This is ascribed to a loss mechanism in the 
quantum wells. 

Very high values of the surface recombination parameter S were required to fit 
some samples. The origin of this problem has not been established. It is suggested 
that it may be due to shortcomings in the way the surface recombination is 
modelled. The presence of a damaged layer near the surface of the cell may 
explain this behaviour. Another candidate is dopant diffusion in MBE AlxGai_x- 
As, since this effect would extend the width of the p layer, thereby reducing the 
QE at short wavelengths. More information is required to resolve this issue. 

A wide range of QWSC and p - i - n designs with graded p layers was then in- 
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vestigated. Samples grown by MOVPE have demonstrated the QE enhancement 
but have shown significant materials problems. One optimised sample in partic- 
ular cannot be modelled despite extensive characterisation. This is ascribed to 
the poor material quality indicated by high background doping levels and a very 
low p layer QE. 

Graded structures grown by MBE have shown a clearer picture. Good fits are 
obtained overall with values of the diffusivity field which are between and 
the Hamaker value. Samples requiring Vx — are interpreted as departing from 
the assumptions we have made regarding the variation of L„ with aluminium 
fraction, and, in some instances, to the reproducibility of the SiN refiectivity. 

An optimised sample which did not initially perform as expected has been 
successfully modelled after supplementary information was provided by SIMS 
and TEM measurements carried out by collaborators. This information allowed 
us to correct the device specifications to allow for dopant diffusion in MBE ma- 
terial. This has cast some doubt over the p layer thicknesses used in modelling 
MBE samples from the same growth system. It was concluded however that the 
influence on the modelled transport parameters is less critical for other samples, 
because they have thicker p layers. This however is reminiscent of the conclusions 
we reached when discussing the surface recombination parameter. 

The optimised MBE sample was repeated on the MOVPE reactor. Experi- 
mental data from this sample cannot be modelled despite similar characterisation. 
The SIMS together with MIV measurements, however, indicate poor material 
quality. This indicates that the samples may be too highly contaminated by 
impurities to operate successfully as solar cells. 

Mirrored samples have been modelled accurately. One sample has confirmed 
that the phase change upon refiection at the back mirror is a variable parame- 
ter. A significant and variable broadening mechanism has been demonstrated by 
comparison of different devices processed from the same wafer. In the case of a 
QT468a device, this can be explained by an accidentally roughened back surface 
mirror and consequent randomising of the relative phase of the incident light. 

The spacing of Fabry- Perot peaks on the wavelength axis has been accurately 
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reproduced. This indicates that QE characterisation of mirror backed samples 
may be of use in estabhshing the optical thickness of samples. Finally, mirror 
backed samples are also useful in determining the precise aluminium composition 
Xb2 of QWSC samples, since the onset of Fabry-Perot oscillations is a sensitive 
function of the former. 

We have demonstrated that both grades and mirrors are effective in increasing 
the QE and hence the short circuit current of QWSC structures in AlxGai_xAs 
with poor minority transport characteristics. 

This information gives us sufficient confidence in the diagnostic and predictive 
capabilities of the model to suggest further optimisations. 

Optimisations using the model together with information from the experi- 
mental data were presented. The projected efficiency for the first of these op- 
timisations is based on separate successful design changes which were verified 
experimentally. This produces a projected Jsc which is significantly greater than 
published good conventional AlxGai_xAs solar cells with comparable aluminium 
fractions. 

The second optimisation makes design changes which are an extrapolation of 
experimentally measured samples, but only exceeds the conservatively optimised 
sample of the previous paragraph slightly. A rough estimate of the efficiencies of 
these projections suggests a value of the order of 19% in AMI. 5 

Finally, these optimisations apply to available material. Further optimisations 
of growth conditions are required for the QWSC application in order to produce 
material with the minority carrier transport efficiencies which have been observed 
in Eindhoven material for example. We have mentioned the work by Amano [75] 
which started with fairly poor material but produced efficient minority carrier 
transport after a detailed growth optimisation study. 

A set of calculations based on good material quality has predicted efficiencies 
of between 20% and 24% in optimised QWSC structures without grades but with 
back surface mirrors. This study shows that the material quality is the prime 
problem. A grade can compensate for poor material quality to a certain extent, 
but is not expected to improve good quality material. 
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We conclude however that the combination of a grade and a back surface 
mirror can bring the Jgc of QWSCs fabricated from poor material quality to 
within a few percent of that achievable with the best material. 



Chapter 9 



Conclusions 



This thesis has adressed the problem of increasing the quantum efficiency (QE) 
of QWSC and p - i - n solar cells in the GaAs/AlxGai_xAs materials system. 
The approach to the problem consisted of joint theoretical and experimental 
investigations. 

Chapter |3] reviewed previous models dealing with compositionally graded 
structures. These do not attempt to reproduce data. We have developed a the- 
ory based on previous work which is capable of modelling a wide range of cells 
with position dependent bandgaps in all regions of the solar cell. In addition, the 
model is designed to model solar cells with back surface mirrors. The internal 
reflections from quantum well-barrier interfaces however was not treated. 

We discussed properties of the analytical and numerical solutions for the QE. 
We concluded that the wavelength dependence of the QE enables us to estimate 
the minority electron diffusion length in certain samples. We further concluded 
that the recombination velocity in the p layer is less significant in graded p layers. 

Consequences of the modelling were discussed in chapter \5\ together with 
an analysis of the ideas underlying the optimisation program. The effect of 
compositional grades in the p type emitter layer, and the effect of adding a back 
mirror were investigated and illustrated for a range of significant cell parameters. 
Theoretical predictions of the QE enhancement were presented separately for the 
compositional grading technique and the mirror design. 
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After an introduction to the experimental techniques used in the project, the 
prehminary characterisation of our samples was presented. These results suggest 
a background doping problem in multiple quantum well material. The problem 
is more severe in MOVPE material. It was suggested that this problem is due to 
a loss of compensation between carbon and oxygen background impurities in the 
MOVPE QWSC structures. 

The effect of the anti reflection coat on solar cell MIV characteristics was 
considered. Parallels were drawn with hydrogen passivated samples examined in 
the course of a collaboration with Eindhoven University of Technology. 

A set of samples with a range of i region compositional profiles including 
quantum wells and grades was characterised and modelled. The high minor- 
ity carrier efficiency predicted in these samples indicated good material quality. 
These samples have confirmed the capacity of the model to reproduce the QE of 
a range of complex graded i layers. They also demonstrate successful modeUing 
of the QE of p and n layers with efficient minority carrier transport. 

We then examined a set of QWSC and p - i - n MBE samples with different 
aluminium fractions. A tabulation of minority electron diffusion lengths L„ in 
the p layer was estimated from the modelling. The MOVPE samples grown for 
this project were also considered and showed values of L„ which are similar to 
the MBE material. These values are significantly lower than those estimated for 
Eindhoven material at a similar composition. Furthermore, the values of L„ arc 
lower than the majority of the published data for low aluminium fractions, and 
are weakly dependent on composition overall. We concluded that the minority 
carrier transport parameters can be improved in both materials systems. 

QWSC and p - i - n samples with compositionally graded p emitter layers 
were designed and characterised. Experimental enhancements were observed in 
all these structures, except one MOVPE sample. The QE measurements were 
modelled subject to the following assumptions. We assumed that the variation 
of Ln with aluminium composition in grades was similar to that estimated from 
QE measurements of ungraded samples. We also assumed that the diffusivity 
decreases exponentially with increasing aluminium fraction. This allowed us to 
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model the solar cells with a limited number of parameters. 

The experimental QE data were successfully modelled for a range of samples 
grown by MOVPE and MBE One device which could not be modelled initially was 
characterised extensively by SIMS and TEM. This additional information allowed 
accurate modelling of the experimental QE with corrected sample specifications. 
Diffusion of the beryllium p type dopant into the i region was identified as a 
problem. With this increased knowledge of the sample structure, a good fit was 
obtained. The companion MOVPE sample proved to suffer from poor material 
quality. Despite similar extensive characterisation, this sample could not be 
modelled. 

The study of ungraded samples has shown that the material used in this work 
suffers from relatively poor minority carrier transport efficiency. 

The QE data for graded structures concluded that significant QE enhance- 
ments can be achieved by compositional grading of the p type emitter. 

A series of mirror backed QWSC samples was examined. Significant enhance- 
ments were seen in all devices. The front reflectivity was used as a fitting pa- 
rameter to adjust the amplitude of the Fabry-Perot peaks observed in the data. 
One sample which shows no Fabry-Perot peaks supports the suggestion that this 
decreased amplitude is due to imperfectly smooth back surface mirrors. 

Two devices from a 50 quantum well sample showed evidence of phase changes 
upon refiection from the back surface which differ by approximately tt. This vari- 
ation supported the use of the back surface phase change as a fitting parameter 
to model the wavelength position of the Fabry-Perot peaks 

We then showed that the QE model for mirrored samples is capable of ac- 
curately verifying both the optical thickness of the sample and the aluminium 
composition of the quantum well barriers. The origin of the overestimation of 
the Fabry-Perot peak amplitude by the model should, however, be verified. A 
calculation of the effect of fight reflection at the quantum well-barrier interfaces 
would contribute to better understanding of this point. One sample matched 
the efficiency (about 14.2% in AMI. 5) of our previous best AlxGai_xAs solar cell 
with benefltted from signiflcantly better minority carrier transport. 
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The final part of tliis work presented two optimised designs incorporating 
tliese two design features. Tlie first is based on optimisations observed sepa- 
rately in samples characterised during the course of the work. This predicted a 
theoretical short circuit current of 22.3mA/cm^. This represents an enhancement 
of 31% over a control without grade or back surface mirror. The projected lower 
limit to the efficiency in this case is 18.6% for AMI. 5. 

The second optimised device extrapolated from the previous design by thin- 
ning the p layer and removing the window layer. In this design, the theoretical 
short circuit current prediction was 23.6mA/cm^. This represents an enhance- 
ment of 20% over a control, and a final projected efficiency of about 19.7%. These 
calculations demonstrate enhancements possible relatively poor material. These 
efficiencies assume an unchanging Vqc and FF, and therefore represent only a 
lower limit. 

Preliminary calculations in section 18.71 compared the highest projected short 
circuit current obtainable from graded QWSCs made from available material, 
and ungraded QWSCs made from the best reported AlxGai_xAs. Both designs 
were coated with a back surface mirror. The QWSC structures made from the 
best reported material have projected short circuit current densities between 
24.2mA/cm^ for the most commonly reported higher limits to minority carrier 
transport, and 25.6mA/cm^ for the highest reported values. The projected effi- 
ciencies ranged from 21% to 24% for the most efficient material, and from 20% to 
23% for the standard, high quality material. These calculations suggest that the 
compositional grading and back surface mirror techniques can enhance the short 
circuit current density and efficiencies of solar cells made from poor material to 
within a few percent of cells made from the best material. All these projections 
are significantly greater than the efficient, independently verified 14.6% achieved 
by Amano et al. ^3] and mentioned in chapter 18.71 

We conclude that both optimisation techniques have been successfully demon- 
strated, and that significant enhancements in short circuit current are possible. 

Future work will investigate the optimised designs presented in chapter EJ 
The light intensity calculation should be generalised to include internal reflec- 
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tions at interfaces other than the front and back surfaces. The optimisation and 
modeUing techniques can furthermore be generahsed to other materials systems 
where material quality is an issue, and where compositional grades are attractive, 
such as CdHgTe or CuInGaSe. 



Appendix A 
Analytical Solutions 



The analytical solutions for the quantum efficiency of the neutral regions were 
described in section |3l This appendix discusses exact and approximate analytical 
solutions to the differential equation for the minority carrier concentration 13.231 
of chapter [31 and uses the same definitions. 

The solutions to equation 13.231 are dependent on three materials and device 
dependent transport parameters, which are the diffusion length L„, the recombi- 
nation velocity Sn and the diffusivity The number of free parameters must 
be reduced to one before experimental data can be modelled with confidence and 
optimised structures designed. 

This appendix examines these analytical solutions in a number of limiting 
cases in order to determine conditions in which the number of free parameters 
may be reduced in order to allow single parameter fits. 

The discussion quotes the results given in section |3l and is similarly restricted 
to the case of a p type layer for brevity. Where illustrative examples are given, 
the calculations apply to a p type layer with an aluminium fraction of 30%, which 
is the design most commonly found in our structures. 

The ffist section briefly restates the solutions obtained in section [3] and states 
the exact general solution subject to the usual boundary conditions. 

The second section deals with an ungraded cell without fields in the neutral 
regions. We examine the behaviour of the analytical solution for the QE in a 
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number of limiting cases. 

We then consider a cell with constant electric field in the neutral regions, but 
again with constant minority carrier transport parameters, and compare these 
results with the results for the cell without fields. 

A.l General Form of the Analytical Solution 

The general solution to equation 13.231 can be written 

n{x) = aiUi + a2U2 + f/3 (A.l) 

where functions Ui, U2 and U3 are described in the general case in section [231 We 
saw in section 13.7.21 that the effective field experienced by the minority electrons 
in the conduction band can be expressed in per metre as 

In order to further simplify the final expression for the photocurrent, we introduce 
a surface boundary condition parameter S which includes both the recombination 
velocity and electric field at the front of the layer 

S = 2Sf-^ (A.3) 

The surface boundary condition is then of the form 

Sn{0) + n{0) = (A.4) 

The boundary condition at the interface between the neutral region and the 
space charge region in the depletion approximation is that the excess minority 
carrier concentration is zero 

n{z) = (A.5) 
where z is the position of the depletion region, measured from the front of the p 
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layer. 

Substituting for the expression for n and n' into these two boundary conditions 
allows us to express the general form of the constants ai and 02, which are 

^ _ -Usiz)iUm + ^^2(0)) + U2iz)iUm + ^^3(0)) .... 

' ui{z){um+su2m-u2{z){um+suM) 



^ U3iz)iU[iO)+SU^m - u^iz){u;,{o)+susm 

' u,{z){um + su^m - U2{z){um + sum) ^ ' 

These expressions define the general solution to equation 13.231 It proves difficult 
however to find an analytical expression for functions Ui and U2 in most cases. 
We now consider the two simplest cases where this is possible. 



A. 2 Ungraded p Layers 



For an ungraded layer, the transport parameters are constant and the pseudo- 
field field is zero. We have seen in section 1X71 that the functions Ui and U2 for 
an ungraded p layer are 

Ui = e+^^ 
Uo = e-^"" 



(A.8) 



where k is 

k = l/Ln (A.9) 
The generation function g defined in section 13.6.21 takes the form 

Oie~'^^ 

g{x,X) = —— (A.IO) 

Function f/3 and its differential U^ are then found from equation 13.291 and take 
the form 

TT — ae~"^ 

^3 - D„{a^-k2) (A.ll) 
TTl a e 

We note that both U3 and f/3 are inversely proportional to the diffusivity Dn- 
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Equations IA.6I and IA.7I for this case yield the following forms for ai and 02: 



-ae 



a2 



'-{S + k) -e+''%a^ 



aS) 



- P) [e+f^'iS -k)- e-^'{S + k)] 



(A.13) 



Equations IA.8I to lA. 131 define the excess minority carrier concentration n. We 
can then write down the quantum efficiency using equations 13.11 and 13.201 from 
chapter |3l We choose to express the QE for a unit incident light ffux in terms of 
the three following functions 



/i 



ake 



+kz 



e~''%S-k)+e~'''{a-S) 



-kz 



(S-k) 



-kz 



(S + k) 



(A.14) 



/2 



ake~ 



k"^ — 



e-'^^S + k) + e+^''{a-S) 



(A.15) 



— k"^ 

such that the fractional quantum efficiency is 



(A.16) 



QE{\) = f, + f, + f, (A.17) 

We note that all three functions go through a singularity at k = a. It can be 
shown that /i is continuous. Although functions /2 and /a are both discontinuous, 
a numerical study shows that the sum /2 + /a is small compared to /i and 
continuous. This is illustrated in figures lA.ll and IA.2I for high and low values of 
the diffusion length and high recombination velocity. 

The ffgures show that /i dominates the general form the QE, whereas the 
sum /2 + /a constitutes a correction which has a varying signiffcance in different 
physical limits. 
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Figure A.l: Components of the quantum efficiency for small diffusion length 
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Figure A. 2: Components of the quantum efficiency for large diffusion length 
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A. 2.1 Case of Dominant Bulk Recombination 

An examination of equations IA.14I to IA.15I shows that the surface boundary 
condition parameter S occurs in terms S ± k and in the term — aS. For 
absorption coefficients greater than ~ 10^, the latter is dominated by a^. This 
is true in particular of the absorption coefficient from 400nm to approximately 
GOOram in the typical 30% AlxGai_xAslayer we are considering as an example. 

We will see in the next section that for longer wavelengths and lower absorp- 
tion coefficients, the effect of the recombination velocity is small if the diffusion 
length is small compared with z. 

We find that the QE will be insensitive to both the recombination velocity and 
the diffusivity if \ S ±k |~ A;. Expressed in in terms of the materials parameters, 



this gives the condition 



f » §^ (A.18) 



This is trivially true if the recombination velocity is very low in which case the 
only free parameter is L„. 

More significantly, this condition will hold if the diffusion length is very small. 
In this case, the number of minority carriers lost at the surface is dwarfed by 
recombination in the bulk, and L„ is the main parameter determining the QE. 
Further discussion of the balance between these parameters is given in chapter [HI 

For S much greater than k and a, functions /i and /2 both tend to limits 
which are independent of S. This is equivalent to a vanishingly small minority 
carrier concentration at the front of the p layer. We conclude that the QE is 
again independent of the magnitude of Sn and Dn for S much larger than k and 
a. 



A. 2. 2 QE Near the Band Edge 

This section investigates simple limits for the behaviour of the quantum efficiency 
near the band edge. In this regime, the absorption coefficient is small compared 
to z and the generation rate is approximately constant. Generation rates in this 
regime are illustrated in chapter O 
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Small Diffusion Length 

If tfie diffusion length is small compared to the p layer thickness z, the product 
kz becomes large. Function Ui is much greater than 1/2'- 

This condition is usually true in the QWSC designs we consider, function Ui 
being typically one or two orders of magnitude greater than f/2. 

Of the three functions contributing to the QE, functions /i and /2 have 
prefactors e^^ and e~^^ respectively. Neglecting the contribution from function 
/2 and terms in e~^^ in function fi gives the following expression for the QE near 
the band edge 

Cancelling hj a + k and assuming a <^ k in the resulting denominator gives 
the following result, expressed in terms of the diffusion length: 

QE ~ L^ae""" (A.21) 

The QE in this regime is independent of the surface boundary condition and 
hence of the diffusivity. We are left with only one fitting parameter, which is the 
diffusion length. 

An identical result is reached if we assume that all carriers generated within 
a diffusion length of the depletion region are collected. The quantum efficiency 
in this picture is determined by the number of carriers generated in a region of 
thickness L„ adjacent to the depletion layer. For low absorption and a constant 
generation rate, this is simply the product of the generation rate and the diffusion 
length, which is the result we obtain in equation IA.21I 

Figure IA.3I illustrates the agreement between this approximation and the ex- 
act result for our standard 30% Al p layer of thickness 0.15/im for a standard 
diffusion length of 0.08/im and a high surface recombination velocity of 500ms~^. 
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Figure A. 3: Comparison of exact and approximate solutions to the QE for low absorption 
and low diffusion length L„ = 0.07/im for x^jp = 0.148/im 

The approximation is remarkably accurate for wavelength within a few tens of 
nanometers of the band-edge. 

For samples obeying this limit, modelling can indicate values of the diffusion 
length from the magnitude of the QE near the band-edge, since the absorption 
coefficient is well known. The QE near the surface can then be modelled in terms 
of the surface recombination velocity. This is a useful diagnostic tool, which has 
already been proposed by Hovel [69] and used by Ludowise [27] . 

At higher energies the generation rate becomes strongly position dependent 
and the expression overestimates the minority carrier concentration at the edge 
of the depletion region, leading to overestimated values for the QE. 

If the diffusion length is not small compared to z, minority carriers can dif- 
fuse to the surface and recombine. The QE across the whole wavelength range 
becomes sensitive to the surface recombination. 
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A. 2. 3 High Absorption 

Figures lA.ll and IA.2I show that function tends to zero at short wavelengths. 
This is because the term e~"^ becomes smalL To a good approximation, the QE 
in this wavelength range is approximately given by 



~ /i + /2 (A.22) 
Neglecting terms in e~"^ in functions /i and /2 then yields 

QE . .,,°''rf', (A.23) 

(/e^ — a ){<-> smh[kz) — k cosh[kz)) 

Low Diffusion Length 

For low diffusion length, sinh(fc2;) ~ cosh.{kz) ~ Equation IA.23I becomes 

The short wavelength QE is sensitive to both the surface recombination velocity 
and the diffusion length 

Since, however, section IA.2.21 has shown that a low Ln can be estimated 
by modelling experimental data in the long wavelength, low absorption regime, 
parameter S can be estimated seperately by modelling the short wavelength 
response. 

High Diffusion Length 

If the diffusion length is large compared with z, smh.{kz) ~ kz and cosh(A;2;) ~ 1. 
Furthermore, k is small, and k. Equation IA.23I takes the simple form 

QE ^ (A.25) 
a^[oz — 1) 

This expression is accurate to within about 5% of the QE for a range of about 
lOnm below 400nm, for typical structures. 
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The QE in this case largely determined by the surface recombination param- 
eter. This means it is sensitive to a combination of recombination velocity and 
minority carrier diffusivity. Recombination losses in the bulk become negligible 
compared to surface losses, and the QE reaches a maximum which is limited by 
surface recombination. 

A. 3 Homogeneous Cells with Electric Fields 

We now examine the hypothetical case of a cell with a non-zero effective field 
in the doped regions but with constant transport parameters. Models of this 
type are frequently used in theoretical investigations of compositionally varying 
structures (see section [3.1.ip . 

The parameters of equation for this type of structure are are b = 2Sf and 
c = —1/Ln, where definitions and solutions are as given in section 13. 7. 2^ but are 
repeated here for convenience and take the form 

ki = +^S] + L-2 - Sf (A.26) 

k2 = -^S} + L-^-Sf (A.27) 

The two functions making up the restricted solution to the minority carrier equa- 
tion were shown to be 

Ui = e'^i" (A.28) 

U2 = e^^^^ (A.29) 

Although the boundary condition at the edge of the depletion region remains 
unchanged, a field term appears in the surface recombination boundary condition, 
as explained in 13.61 The surface recombination parameter S of equation IA.6I is 
replaced with 

S = 2£f- SjDn (A.30) 
In a well designed cell, the field is negative, causing carriers to drift towards the 
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space charge region. This study deals with fields E ~ —10^Vm~^ , such that 
Sf ~ — 10^m~^ at room temperatures. This corresponds to a moderate grade 
from 40% Al at the front surface of the layer to 30% at the top interface of the 
depletion layer. 

The presence of a field removes some of the symmetry we found in the solutions 
for an ungraded layer. Function U2 is no longer the inverse of f/i. It becomes 
necessary to consider the relative magnitudes of Sf, Ln and Sn/Dn- 

The presence of a field term in the surface boundary condition increases the 
magnitude of the surface boundary parameter S. This indicates a higher gradient 
in the minority carrier concentration n'(0) (subject to the magnitude of ^(0)). 
This seemingly contradictory increase in the value of S is due to the fact that 
minority carriers are removed away from the surface by the effective field and 
swept towards the depletion layer. The minority carrier profile as a whole is 
shifted in a positive direction, increasing the gradient and hence the QE at the 
edge of the depletion layer. This is illustrated in chapter [51 

The solution for the QE in this case is of a similar form to that for an ungraded 
cell, but with the new expressions 13.421 to IA.30I substituted. Accordingly, we 
rewrite equations IA.14I to IA.16I in terms of three new functions: 



ei = 



(a + ki){a + k2) 



a 



k^e^i^{e-'''{S + fca) + e^^^{a - S)) 



(A.31) 



62 = 



(a + ki){a + /C2) 



a 



fc2e^2^(e-"^(5 + ki) + e'=i^(a - S)) 



(A.32) 



az 



(a + ki){a + ^2) 



(A.33) 



such that the fractional quantum efficiency is 



QE{\) = 61 + 62 + 63 



(A.34) 
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A. 3.1 Surface Recombination Effect with Field 

In this section we consider analogue to the balance between bulk and surface 
recombination which was discussed in section IA.2.11 The discussion of this case 
is more restricted than for an ungraded cell because the spatial variation of 
materials parameters forbids many of the simplifications we were able to make 
in the previous section. We restrict the discussion principally to the effect of the 
surface recombination. 

In section lA. 2. II we saw that the relative magnitudes of k and S determine the 
dominant recombination effect. Equations IA.31I and IA.32] show that the QE in 
this case is insensitive to the surface recombination velocity if ki+ S and k2 + S 
are both weakly dependent on Sn/Dn 

If the diffusion length is very small such that £f ^ L^"^, we return to the 
situation described in section IA.2.11 A calculation shows that this limit applies 
only for diffusion lengths below O.Olyum for the field strengths mentioned above. 

More significantly, functions eiand 62 will be insensitive to the surface recom- 
bination velocity if S is dominated by 2£f, that is, if 

I |» A. (A.35) 

This condition is less stringent than equation IA.18I of section IA.2.11 and is 
fulfilled even the moderate values of Sf of ~ — 10''m^^. 

A. 3. 2 QE Near the Band Edge for Small Diffusion Length 
with Field 

Because £f is negative, we see that the positive ki is greater than k whereas the 
magnitude of /c2 is smaller than k. The condition of section lA. 2. II is more readily 
fulfilled. That is, 

6^=1" > e*=^" (A.36) 

This remains true even for very large diffusion lengths, for which A;i ~ + | £j | 
and k2 — — \ Sf \. 
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We follow a similar argument to that given in section IA.2.21 Neglecting terms 
in e^^^ and assuming a <^| ^2 | gives 

QE ~ , (A.37) 

This expression is in fact a restatement of equation IA.21I in the presence of a 
field. 



A. 4 Conclusions 

A. 4.1 Ungraded Solar Cells 

In section IA.2I we established that the analytical expression for the QE of a p 
layer can be reduced to simple analytical expressions in a number of limiting 
cases. 

For high absorption coefficient at short wavelengths, the QE is in general 
sensitive to both surface recombination and diffusion length, and tends to a 
constant value. 

For low absorption whithin a few tens of nm of the bandedge, the QE is 
independent of the surface recombination for the structures we are studying, 
because the diffusion length is smaller than the width of the p layer. The only 
free parameter in this regime is the diffusion length. 

Values of the diffusion length can therefore be estimated from the QE at long 
wavelengths. The short wavelength response is then modelled independently in 
terms of the recombination velocity. 

Cases where very large values of S are required indicate a poor surface and 
the existance of a dead layer and a zero minority carrier concentration at the 
front of the p layer. We have seen however that the QE in the low absorption 
regime remains principally determined by the diffusion length. 
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A.4.2 Solar Cells With Fields but Constant Materials Pa- 
rameters 

A first approximation to the QE for a compositionaiiy graded ceii is tlie exampie 
of an ungraded cell where an effective field is present in the doped regions. Fewer 
conclusions can drawn from this case than for an ungraded cell because the spatial 
variation of materials parameters is not negligible. 

The main conslusion wc reach is that even moderate grades generate fields of 
sufficient magnitude to reduce the effect of the recombination velocity to negli- 
gible levels. 

This model however overestimates the QE for a graded structure for two 
reasons. The constant absorption coefficient overestimates the amount of light 
absorbed, particularly at long wavelengths. Secondly, comparison between nu- 
merical and analytical calculations shows that the gradient in the transport pa- 
rameters is important. The gradient in Dn, and hence the minority carrier mo- 
bility, contributes a field-like term which is opposite to the band gap gradient 
field. 

The QE in a real cell therefore remains insensitive to the magnitude of the 
recombination velocity for realistic values of the latter, but is significantly lower 
than the predictions of the simplified analytical model. 
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